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EDITORIAL 


With this issue we begin our second year of operation of PHYSICAL 
REVIEW LETTERS. The journal has received general approval and seems to 
be meeting a real need. The subscriptions, although slower than expected at 
the start, have increased rapidiy andthe financial situation is now reasonably 
healthy. 

Unfortunately, the number of contributions has also increased. As we 
have pointed out before, the journal would lose its usefulness if it became so 
voluminous that its contents could not be scannedeasily. Moreover, our facil- 
ities are limited and we cannot handle more than 15 to 20 Letters per issue. 
(The size of the current issue, which is 50% larger, cannot be maintained on 
a continuous basis.) We therefore must apply quite stringent criteria of im- 
portance and urgency, andrequire the full cooperation of contributors and ref- 
erees in restricting the Letters to those whose rapid publication can be of 


greatest value. 


There are two innovations on the cover of this issue of PHYSICAL 
REVIEW LETTERS, both designed to facilitate storage and use. The three 
heavy bars at the top designate the volume number, and the thin line near the 
bottom designates the issue number within the volume. Also, the back cover 
near the fold carries the volume number, page numbers, 
and year. Subscribers who do not bind their completed vol- 
umes may find it useful to mark the covers of volumes 1 and 


2 correspondingly. (The spacing of the issue number lines 








is illustrated for an issue No. 7.) 


S. A. Goudsmit and S. Pasternack 























" u Pe Ge BOT we Oo Oe Bette DH Shh Hw we ree De 

















PHYSICAL REVIEW 





LETTERS 





eee 





VOLUME 3 





JULY 1, 1959 _ 


NuMBER 1 





FIELD MIXING AND ASSOCIATED NEUTRON PRODUCTION IN A PLASMA* 


A. C. Kolb and C. B. Dobbie 
United States Naval Research Laboratory, Washington, D. C. 


and 


H. R. Griem 
University of Maryland, College Park, Maryland, and 
United States Naval Research Laboratory, Washington, D. C. 
(Received June 4, 1959) 


A series of experiments has been performed 
to determine the magnitude of the internal elec- 
tric fields in a preheated deuterium plasma with 
a trapped internal axial magnetic field H;, which 
is compressed by a large externally generated 
axial field H,. It has been shown that Ep fields 
of several hundred volts/cm can be readily pro- 
duced by the rapidly changing internal magnetic 
fields and result in the production of neutrons 
from D-D reactions. 

The electric field strength in the plasma can 
be controlled by varying separately (dH,/dt), Hj, 
and the relative polarity of the internal and ex- 
ternal fields. The highest Ep fields are observed 
when the two magnetic fields, which are gener- 
ated by separate condenser banks (the “preheater” 
and “main” banks) have opposite polarity. Dur- 
ing the early stages of the compression, on the 
second half-cycle of the preheater and the first 
half-cycle of the main bank, the trapped mag- 
netic field first increases in magnitude as the 
plasma radius diminishes and then decreases 
rapidly, reversing its direction as the two fields 
interpenetrate [Figs. 2(C), (E)]. Shortly after 
the internal magnetic field passes through zero, 
with dH ;/dt~ 10" gauss/sec, a burst of neutrons 
is observed. For a trapped field of ~+10000 
gauss, which is reversed to -10000 gauss in 
~0.2 usec, e.g., see Fig. 2(E), a kev deuteron 


can double its energy in about one Larmor period. 
Particles with energies sufficient to cause a D-D 
reaction can easily be produced in times of 

~1 psec by the induced electric field which trans- 
fers energy from the reversed trapped field to 
the plasma. The particles are confined by the 
external field when the internal field passes 
through zero since the external field from the 
second half-cycle of the preheater and the first 
half-cycle of the main bank are in the same di- 
rection and do not cancel one another. 

The cancellation of a trapped reversed axial 
magnetic field from the first half-cycle by the 
penetration of an external field generated during 
the second half-cycle of a fast condenser dis- 
charge was first suggested by Furth’»* as the 
mechanism responsible for neutron production in 
other experiments.*»‘ (In these cases there was 
no trapped reversed field during the first half- 
cycle since there was no preheater field as in the 
present experiment.) This suggestion was based 
on earlier observations of Colgate® with rapidly 
oscillating axial fields (the “collapse” experi - 
ment), on field measurements of Anderson et al.” 
with inverse pinches, and also on experiments 
performed with the triaxial pinch.* The latter is 
a tubular pinch compressed between magnetic 
fields of opposite sign and neutrons are emitted 
near a current maximum with an isotropic en- 
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when there are rapidly cancelling reverse fields. 
Our experiment with axial confining fields in 
which the internal trapped field can be controlled 
lends support to the hypothesis of Furth for the 
trapped field origin of “second half-cycle neu- 
trons” and does not appear to be consistent with 
the suggestion’ that the energetic deuterons are 
produced by simple Joule and shock preheating 
with subsequent adiabatic compression by an ex- 
ternal field. The cancellation of the trapped field 
apparently plays an essential role. 

Also in this experiment the trapped field from 
the first half-cycle of the main condenser bank 
is observed to be trapped at the beginning of the 
second half-cycle [e.g., Fig. 2(H)] and then com- 
pressed by the rising external reversed field. 
This trapped field appears to be responsible for 
the excitation of large amplitude radial plasma 
oscillations observed previously® as well as 
providing a plausible mechanism for the neutron 
generation reported elsewhere.*’* 

The experimental parameters employed for 
most of this investigation were: voltage 15 kv, 
capacitance 130-400 yf for the main bank and 


Hi (4000 GAUSS/DIVISION) V (10 KV/ DIVISION) 





FIG. 1. Internal magnetic field Hj and voltage V 
across the coil for the preheater discharge: A, B one 
atmos air (no gas current); C, D 100 microns D, (no 
rf pre-excitation), C shows the trapped field on the 
third half-cycle; E, F 100 microns D, (rf pre-excita- 
tion), E shows the trapped field on the second half- 
cycle. The probe is positioned at the central plane of 
the coil. The time scale is 1 psec/division. The dots 
serve to improve the clarity in the oscillograms. 
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ergy distribution at a time of enhanced resistivity 





20 kv, 1 uf for the preheater discharge; maxi- 
mum current, 1.1x10*° amp; maximum field, 
100 000 gauss; total inductance 0.05 yuh; coil 
inductance 0.04 wh; coil length 10 cm; coil diam- 
eter 5.7 cm i.d.; Pyrex tube 5.0 cmi.d. The 
electrical system is described more completely 
elsewhere.* No magnetic mirrors were used so 
that deuterons moving axially with energies 
above one kev could escape from the confining 
field in a few tenths of a microsecond. There- 
fore, the neutron generation is probably caused 
principally by a highly nonisotropic distribution 
of deuterons moving in a plane normal to the 
magnetic axis. The observation of neutrons with 
no mirror fields suggests that the presence of 
magnetic mirrors is not the dominant factor in 
the initial heating since deuterons are apparently 
raised to high energies in a short period of time 
by an induced Eg field in the plasma. 

Typical observations are shown in Figs. 1 and 
2. The magnetic field from the second half- 
cycle of the preheater is trapped on its third 
half-cycle if there is no rf (1 kw, 27 Mc/sec) 


660 KEV CS-7 








FIG. 2. Internal magnetic field Hj with internal and 
external fields in opposite directions (+-) C,E and in 
the same direction (++) A, G, H showing the compres- 
sion and rapid reversal of the trapped preheater field 
(C, E) and the trapped field from the first half-cycle 
of the main bank (G, H). The voltage trace F corre- 
sponds to the field E; the modulation is due to the 
coupled preheater and main banks. The time of neu- 
tron emission is shown in D (corresponding to C) and 
the scintillations due to cesium y-rays are shown in B 
for comparison. The scale for H; is 15000 gauss/divi- 
sion and the time scale is 1 ysec/division. The probe 
was positioned on the magnetic axis at the end of the 
coil for A and C and at the central plane for E, G, and 
H, 
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pre-excitation [Fig. 1(C)]. With rf pre-excitation 
the field from the first half-cycle is trapped and 
on the second half-cycle the internal field [Fig. 
1(E)] has a direction opposite to the field ex- 
ternal to the plasma. The magnetic probe signal 
with the tube filled with air at atmospheric pres- 
sure (no ionization) is shown in Fig. 1(A) for 
comparison. The difference in the voltages 
across the coil with air (1 atmos) [Fig. 1(B)] and 
deuterium (100 microns) [Figs. 1(D), (F)] shows 
a marked effect when a plasma is created and 
the load inductance is reduced accordingly. 

The neutrons were identified with a BF, pro- 
portional counter surrounded by paraffin and the 
time of emission was determined with a lead- 
shielded (1 cm) plastic scintillator. Additional 
shielding is provided by the 7.5-cm wall of the 
steel coil. Figure 2(B) shows the pulse shape 
for individual scintillations obtained using 0.66- 
Mev y-rays from a cesium source, which have 
about the same average pulse height as 2.5-Mev 
neutrons. The decay time for a single scintilla- 
tion is determined by the amplifier. Neutrons 
are observed [Fig. 2(D)] if the main bank of con- 
densers is discharged during the second half- 
cycle of the preheater when the trapped field has 
the correct reverse polarity.° The neutron emis- 
sion lasts for ~0.2-2 usec which, for densities 
of ~10** - 10" cm=*, corresponds to the expected 
scattering time for particles out the ends. If the 
internal magnetic field is in the same direction 
as the main external field [Figs. 2(A), (G)] it 
rises slowly in the plasma as compared to the 
case when the trapped field is in the opposite di- 
rection [Figs. 2(C), (E)]. After about 2 psec the 
plasma radius from streak camera observations 
is comparable to the probe radius, and the field 
measured at later times is essentially the vac- 
uum field. 

The experiments have been repeated using a 
somewhat slower preheater discharge with the 


following parameters: 25-30 kv, 2 uf, 330 kc/sec. 


The highest neutron yield was observed if the 
main bank was fired near the voltage maximum 
at the beginning of the second half-cycle of the 
preheater. Hard x-rays (several hundred kev) 
were also recorded with the scintillation counter 
when hydrogen was substituted for deuterium. 
The x-ray emission occurred during four half- 
cycles of the main bank when there were large 
voltages across the coil. The neutrons appeared 
at lower coil voltages. The polarity of the trap- 


ped preheater field did not seem to be a critical 
factor in the x-ray production. 

Discussions with Dr. Harold Grad and Dr. 
Harold Furth have been of considerable value in 
the interpretation of these experiments which 
were carried out with the technical assistance of 
Mr. L. J. Melhart and Mr. T. H. DeRieux. 





* Jointly supported by the U. S. Atomic Energy Com- 
mission and the Office of Naval Research. 

1H. P. Furth (private communication, Geneva Con- 
ference, 1958). 

2 Anderson, Furth, Stone, and Wright, Phys. Fluids 
1, 489 (1959). 

3Elmore, Little, and Quinn, Phys. Rev. Letters 1, 
32 (1958); Proceedings of the Second United Nations In- 
ternational Conference on the Peaceful Uses of Atomic 
Energy, Geneva, 1958 (United Nations, Geneva, 

1958), Vol. 32, p. 337, P/356. 

‘kK. W. Allen and B. Niblett, Conference on Con- 
trolled Thermonuclear Reactions, Berkeley, 1959 (un- 
published); W. Millar (unpublished) communicated by 
P. C. Thonemann; H. Hurwitz (private communication). 
The hard radiation reported by 8. M. Osovets et al., 
Proceedings of the Second United Nations International 
Conference on the Peaceful Uses of Atomic Energy, 
Geneva, 1958 (United Nations, Geneva, 1958), Vol. 32, 
p. 311, P/2225 (see Fig. 24), may originate from the 
same mechanism. A trapped field from an auxiliary 
coil was rapidly compressed and cancelled during the 
collapse of a plasma loop. 

5S. A. Colgate and R. E. Wright, Proceedings of the 
Second United Nations International Conference on the 
Peaceful Uses of Atomic Energy, Geneva, 1958 (United 
Nations, Geneva, 1958), Vol. 32, p. 145, P/368. 

® Anderson, Baker, Ise, Kunkel, Pyle, and Stone, 
Bull. Am. Phys. Soc. 4, 119 (1955). 

TBoyer, Little, Quinn, Sawyer, and Stratton, Phys. 
Rev. Letters 2, 279 (1959). 

*A. C. Kolb, Proceedings of the Second United Na- 
tions International Conference on the Peaceful Uses of 
Atomic Energy, Geneva, 1958 (United Nations, Geneva, 
1958), Vol. 31, p. 328, P/345. 

*In earlier experiments [Griem, Kolb, and Faust, 
Phys. Rev. Letters 2, 281 (1959)] the existence of 
electron temperatures of ~ 8 x 10K were inferred 
from an analysis of the bremsstrahlung radiation in the 
soft x-ray region emitted from a dense deuterium 
plasma confined by an externally generated magnetic 
field. No neutrons were observed with the scintillation 
counter which has a sensitivity of 1.5 x 10-4 count/neu- 
tron. This negative result is consistent with the meas- 
ured density, temperature, and volume of the plasma. 
No trapped reverse field was present in this case and 
the induced E, fields were too small (a few volts/cm) 
to generate neutrons by the mechanism discussed here. 
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NEW APPROACH IN THE THEORY OF SATELLITE ORBITS* 
John P. Vinti 


National Bureau of Standards, Washington, D. C. 
(Received April 27, 1959) 


If O is a center for the earth and if Oz points 
along the polar axis and Ox towards the vernal 
equinox, the oblate spheroidal coordinates £é, 7, 
@, defined by 


x =c[(#*+1)(1 - 1°)” cosd, (0<t<@) 
y =c[(é+1)(1 - 7)” sing, (-1<n<1) 
2 =cén, 
turn out to be a useful system of coordinates for 
investigating the motion of an earth satellite. 
In this coordinate system any axially symmetric 
potential V(é, 7) that results in separability of the 


Hamilton-Jacobi equation must have the general 
form 


V(E, m) = (e+ 7) “*[7(E) + gn) }. 


The most general solution of Laplace’s equation 
compatible with such a form is 


V =C +b, Re(é+in)~* + b, Im(é+in)™ 


+69 fastan'tenin( +), 





1-9 
where the logarithmic term has a singularity 
everywhere on the z axis. We therefore place 
b, = 0 and we also place C =0 to make V vanish at 
infinity. 

Expansion of (¢+in)~* in spherical harmonics 
leads naturally to the following choices for the 
adjustable constants 5, and c: 


— ba. -4 
bc =-GM, c U, 1)/M, 


where G is the gravitational constant, M is the 
mass of the earth, and Ip and J; are, respect- 
ively, the polar and transverse moments of in- 
ertia of the earth. 

The above V then everywhere gives the r™* 
term and the second harmonic exactly. For the 
fourth harmonic it yields an amplitude which has 
about half of Jeffreys’ value’ or about 86% of the 
value obtained by King-Hele and Merson,” from 
analysis of data on satellite orbits. Since these 
authors indicate a standard deviation of 16%, the 
resulting value for the fourth harmonic agrees 


with their value within the limits of observational 
error. Even the sixth harmonic, of amplitude 
1.3x10°*, falls within the (admittedly wide) am- 
plitude limits, (0.1+1.5)x10~*, given by King- 
Hele and Merson. 

[The term b, Im(é+in)~* leads to odd harmonics, 
with amplitudes all proportional to the first 
harmonic amplitude 5/R, where R is the equa- 
torial radius and 6 is the distance OC from the 
coordinate center O to the center of mass C. 
Thus if O is taken as coincident with C, all the 
odd harmonics drop out, so that the potential 
V(é,n) here derived always has symmetry with 
respect to a plane through C, perpendicular to 
the polar axis. It remains an open question 
whether use of a displaced center O, noncoinci- 
dent with C, which would give rise to first, third, 
..., harmonics with amplitudes 5/R, 0.001095/R, 
..-, could usefully represent equatorial asym- 
metry, as found by O’Keefe, Eckels, and 
Squires.*] 

The net result of the present work is a reduc- 
tion of the problem of satellite motion to quad- 
ratures, for a potential whose even harmonics 
are exact through the second, and whose fourth 
harmonic is comparable with empirical values. 
The kinetic equations of motion have been set up 
formally, in terms of certain integrals, which 
are now being evaluated. The effects of drag, of 
axial asymmetry, and of higher odd harmonics 
remain to be handled by perturbation theory. 

This paper had its inspiration in the earlier 
work of Sterne* and of Garfinkel,® who also re- 
duced the problem to quadratures, for potentials 
with approximate second harmonics. 





“This work was supported by the U. S. Air Force, 
through the Office of Scientific Research of the Air Re- 
search and Development Command. 

‘Harold Jeffreys, The Earth (Cambridge University 
Press, Cambridge, 1952), second edition, Chap. IV. 

?D. G. King-Hele and R. H. Merson, Nature 183, 
881 (1959). 

30’Keefe, Eckels, and Squires, Science 129, 565 
(1959). 

‘T. E. Sterne, Astron. J. 62, 96 (1957). 

*B. Garfinkel, Astron. J. 63, 88 (1958). 
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PROTONS IN THE EARTH’S MAGNETIC FIELD* 


Stanley C. Freden and R. Stephen White 
Lawrence Radiation Laboratory, University of California, Livermore, California 
(Received May 28, 1959) 


The identity, flux, and energy distribution of 
particles trapped in the earth’s magnetic field 
have been determined for those particles which 
penetrate more than 6 g/cm? of material. 

A small stack of nuclear emulsions has been 
exposed to the particles in the lower Van Allen 
radiation belt." The stack consisted of 10 sheets 
of Ilford K.5 emulsion, 1% in. x1} in. x600 y, 
enclosed in a water-tight stainless steel box 
which was mounted on the instrument board 13 in. 
from the side wall of the nose cone of a Thor- 
Able ballistic missile. The normal to this wall 
was parallel to the plane of the emulsions. 

The missile was launched from Cape Canaveral 
on April 7. It reached a maximum altitude of 
1230 km and spent about 15 minutes above 1000 
km between the latitudes of 20°N and 3°N. The 
nose cone was recovered in the South Atlantic 
and returned for analysis. The emulsions were 
in excellent condition. 

The particles that penetrated 5 g/cm? of nose 
cone wall, 0.032 in. of stainless steel emulsion 
box, and a minimum of 1 mm of nuclear emulsion 
have been identified and their numbers and en- 
ergies measured; these correspond to protons 
and electrons with energies greater than 75 Mev 
and 12 Mev, respectively. The separation of 
protons from electrons was done unambiguously 
for proton energies up to 700 Mev by measuring 
the ionization of their tracks, only. The back- 
ground of 7 mesons and other particles that arose 
from trapped proton or cosmic-ray interactions 
in the walls of the nose cone was negligible. No 
attempt was made to determine the angular dis- 
tribution of the charged particles in the belt be- 
cause the nose cone of the missile was not orien- 
ted in space. Therefore, in all the results which 
follow, an isotropic distribution has been as- 
sumed. 

All proton energies were measured using the 
g* method of ionization.? Normalization was ob- 
tained for g*’s from 5 to 7 using the tracks of 
protons which had residual ranges of 8 to 14 mm 
in the emulsions. The energies of these stopping 
protons were determined from their ranges. The 
scanning for tracks was done 5 mm in from the 
edge of the emulsion that was nearest the nose 
cone wall where accurate g* measurements could 


be taken. In addition, scanning was done 1 mm > 
from the edge to pick up low-energy protons 
which did not penetrate to 5 mm. All particles 
whose vertical angles were within +15° of the 
emulsion plane and whose horizontal angles, in 
the plane of the emulsion, were within + 15° of 
the normal to the wall of the nose cone were ac- 
cepted for measurement. The estimated absorber 
thickness in the backward direction was used to 
correct the proton spectrum for the number of 
protons which crossed the scan lines in the wrong 
direction. The calculated correction varied con- 
tinuously from 10% at 75 Mev to 25% at 700 Mev. 
This was checked experimentally at 100 Mev. 

Emulsions which underwent the same history 
as those exposed to the Van Allen belt except for 
being aboard the missile were scanned for back- 
ground tracks. The proton background was small. 
After the background of tracks of plateau ioniza- 
tion was subtracted out, no electron signal re- 
mained. An upper limit of 1% is placed on the 
ratio of the number of electrons to protons that 
can penetrate 6 g/cm’ of material. 

The resulting value, N, the absolute number of 
protons Mev™ sterad™' sec™* cm~?, outside the 
nose cone at an altitude of about 1200 km is plot- 
ted versus the proton energy T in Fig. 1 ona 
log-log plot. The data approximate a straight 
line from the minimum detectable energy of 75 
Mev to the maximum observed energy of 700 Mev. 
The data are represented by 


N=N,T™", (1) 
where a least-squares fit to the data gives n=1.84 


+ 0.08 and N, = (2.139) x10° protons Mev™ sterad™* 
sec" cm™*. The data also fit an exponential 


N=N,e72/To , (2) 


with T,=120+5 Mev, and N,=0.84+0.07 proton 
Mev™ sterad™ sec”! cm~?. 
At least two theories have received consider- 


able attention for the origin of the trapped par- 
ticles in the earth’s magnetic field: (1) Solar in- 
jection.’ Since a solar origin would give protons 
of a few Mev, at most, it seems unable to ac- 
count for the high-energy protons found in this 
experiment. (2) Earth neutron albedo. In this 
theory neutrons from cosmic-ray interactions in 
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PROTON KINETIC ENERGY (Mev) 


FIG. 1. The absolute energy spectrum of protons 
above 75 Mev in the lower Van Allen radiation belt at 
an altitude of 1200 km. The line is a least-squares 
fit of the equation N=N,7~™ to the data and gives 
N, =(2. 17}: % = 10° protons Mev sterad~! sec”! cm=* 
and n=1.8420.08. 


the upper atmosphere travel to the region of 
trapping and there decay into protons and elec- 
trons.*~* The earth neutron albedo can explain 
the data since high-energy neutrons are formed 
in cosmic-ray stars in the atmosphere in suffi- 
cient abundance to give the number of protons 
observed here. 

Singer’ predicted a velocity distribution that 
was proportional to v-°**. This distribution is 
transformed to give an energy distribution pro- 


portional to T~°*. Hess* made use of the neutron 


energy distribution that he measured high in the 
atmosphere to predict the neutron flux in the 
trapping region. From that he obtained a proton 
energy distribution proportional to T~*** for 

T, 2100 Mev (see the accompanying Letter). 
Neither of these distributions is as steep as the 
experimental data. If only the data for proton 


10 


energies <300 Mev are used, m=1.56+0.12 and 
N, = (5.274-9) x1 protons Mev™ sterad™* sec™ 

cm”. This difference in m can be explained by 
the larger orbits and the subsequently greater 

loss of the high-energy protons. 

The omnidirectional proton flux, J,,, i.e., the 
number of protons above 75 Mev which pass 
through a sphere of 1 cm? cross section per sec, 
is found by multiplying Eq. (1) by 4m and inte- 
grating over all proton energies greater than 75 
Mev. The resulting value of J,, is 800 +200 pro- 
tons cm™ sec™ If the power law with n=1.56 
is used to integrate down to a proton energy of 
10 Mev, the resulting omnidirectional flux is 
3 x10° protons cm™ sec” which is not in agree- 
ment with the value of 100 to 1000 that was ob- 
tained from previous data.’ It appears that a 
still smaller value of m should be used at low en- 


ergies. Better agreement with this previous result 


is obtained by using the exponential fit which 
gives an omnidirectional flux of 1100 protons 
cm~* sec™* 

The radiation level at the position of the emul- 
sion stack inside the nose cone due to electrons 
and electromagnetic radiation is less than 50 
mr/hr. Using the exponential extrapolation of 
the proton data, the radiation level in space at 
1200 km due to protons is estimated to be 1 r/hr. 

It is a pleasure to acknowledge the assistance 
and cooperation of personnel in the Space Tech- 
nology Laboratories both at Los Angeles and at 
Cape Canaveral, General Electric— Missile and 
Space Vehicle Department, and Avco— Research 
and Advanced Development Division. We thank 
John Lindner for coordinating the experiment 
with the missile requirements. We are grateful 
to Al Oliver for assembly of the stacks. We es- 
pecially thank Miss Patricia Banks, Mrs. Irene 
Brown, Miss Claudia Crane, Mrs. Theodora 
Hilliard, Mrs. Beverly Lagiss, Mrs. Bobby Saw- 
yer, Mrs. Nancy Stratton, Mrs. Sophie Ward, 
and Mrs. Elizabeth Wilson for locating the proton 
tracks and for taking measurements on them. 
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Mcliwain, and Ludwig, J. Geophys. Research 64, 271 
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VAN ALLEN BELT PROTONS FROM COSMIC-RAY NEUTRON LEAKAGE 


Wilmot N. Hess 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received June 1, 1959) 


One of the sources of particles for the Van 
Allen radiation belt is the decay of cosmic-ray 
neutrons leaking out of the atmosphere of the 
earth. Recently the cosmic-ray neutron energy 
spectrum has been measured,’ and from this, by 
use of a multigroup neutron-diffusion calculation,’ 
the neutron leakage has been calculated (Fig. 1). 
Some of the neutrons leaking out of the atmos- 
phere decay in the earth’s magnetic field and are 
trapped. From this leakage we can calculate the 
equilibrium proton energy spectrum in the inner 
Van Allen belt. 

The calculated neutron leakage is given by 
o(E) = 0.8E~? neutron/cm? Mev sec in the region 
from 10 Mev to 1 Bev. The decay density of 
neutrons, dn/dV, is given approximately by 


dn 1 R 
ave baer (| 


‘R 2 
x o(E \(=2) neutron decays/cm® Mev sec, 
where L =neutron mean life, R, = radius of earth, 

v=neutron velocity, R =distance from earth’s 
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FIG. 1. The energy spectrum of the cosmic-ray 
neutrons that leak out of the atmosphere of the earth. 


center, and y is the time-dilation factor. Pro- 
tons resulting from neutron decays have very 
nearly the energy of the parent neutron. This 
gives us, for the proton source, 


S(E) =k, [o(E)/By]. 


Whereas the loss mechanism for electrons in the 
Van Allen belt is multiple small-angle Coulomb 
scattering, the loss mechanism for protons is 
slowing down by collisions with bound electrons.* 

In order to calculate the equilibrium proton en- 
ergy spectrum we must consider the continuity 
equation for the slowing-down process. Follow- 
ing Singer* we can write 


dN(E) 


d dE 


for the steady state, where N(E) is the equilibri- 
um proton energy spectrum. This continuity 
equation follows the flow of protons along an en- 
ergy axis. Now we can write dE/dt = (dE/dx) 

x (dx /dt) = (dE /dx)Bc; then, by substituting into 
the continuity equation, we get 


k,g(E) a ‘ dE 
By ale "ls “| 


where we have taken N(E) =k,E~”. Approxima- 
ting 8, By, and dE/dx as functions of E, we can 
write* 


B=0.0885 E344 
By =0.0393 E°-54, 
dE /dx =1.16 E~°®™ Mev/cm of NTP air. 


These are all accurate to 5% or less in the en- 
ergy region 80 Mev to 700 Mev. The continuity 
equation now becomes 


bE -2.545 flee] =k, -M-2-242 
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This gives, for the desired equilibrium proton 
energy spectrum, 


NE) =k,E™" =k, E “1.30 


This is shown plotted in Fig. 2, extended down 

to 1 Mev by refitting 8, By, and dE/dx and re- 
peating the process above. The peak in the 
spectrum is at about 1 Mev. There are two 
reasons for this; first, this is about the position 
of the peak in the neutron leakage spectrum (see 
Fig. 1), and secondly, below this energy, charge- 
exchange reactions on neutral hydrogen remove 
the particles quite rapidly.* One can see that 

the equilibrium spectrum is flatter than the 
source spectrum. This says that the proton ef- 
fective lifetime increases with the proton energy, 
as might be expected from the nature of the loss 
mechanism. We can find the lifetime here by 
using the “leaking bucket” equation as Van Allen® 
has done: 


Input = output = contents/r. 


The mean life 7 in this problem is given by 


7 =contents/input =k,[T *-**/T-?-**]=k,T***. 


Freden and White, in the accompanying Letter, 
report measurement of the proton spectrum of 
the inner Van Allen belt.” The shape of our spec- 
trum agrees well with theirs from 90 Mev to 
about 200 Mev. They find considerably fewer 
protons above 200 Mev than are predicted by the 
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FIG. 2. The equilibrium energy spectrum of protons 
in the Van Allen radiation belt. A, the spectrum cal- 
culated in this report. B, the measured spectrum of 
Freden and White, fitted by an exponential law N(E) 
=k exp(-E/125) and normalized at 100 Mev. 
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present analysis. There are at least three pos- 
sible causes of this difference. (a) The neutron- 
leakage angular distribution is a function of en- 
ergy.” This has been neglected, and although it 
is not obvious how it would, this might change 
the spectrum. (b) A 500-Mev proton at a dis- 
tance of one earth radius from the surface has a 
radius of gyration of 600 km. This is so large } 
that this particle will reach appreciably further , 
into the atmosphere of the earth at its mirror : 
point than lower energy particles, and thus its 
lifetime will be shorter. (c) This large radius ) 
of gyration means that this particle will be quite 
susceptible to nonadiabatic processes, which 
will lower its energy. For example, time fluc- 
tuations in the earth’s magnetic field® would act 
on this particle (to remove it) more readily than 
on lower energy particles. 

If we evaluate the constant in the expression 
for the proton lifetime, we get 


7=2.1x10"'° Me™ sec, 





or 


where M is the thickness of air (in cm) that must 
be traveled by the particle in its trajectory to 
equal 1 cm of air at NTP: 


7 atoms/cm at NTP , 6 x10" 
~ atoms/cm® along real path Y 





Following Christofilos’ analysis, ° we get 


Y= (No/p)(h/r.)”, 


where N,=mirror-point density, 7x7, =length of 
magnetic line from mirror point to equator, 
Y,=mirror-point radius, and h = scale height of 
atmosphere. If we take N,=10° atoms/cm (for 
about 1100 km, according to Johnson’®) and y =2 
and h=100 km, we get Y=5.7x10° atoms/cm*. 
This gives t=6X10* sec for a 100-Mev proton. 
Now if we use the neutron decay density dn/dV 
=(1/vyL) o(E), we can evaluate the absolute 
value of the proton flux, F, to be expected in the 
inner Van Allen belt by using the leaking-bucket 
formula, 1=C/t. We get 


F =Cv=Itv=(dn/dV)rtv, 
F= (1/yL)ér. 








Integrating ¢(£) from 40 Mev up, we get the 
total neutron flux @=0.019 neutron/cm? sec. 
This gives, for the proton flux, F =1 x10* pro- 
tons/cm’* sec. This agrees well with what Van | 
Allen finds for the penetrating component of the 
inner belt,** which, according to the neutron- 
leakage source model, will be all fast protons. 
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This calculation is not an exact one, but it does 
show that the lifetime and flux predicted by the 
neutron-leakage source have reasonable values. 

Detailed calculations on the neutron- leakage 
source are now under way and will be published 
later. 
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POLARIZATION OF THE Al*” NUCLEI IN RUBY* 


J. A. Cowen, W. R. Schafer, and R. D. Spence 
Department of Physics, Michigan State University, East Lansing, Michigan 
(Received May 14, 1959; revised manuscript received June 1, 1959) 


In the course of a series of studies of the re- 
laxation process of the aluminum nucleus in ruby 
we have found that it is possible to observe either 
enhanced nuclear absorption or emission signals 
by driving a microwave transition of the chro- 
mium impurity. The crystal, containing 0.1% 
chromium, was cooled to 4.2°K and was irradiated 
with microwave power at 9300 Mc/sec. The dc 
magnetic field and the frequency of the nuclear 
resonance spectrometer were adjusted so that 
one of the transitions of the chromium ion (S=$) 
and at least one of the transitions of the Al?’ nu- 
cleus (J = 3) were simultaneously observed on the 
oscilloscope. The aluminum nuclear resonance 
signal in ruby consists of five components aris- 
ing from the interaction of the quadrupole mo- 
ment of the nucleus with the crystalline electric 
field. 

With the magnetic field adjusted so that the 
chromium spin resonance occurred at the low- 
field end of the magnetic field sweep, all the 
components of the aluminum nuclear resonance 
showed an enhancement [Fig. 1(a)] if sufficient 
microwave power was applied. With the mag- 
netic field adjusted so that the chromium spin 
resonance occurred at the high-field end of the 
Sweep, all of the components of the nuclear re- 
Sonance gave strong emission signals [Fig. 1(b)] 
if sufficient microwave power was applied. The 
exact position in the magnetic field sweep at 


which the largest enhanced or emission signals 
occur depends on the orientation of the crystal 
and on the microwave transition being excited. 
The dependence of the emission signal on the 
microwave power level is shown in Fig. 2. The 
maximum enhancement observed thus far is 
about 30. The effects are strongly temperature 
dependent as indicated by the fact that they dis- 
appear rapidly as the system warms up from 
helium temperature. No effects have been ob- 
served at liquid air temperature. When the 





(a) (b) 


FIG. 1. 
resonance signal in Al?’. 


One component of the nuclear magnetic 
(a) Top: enhanced signal; 
bottom: signal with no microwave power. (b) Top: 
signal with no microwave power; bottom: emission 
signal. 
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FIG. 2. The amplitude of the inverted nuclear res- 


onance signal in arbitrary units as a function of micro- 
wave power. 


microwave power is applied to the system, the 
enhanced or emission signals build up exponen- 
tially with a characteristic time which varies 
from less than a second to a few seconds depend- 
ing on the orientation of the crystal. 

These results can be explained in terms of the 
dynamic polarization process discussed by 
Abragam and Proctor’ and others.? They con- 
sider a system in which the nuclear spin-lattice 
relaxation time is much greater than the elec- 
tronic spin-lattice relaxation time. A mutual 
spin flip may take place inducing either Am,=+1 
or 4m,=-1 transitions if the frequency at which 
the electronic system is irradiated is the sum or 
difference of the electronic and nuclear resonance 
frequencies. If H,, is the amplitude of the mo- 


dulation field and H, the dc magnetic field, then 
4m] =+1 transitions are induced when the fre- 


quency of the microwave source satisfies 
v=(y, -ypy -H,,), a) 


and Am y= transitions are induced when 


v=ly +y)H)+H, ). (2) 


The resulting change in population of the nuclear 
levels produces an enhanced absorption signal 
for case (1) and a stimulated emission signal for 
case (2). It does not appear that in the absence 
of the microwave power the spin-lattice relaxa- 
tion times of the system satisfy the condition 
imposed by Abragam and Proctor, but there is 
definite evidence that the nuclear relaxation 
time increases appreciably in the presence of 
the microwave power. The energy levels of the 
chromium ion in ruby are complicated by large 
zero-field splitting and as a result they are 
strongly angle dependent and quite impure except 
at fields higher than those used in our experiment 
(~3000 gauss). Nevertheless, all orientations 
and all transitions examined thus far have given 
at least a small amount of polarization. 

We wish to thank Dr. Chihiro Kikuchi and Dr. 
John Lambe of the University of Michigan 
Willow Run Laboratories for their interest, as- 
sistance, and encouragement. We also wish to 
thank Dr. N. Bloembergen of Harvard University 
for several helpful comments. 
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EVIDENCE FOR ANISOTROPY OF THE SUPERCONDUCTING ENERGY GAP 
FROM ULTRASONIC ATTENUATION* 


R. W. Morse, T. Olsen, Tt and J. D. Gavenda 
Physics Department, Brown University, Providence, Rhode Island 
(Received May 15, 1959) 


The Bardeen, Cooper, and Schrieffer theory of 
superconductivity has shown remarkably good 
agreement with most experiments considering 
that a single isotropic energy gap has been used.’ 
However, there are disagreements, particularly 
with respect to specific heats, as was emphasized 
by Boorse in a recent Letter.” A likely explana- 
tion of such discrepancies is that the energy gap 
in most superconductors is appreciably aniso- 
tropic. Thus the specific heat variation with T 
would be determined near the transition temper - 
ature (7.) by some average value of the gap, while 
at low temperatures the smaller values would 
dominate. 

In order to demonstrate whether or not the en- 
ergy gap depends upon the electron’s position on 
the Fermi surface, an experiment must select 
electrons from a fairly limited region of the sur- 
face. Such is the case with the ultrasonic attenu- 
ation. If the electronic mean free path (J) is large 
compared to the wavelength (A), then a sound 
wave having vector q scatters electrons from 
state k to k’, where k’=k+q. In addition, a pho- 
non is absorbed or emitted; i.e., E(k’) =E(k) +%w, 
where w is the angular frequency of the sound 
wave. These conditions of momentum and energy 
conservation are satisfied only by a particular 
group (or groups) of electrons on the energy sur- 
face. Since Aw is small, then aE =(ak)-(V,E), or 
hw =q- (V,E), where V, is the gradient operator 
ink space. Since the electronic group velocity is 
¥ =h" VE, the required electrons are those for 
which v,cosé@ =v, where v, is the velocity of 
sound and 6 is the angle between q and ¥,. Thus 
the scattered electrons are ones which drift in 
the direction of the wave with the speed of sound 
and so remain in a constant phase of the wave. 
For a spherical energy surface these electrons 
lie on a ring perpendicular to the direction of the 
sound propagation, and since v, >v., this ring 
of electrons is nearly an equatorial one.* 

As has been shown previously,® the tempera- 
ture dependence of the longitudinal wave attenua- 
tion in superconductors is determined in the 
Bardeen-Cooper -Schrieffer theory by the energy- 
gap temperature variation. By this theory the 
ratio of the superconducting attenuation coefficient 


(a,) to that in the normal state (a,, @ constant) 


is given by 


€/kT +1), 


a. /a, =2(e (1) 
where 2e€ is the temperature-dependent energy 
gap. Thus measurements of a,/a, as a function 
of T for waves in various directions in a single 
crystal, by selecting different groups of electrons, 
should show up any significant anisotropy of €. 

We have made measurements of longitudinal 
wave attenuation at frequencies up to 80 Mc/sec 
in oriented single crystals of tin at temperatures 
down to 1.00°K. From these observations it 
appears that 2¢,, the energy gap as T-0, varies 
at least between 3.2k7,. and 4.3kT,.. The sam- 
ples were oriented such that propagation was 
along [001], [100], and [110]. Figure 1 shows 
the measurements along two of these directions 
and displays the extreme behaviors. Observations 
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FIG. 1. Normalized longitudinal wave attenuation in 
tin. Propagation along [001] is at 72 Mc/sec; that along 
[100] is at 42 Mc/sec. 
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reported earlier by Bohm and one of us,° which 
were made on a different [001] crystal with dif- 
ferent equipment, agree closely with the [001] 
measurements shown here. It should be noted 
that in all the measurements made, the condition 
>» 1 was amply satisfied (q/ is estimated to be 
20 or more). 

The energy gap anisotropy is most effectively 
demonstrated if log(a s/%,,) is plotted versus 
T/T since this curve should approach a straight 
line for large T,/T, the slope being determined 
by €,. In order to do this, however, an extra- 
polation of a , toT =0 must be made in order to 
find a,, and to eliminate attenuations of nonelec- 
tronic origin. This was done by choosing that ex- 
trapolation which gave the best straight-line 
limit when log(a ,) is plotted vs T/T. The meas- 
urements for the three orientations when plotted 
in this way are shown in Fig. 2. The estimates 
of 2€,, determined by the straight-line portions 
of the curves, are as follows: [001], (3.2 +0.1)kT,; 
[100], (4.3+0.2)k% ,; and [110], (3.8 +0.1)kT... 
Naturally, the attenuation in any one direction 
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FIG. 2. Plots of log(as/am) vs T/T for longitudinal 
waves along different crystalline axes in tin. Propaga- 
tion along [110] is at 76 Mc/sec. Measurements for 
the other directions are those shown in Fig. 1. 






would not necessarily reflect a single value for 
2¢€, since the scattering involves an average 
around a ring of the Fermi surface (or perhaps 
more than one in tin). Thus the values for 2e, 
given above probably represent only approximate 
values for dominant parts of the surface. How- , 
ever, the large qualitative differences observed 
seem most reasonably explained by energy gap 
variations of at least the magnitude indicated 
above. Such differences cannot be explained by 
variations among the samples since the same 
sample was used for measurements along [001] _ * 
and [110]. Furthermore, measurements in one 
crystalline direction at different frequencies and 

on different samples gave very similar results. 

It should be noted that Eq. (1) does not give a 
very good fit to the experimental points shown 
in Fig. 1 when the single energy gaps quoted 
above are used. In particular, the experimental 
decrease just below T e is steeper than indicated 
by these energy gaps. This suggests that in each 
direction of propagation there are some electrons 
involved which have an energy gap higher than 
the value estimated by the slope of Fig. 2. Thus 
one obtains a good fit to the experimental curve 
for [001], for example, by assuming that 80 % of 
the electrons have 3.2kT , and 20 % have 4.3kT,, 
for 2e,. 

Similar measurements are being made with 
shear waves. Because of the introduction of a 
plane of polarization, these should be more selec- 
tive in choosing electrons on the Fermi surface 
than are longitudinal waves. 








‘Work supported by the U.S. Air Force Office of 
Scientific Research. 

tOn leave from the University of Oslo, Oslo, Norway. 
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‘This possiblity, which has been suggested by a 
number of persons, was discussed at the International 
Conference on the Electronic Properties of Metals at 
Low Temperatures, Geneva, New York, August 25-29, 
1958 (unpublished). 

‘It is probable that the group velocity of a given group 
of electrons becomes modified in the superconducting f 
state [see Bardeen, Rickayzen, and Tewordt, Phys. 
Rev. 113, 982 (1959)]. This effect, however, does 
not alter the argument that propagation in different 
directions involves different electrons. 

*R, W. Morse and H. V. Bohm, Phys. Rev. 108, 

1094 (1957). 
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SPECIFIC HEAT MEASUREMENTS AND THE ENERGY GAP IN SUPERCONDUCTORS* 


Leon N. Cooper 
Brown University, Providence, Rhode Island 
(Received May 13, 1959) 


It has recently been pointed out by Boorse’ that 
the low-temperature (7,,/T > 4) specific heats of 
various superconductors seem to deviate from 
the theoretical curve calculated by Bardeen, 
Schrieffer, and the author.? We should like to 
suggest that the theory of superconductivity pro- 
posed in reference 2 contains a possible explana- 
tion of these deviations. 

If instead of the constant electron-electron in- 
teraction used in reference 2 for the detailed 
calculation of electrodynamic and thermal proper - 
ties, one uses the electron-electron interaction 
Vik, k’) which is nonconstant and anisotropic as 
might be expected in actual metals, one finds the 
results of reference 2 altered only in that the 
square of the energy gap ¢€,7(T) is everywhere 
replaced by 55*(k, T), while the criterion for the 
existence of the superconducting state is that the 
integral equation 


ae =,,, _ tanh[(8/2)E(k;7)] 
§(k) = 27 5 Vie — eT) (1) 





nave a nonzero solution for ¢. Here 8 =1/kT, 
E(k‘T) =(e7 + g5*)¥2, and Vek = Veg due to argu- 
ments of time-reversal invariance. 

It is seen that s5*, which now plays the role of 
the square of the energy gap, is a function of k 
as well as of 7; thus the energy of a single- 
particle excitation depends upon the direction as 
well as the magnitude of its momentum vector.* 
This might result in the appearance of single- 
particle excitations which seem [e.g., in absorp- 
tion experiments] to lie below the “energy gap.” 
If, for example, in a particular material there 
was a small region of the Fermi surface over 
which the energy gap was smaller than that over 
the rest of the surface, the single-particle excita- 
tions in this region might appear to lie below the 
“energy gap” determined by the gross behavior 
of the entire surface. 

Further one would expect anisotropies in meas- 
urements which pick out a direction in a crystal 
and which are sensitive to excitations from selected 


portions of the Fermi surface. One such would be 
the absorption of ultrasound in a superconductor. 
Recently Morse, Olsen, and Gavenda‘ have in 
fact observed a variation of the rate of ultrasonic 





attenuation with direction in a single crystal of 
tin. 

With regard to the low-temperature specific 
heat measurements, the “upward curvature” of 
the specific heat of aluminum pointed out by 
Boorse in the region T ft 24 would be another 
consequence of a variation in the energy gap over 
the Fermi surface. To provide a very simple 
illustration: if the energy gap had a value €,, 
over one region of the Fermi surface and a value 
of €,, over the rest of the surface, the logarithm 
of the electronic specific heat of the supercon- 
ductor would be expected to assume roughly the 
following form at low temperatures: 


C,s (3) 
aden Tag 3)ni —< 
in(-2 constant + 5 1n r 


c 


€,, T €o2 T 
c c 
vin, exe(— 7 Z)+a08 “or T/t (2) 
c c 


If a, >a, and €,,/kT > €,/kT,, the first term in 
the square brackets would dominate the second 
for T-/T small, while the second would dominate 
for T,/T large. The logarithm of the specific 
heat would then change its slope in going from 
the region of dominance of the first term to that 
of the second, producing an upward curvature. 

One test of these ideas could be made by a 
direct comparison of the variation of the energy 
gap in a given material as determined by the 
ultrasonic attenuation and by the low-tempera- 
ture electronic specific heat. 








This work was supported in part by the U.S. Atomic 
Energy Commission. 

'H, A, Boorse, Phys. Rev. Letters 2, 391 (1959). 
See also for references to the measurements. 

*Bardeen, Cooper, and Schrieffer, Phys. Rev. 108, 
1175 (1957). 

‘This fact, which has been recognized by several 
other authors, among them P, W, Anderson, A, B, 
Pippard, and V. Heine, was discussed at the Inter- 
national Conference on the Electronic Properties of 
Metals at Low Temperatures, Geneva, New York, 
August 25-29, 1958 (unpublished). Anderson has also 
pointed out the relevance of a possible anisotropy of 
the energy gap to specific heat measurements. 

‘Morse, Olsen, and Gavenda, preceding Letter [Phys. 
Rev. Letters 3, 15'(1959)]. 
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NUCLEAR QUADRUPOLE RESONANCE IN METALLIC INDIUM* 


R. R. Hewitt? and W. D. Knight 
Department of Physics, University of California, Berkeley, California 
(Received May 28, 1959) 


Subsequent to the observation of the nuclear 
quadrupole resonance in gallium,’ a resonance 
was also found in indium. Since the nuclear spin 


of In‘** is 9/2, one expects to find four resonances 


with frequencies in the ratio of 1:2:3:4 in a tetra- 
gonal crystal. 

Following a suggestion of R. H. Hammond that 
the lattice of the indium sample was probably 
strained badly, another sample was prepared in 


which dry, annealed indium spheres (2-4 microns 


diameter) were mixed loosely in quartz powder. 
The other lines for In'"*, and also the weaker 
ones for the less abundant isotope, In'“*, were 
then successfully detected. The quadrupole 
coupling constants for In**® and In’"* are 45.24 
+0.02 Mc/sec and 44.60+0.1 Mc/sec, respect- 
ively, at 4.2°K. The ratio of these frequencies 
agrees with the corresponding ratio measured by 
a different method.” 

The average temperature coefficient of the 
quadrupole coupling constant between 4.2°K and 
77°K is 5x10-*/°K. This is to be compared with 
the value 9x10~°/°K for gallium, and 6x10-°/°K 
for niobium in ferroelectric KNbO, in the same 


temperature range. 

The result is of interest for several reasons. 
First, the nuclear quadrupole resonance may be 
studied in the superconductor. Second, the con- 
tribution of the quadrupole energy to the low- 
temperature specific heat* may be calculated. 
Finally, the temperature dependence of the quad- 
rupole coupling may be compared with that ob- 
tained‘ for In’™ in the experiment on y-ray asym- 
metry. 

We wish to thank M. Pomerantz and R. Ham- 
mond for helpful discussions. 





*This research was partly assisted by the Alfred P. 
Sloan Foundation and the Office of Naval Research. 

T Alfred P. Sloan Postdoctoral Fellow on leave from 
the University of California at Riverside, California. 

1Knight, Hewitt, and Pomerantz, Phys. Rev. 104, 
271 (1956). 

?P. Kusch and T. G. Eck, Phys. Rev. 94, 1799 
(1954). 

3G. Seidel and P. H. Keesom, Phys. Rev. Letters 
2, 261 (1959). 

‘R. M. Steffen, Bull. Am. Phys. Soc. 2, 231 (1957). 





THEORY OF NEGATIVE-MASS CYCLOTRON RESONANCE 


D. C. Mattis and M. J. Stevenson 
Research Laboratory, International Business Machines Corporation, Poughkeepsie, New York 
(Received May 19, 1959) 


In recent communications’»* Dousmanis et al. 
suggest that the re-entrant energy contours of 
the heavy-hole bands in Ge and Si® could contri- 
bute a negative resistivity component to the over- 
all resistivity of these materials, possibly re- 
lated to a nonequilibrium density of carriers in 
these re-entrant states. Such a phenomenon 
could be observed under cyclotron resonance 
conditions (#4, along [100] axis, electric field 
transverse to Hac) as an emission spectrum cen- 
tered about the cyclotron frequency of these 
states, or as a reduction in the strength of the 
background absorption at that frequency. 

We have examined this idea in some detail, and 
find that except under very restrictive conditions, 
these states will in general add a normal, posi- 
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tive resistivity contribution to the material in- 
dependently of their population. Therefore, in- 
disputable experimental evidence of any emissive 
properties at the cyclotron frequency of the car- 
riers in these “dimple” states will be shown to 
reveal more about the nature of conduction in the 
valence band of Ge and Si than merely the exist- 
ence of these re-entrant energy contours, or de- 
viations from thermal equilibrium, or a combi- 
nation of these. 

As an introduction to the analysis consider the 
expression for the integrated absorption coeffi - 
cient A, calculated for a hypothetical system of 
carriers of constant, isotropic effective mass 
m . Under crossed dc magnetic and rf electric 
fields, the carriers will be stimulated by the 
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latter to make transitions between Landau levels 
i (denoted by quantum number z). We obtain 


’ A=(ne?/2m*)> (n+ 1)[f(n, E) - f(n+1,8)] 
= (ne*/2m*)>> f(n, &) = 1e7N/2m*. (1) 


' The sum in Eq. (1) is over all quantum numbers, 
only ” changing in a transition. The result we 
obtain is a sum rule independent of the magnetic 
field, namely, that the integrated absorption de- 
pends only upon the total density of carriers, N. 

Were it possible to isolate carriers in a region 
of opposite band curvature where the energy is 
given by E=E, -p*/2m*, such that over the whole 
range where carriers exist the effective mass 
-m is constant, isotropic, and negative, the in- 
tegrated absorption coefficient would be 


A=(ne?/2\m*|)> (n+1)[f(+1, © -f(n, BE] 
= - 7e?N/2|m* |. (2) 


This system of negative masses is emissive, 
just as the system of positive masses was ab- 

sorptive. Now, it might be argued that the sign 
of the mass determines whether a system of par- 
ticles is emissive or absorptive. We believe this 
point of view to be unsound. Our arguments, to 
be developed below, will indicate that in the case 
of distributions which can be characterized by a 

function of energy alone, emission will occur 
only if there is a maximum energy cutoff on the 
states of interest (such as E, in the case of nega- 
tive masses discussed above) combined with an 

‘ unbalanced “inverted” distribution. These argu- 
ments do not hold when the distribution is aniso- 





ow) = 


tropic, as is the case for the NEMAG (negative- 
effective mass amplifiers and generators) sug- 
gested by Kromer.* However, in experiments 
where holes are optically excited, thermal scat- 
tering of the holes along a surface of constant 
energy is presumably rapid enough that the po- 
pulation of states along such a contour may be 
taken to be constant. We idealize the experiments 
of Dousmanis et al. by means of the following 
assumptions: 

1. Holes are excited up to the optical phonon 
energy Eopt- (Those created at higher energies 
will rapidly emit optical phonons and fall to this 
energy or below.) 

2. A negative transverse mass -m and posi- 
tive longitudinal mass m, will be assumed con- 
stant for carriers within narrow cones of half- 
angle a whose axes lie along the [100] axes of 
Ge and Si. Beyond this cone, the heavy holes 
under consideration have a positive effective 
mass mp, whose magnitude is sufficiently differ- 
ent from mr that the conductivity at frequencies 
near w7=eH/m-rc will be determined largely by 
the properties of the carriers within the cones. 
We shall consider the possibility that m7, m7, 
or a vary with energy. 

3. Steady state is described by a distribution 
function f(E) dependent only on energy. The 
microwave fields are a small perturbation, whose 
effects we shall describe by the Boltzmann equa- 
tion. Thus, under cyclotron resonance conditions, 
with the magnetic field along the axis of one of the 
cones, we find a frequency-dependent conductivity 
given by* 


-Ne* | “do sin?6 {dp p*v*[df(E) /dE \(r + iwt”)/[wt? + (1+ iw] 





2{"d0 sind {dp p*f(E) 


-Ne*a* [dE E™*(df/dE)(1 + iwt*)/[(w pt)? + (1 +iw)?] 


(3) 





The second expression follows if a*<<1. The 
total density of carriers within the cone is taken 
to be N, and the momentum relaxation time 7 

i obeys wT >> 1 for a well-defined resonance. A 
dispersion relation valid near the axis of the 
cone is 


E= 2*(cos*8/m | - sin’6/m_,) £ 3p*/m >, 


=20*(m , cos*é -m,,sin*é)=30*m_ , (4) 






(a* << 1) 








2m , [dE E™* f(E) 





and was used to express the variables of Eq. (3) 
in terms of the energy. 

If w7 and 7 vary slowly with energy compared 
with f(E), the numerator of Eq. (3) can be inte- 
grated by parts, and we obtain 
3Ne*a* T+iwt* 

—_—— Q (5) 
4m w 7 + +iwt)" 


L 


It is easily verified that Re[o(w)]>0, so that we 
have a direct check that under reasonable assump- 
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tions the conductivity, and hence the absorption 
for the system under consideration is positive. 

If all the hypotheses we stated are not valid, 
there will exist possibilities for these negative 
masses to amplify resonant radiation. We en- 
visage two distinct situations: the first involves 
a cutoff energy E-<Epopt, such that for E>E, 
the shape of the energy contours has changed 
sufficiently that no particles® above this cutoff 
energy can contribute to the resonance at wy. If 
such is true for the energy band structure of the 
heavy holes in Ge or Si, emission would be ob- 
served at wy provided 


E, 
, GEE df/dE >0. 


This last condition is tantamount to an inversion 
of the population of these low-lying energy levels, 
which could be achieved by a mechanism such as 
selective recombination of slow carriers. 

If the band structure does not permit such a 
drastic change in m7, a second situation is con- 
ceivable whereby the distribution function while 
still of cubic symmetry shows some anisotropy. 


If, for example, 


f=f,(E)+f,(151), (6) 


those particles distributed according to /, will 
absorb radiation, whereas the contribution of the 
carriers with negative transverse mass -m 7 
distributed according to f, can be shown to be 
emissive, and peaked about the frequency w7. 
Whether emission or absorption predominates at 












that frequency will depend on the relative mag- 
nitudes of f, and /,. 

In this connection we propose that light may not } 
provide the best mechanism for obtaining a se- 
lective excitation of such carriers. A possible 
alternative excitation technique could utilize 
microwave impact ionization of impurity centers 
in p-type Ge by resonant absorption at the cyclo- = 
tron frequency wp of the positive-mass heavy 
holes. With a double resonant cavity system, the 
changes in the rf conductivity may then be studied 
with a low-power probing signal at the resonance 
frequency of the negative-mass holes, wy. Such 
a technique would result in the presence of posi- * 
tive-charge carriers only, and with proper cir- 
cular polarization the effects of the negative- 
mass holes can be separated from those of the 
positive-mass holes. 

In concluding, we express our thanks to col- 
leagues of this laboratory—to Mr. W. Kopka for 
computational work on the IBM 704 computer, 
and to Dr. S. Koenig, Dr. R. Landauer, and Dr. 

E. Adams for their interest and helpful sugges- 
tions. 








1G. C. Dousmanis, Phys. Rev. Letters 1, 55 (1958). 

*G. C. Dousmanis et al., Phys. Rev. Letters 1, 404 
(1958). 

‘Dresselhaus, Kip, and Kittel, Phys. Rev. 98, 368 
(1955). 

‘H, Kromer, Proc. Inst. Radio Engrs. 47, 397 (1959). 

5For method of derivation, see, for example, D. 
Mattis and G. Dresselhaus, Phys. Rev. 111, 403 (1958). 

®This is an idealization. In reality we would not nec- 
essarily expect a sharp cutoff and there would be over- | 
lap in the spectra of states above and below E c 





ROLE OF NEGATIVE EFFECTIVE MASS IN NEGATIVE RESISTANCE 


Peter Kaus 
University of Southern California, Los Angeles, California 
(Received May 21, 1959) 


Experiments on cyclotron resonance of nega- 
tive-mass carriers in germanium have been re- 
ported.’ The distribution of carriers (heavy 
holes) is not an equilibrium distribution, since 
the carriers were produced by optical excitation 
partly from acceptors and partly across the gap. 
The authors find a new resonance for the oppo- 
site sense of circular polarization to the positive- 
mass resonance line and the resonance corre- 
sponds to a decreased absorption. It has also 
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been proposed to construct amplifiers using ne- 
gative-effective-mass carriers in a nonequilib- 
rium distribution maintained by auxiliary power 
sources.” 

The decreased cyclotron resonance absorption 
as well as the proposed amplifiers depend on the 
negative resistance aspect of negative-effective- 
mass carriers. It is therefore of interest to in- 
vestigate the relationship between negative effec- ‘ 
tive mass, carrier distribution, and negative 
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resistance. 
The stationary states in a crygtal are described 
by Bloch functions, ¥,,(K, F)=e**-Tu, (K,), with 
energy eigenvalues E,,(k). Dropping the band 
subscript m, the corresponding group velocity is 


v(k) = (1/n)V,E(). (1) 


Within well-known limits of validity the motion 
of an electron under the influence of a force F is 
given by 


dk/dt =(1/n)F. (2) 


Differentiating (1) with respect to the time and 
combining with (2), the equation of motion be- 
comes 


*x 
dv /at=2,(1 /m dF > (3) 


where a component of the reciprocal effective- 
mass tensor 


(1 /m"),, = (1 /WE)@PE / 8k 8h (4) 


may be positive or negative. 

A one-dimensional example is given in Fig. 1. 
Since every band has a lower and upper bound in 
energy, there is a positive and negative effective- 
mass region in every band. Consider a carrier 
in the negative effective-mass region [Fig. J (a)]. 
If a force in the positive x direction is applied, 
the carrier moves in the positive k direction re- 
gardless of the energy band picture. It therefore 
moves to a region of smaller velocity or it accel- 
erates in the negative x direction. If we think of 
a current associated with this motion, then the 
change in current is in the negative x direction 
or opposite to the applied force. A single car- 
rier with negative effective mass, therefore, 
contributes to negative resistance. 

For the same band consider now an equilibrium 
distribution at low temperature. Assume that the 
Fermi level is high up in the band [Fig. 1(b)]. 

The inflection point has been drawn at a low k 
value, so that many more negative-mass states 
are occupied than positive-mass states. At first 
glance it may appear that consequently this dis- 
tribution would exhibit negative resistance to all 
applied forces at thermal equilibrium, or that it 
would emit rather than absorb. 

This conclusion would be in obvious contradic - 
tion to the second law of thermodynamics. The 
sides of the crystal would charge up, producing 
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FIG. 1. Schematic representation of energy and 
velocity as a function of k for a simple one-dimensiona 
band. The black dots represent occupied states. The 
black dots with a vertical line are states occupied at 
t=0, but unoccupied at t= At; the black dots with a 
horizontal line are states unoccupied at t=0, but occu- 
pied at ¢=Af, 


a field which would enhance the current. The 
original field is then seen to be unnecessary to 
produce a current which can be made to do work 

In the present example it is easy to see that th: 
resistance of the distribution in Fig. 1(b) is in 
fact positive. Every carrier moves in the posi- 
tive k direction at the same rate, according to 
(2). The total effect, therefore, is to add posi- 
tive velocity to the total distribution of velocity. 
This conclusion is evidently seen to be independ- 
ent of the shape of the energy band or the positio: 
of the Fermi level. However, this simple pictur: 
holds only at very low temperatures. 

In the next section some necessary conditions 
for negative resistance are developed for a 
general distribution at a finite temperature. 

The contribution of a quantum state to the cur- 
rent density is given by 


TR) = (1/L*)(e/n)¥,E®. (5) 


The time derivative of the current contribution 
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(5) is obtained using (1) and (2): 


4,® ~ ii)ari*; ” 


Let g(k) be a general steady nonequilibrium 
distribution function. The function is maintained 
by the application of some auxiliary forces not 
contained in F;. We now want to test the response 
of the distribution to a small incremental force 
F;. 

Jone change in the total current density is ob- 
tained by integrating (6) over the distribution 


at ~ \Z5 jJ ee) 


If in particular we chose an electric field & 
for the incremental field, the derivative of the 
current density becomes 


a 


In the following considerations nothing new 
emerges if the relaxation time 7 is introduced. 
The derivative dJ;/dt is then to be replaced for 
small manemiitel fields é by the expression 
4J;/T, where AJ; is the increment of current 
density over the steady-state current density. 

The sign of the diagonal conductivity term, 

0;;, is therefore given by the sign of the integral 


8v {© 
S..= “a HOB (9) 


ii 
We can write the integral (9) in the form 


54° /f 5.0 k .)dk dk 


where 


i i’ ,»k 
Su Sj%*) = Paget z a aN yh ak. 


The corresponding resistance is certainly posi- 
tive if S; ii(2j,}) is positive for all k; and kj. 

nausea (10) by parts and taking advantage 
of the equality of v; at opposite edges of the zone, 


(10) 
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we obtain 
og 
= 11 
Si ke, k)=- v; me (11) 
or 
= og 
jh = - Be Bes (12) 


From (12) it is clear that if the distribution 
function always decreases with increasing energy 
along a line (kj, 7) constant, then the correspond- 
ing contribution 5;; (k;,k) to the conductivity is 
positive. If this is the case for all lines (k;,k ? 
constant, then S;; is positive and the corréspond- 
ing resistance is positive. 

A necessary condition for negative resistance 
is therefore the following: “A distribution g(k) 
cannot correspond to negative resistance in the 
i direction unless 8g/8E is positive along some 
line (kj »k)) constant.” 

The Teversal of population demanded by this 
necessary condition depends in no way on the 
band structure. An obvious conclusion is that 
the equilibrium distribution [Fig. 1(b)] corre- 
sponds to positive resistance regardless of the 
number of negative-mass carriers or the posi- 
tion of the Fermi surface. 

It would be misleading to conclude that negative 
mass has therefore no bearing on negative re- 
sistance. We can put Eq. (9) into the form 


Si,=h i [e(k)/m *,™ 1 (13) 


Singe g(k) is positive definite, another neces- 
sary condition for negative resistance is seen to 
be: “A distribution g(k) cannot correspond to 
negative resistance in the i direction unless some 
states corresponding to negative effective mass ' 
m*, are occupied.” 

A distribution consisting entirely of negative- 
mass carriers [Fig. 1(a)] always corresponds to 
negative resistance. This is of course consistent 
with the first condition, since the lowest energy 
states along a line in k space always correspond 
to positive effective mass in the direction of the 
line. If only negative-mass states are occupied, 
the necessary reversal of population must have 
taken place. 

It must be pointed out here that the reversal of 
population demanded in the first necessary con- 
dition is not a reversal of the total population at 
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one energy. It is not necessary that 


} dg ag *, 
aE-™ | me te 
$s Jj 


with the integral performed over a surface of 
constant energy, be positive for some such sur- 
face. The necessary condition for negative re- 
sistance in the i direction only demands that 
g/8E along some line in the i direction be posi- 
tive. In Ge and Si with warped heavy-hole bands, 
the necessary reversal may only correspond to 


a redistribution from the thermal equilibrium 
distribution along constant-energy surfaces. 

Most of these conclusions were simultaneously 
and independently reached by Dr. C. Kittel. Iam 
very grateful to him for making his results avail- 
able. I would also like to thank Dr. F. Herman | 
for several discussions. 





‘Dousmanis, Duncan, Thomas, and Williams, Phys. 
Rev. Letters 1, 404 (1958). 

7H. Krémer, Phys. Rev. 109, 1856 (1958); Proc. 
Inst. Radio Engrs. 47, 397 (1959). 





BIMOLECULAR ELECTROLUMINESCENT TRANSITIONS IN GaP 


E. E. Loebner and E. W. Poor, Jr. 
RCA Laboratories, Princeton, New Jersey 
(Received March 19, 1959; revised manuscript received June 1, 1959) 


In this Letter we describe direct observations 
of bimolecular band-to-band transitions at recti- 
fying junctions in single crystals of gallium phos- 
phide and indium phosphide. 

The presence of bimolecular recombination 
kinetics is deduced from the nearly square de- 
pendence of the emitted radiation on the injec- 
tion current. Additional evidence for band-to- 
band electroluminescence is supplied by an anal- 
ysis of the spectral distributions of the absorption 
coefficient, the electroluminous emittance, and 
the photovoltaic response of gallium phosphide 
crystal rectifiers. These band-to-band radiative 
transitions differ in a number of ways from those 
occurring via localized levels. 

The above-mentioned square dependence has 
been discussed by van Roosbroeck and Shockley” 
and looked for without success in germanium 

| by Newman® who found a linear dependence. So 

far only Aigrain and Benoit* observed a square 

dependence in germanium, when a very sensitive 
detection technique was used. Braunstein,® who 
investigated band-to-band radiative transitions 

of injected carriers in Ge-Si alloys, GaAs, GaSb, 

and InP, found a linear dependence in all cases, 

| in agreement with Newman’s observations.* 

Our GaP samples were grown by I. Hegyi of 
these Laboratories from gallium solution by the 
method of Wolff, Keck, and Broder® using a 5- 
day programmed heat-cool cycle with a peak 

temperature of about 1200°C. The crystals were 

flat hexagonal platelets with major faces parallel 






to (111) planes and an area of a few square milli- 
meters. They were a few tenths of a mm thick. 
The crystal structure, determined by J. White, 
was that of sphalerite with a lattice constant of 
5.449 A. The impurity with the highest concen- 
tration was copper with approximately 5-50 parts 
per million. Care was taken to avoid using crys- 
tals that had inclusions of free gallium or free 
phosphorus. 

Undoped GaP crystals appear to contain grown 
junctions beneath the surface, which resemble 
somewhat those found in single crystals of silicon 
carbide grown by commercial processes.” A 
highly n-type conducting layer is separated from 
the bulk by a rectifying junction and most likely 
also by a nearly intrinsic layer. Ohmic contacts 
to the crystal bulk can be made by alloying 
indium through the junction or by electrical burn- 
out of the junction layer. 

The opto-electronic characteristics of the rec- 
tifiers were displayed directly on an oscilloscope. 
The radiation was detected by RCA 1P21 and 
7102 multiplier phototubes after having been spec- 
trally analyzed by a Perkin Elmer monochroma- 
tor. Transmission data were obtained using a 
model 14M Cary recording spectrophotometer. 

Figure 1 is a photograph of five oscilloscope 
traces displaying the mutual and time dependence 
of the GaP rectifier parameters light (L), cur- 
rent (J) and voltage (V) for a 60-cps sinusoidal 
8-volt rms excitation. Trace A is a plot of L 
vs I, which shows essentially a square depen- 
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FIG. 1. 


Photographs of oscilloscope signals obtains: 
from an electroluminescing GaP crystal rectifier. 
A-electroluminous emittance vs current; B-current 
vs voltage; C—electroluminous emittance vs voltage; 
D-current ys time; E —electroluminous emittance vs 
time. 


dence of light on current. Traces B and C show 
the J/-V and L-V characteristics, respectively. 
It is seen that current and light appear only on 
the forward bias portion of the cycle. Traces 
D and E show the time dependence of / and L, 
respectively. The light at all times is perfectly 
in phase with both the current and the voltage. 
Microscopic observations show the light to be 
emitted uniformly over an area of about 10 cm 
from a depth of several microns below the sur- 
face. In the cases on hand, only the light clearing 
the electrode whisker participates in the data 
contribution. 
Various results pertaining to the spectral char- 
acteristics of our samples are seen in Fig. 2. 
Our optical absorption measurements lead us to 
a room temperature estimate of an energy gap 
value for GaP of 2.20 to 2.21 ev. This estimate 
results from a fit of the Macfarlane-Roberts ex- 


pression® to an a vs h plot for values of the absorp- 


tion coefficient, a, between 10 and 150 cm™ at 
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FIG. 2. Spectral opto-electronic responses of GaP 
and InP. Radiance indicated by - Short-circuit 
photovoltaic current indicated by ----. Letters in 
parenthesis indicate the following workers and refer- 
ences: O-—Oswald (reference 13); N—Newman (refer- 
ence 13); FO-—Folberth and Oswald (reference 9); LP- 
Loebner and Poor (this Letter). 





several temperatures. Assuming the presence 
of a single phonon assisting in the optical transi- 
tion, its energy is estimated to be 0.026 ev, 
which gives an equivalent temperature of about 
300°K. Excess absorption at lower energies and 
higher temperatures is probably due to optical 
transitions at dislocations or some other crystal 
imperfections. As seen from Fig. 2, the energy 
gap value obtained by us is in good experimental 
agreement with that quoted in the literature.® 
The dashed curve on Fig. 2 represents the spec- 
tral response of the short-circuit photovoltaic 
current in the same junction for which the electro- 
luminous spectrum has also been obtained. The 
photoelectric current flows in the direction op- 
posite to that producing electroluminescence, 
while the photovoltage has the same polarity as 
that producing electroluminescence, in agreement 
with the cases of SiC and CdTe.” The analysis of 
the photocurrent spectral response is identical to 
that of Wysocki, Loferski, and Rappaport." For 
hv just greater than E_, In(Q, max - Q)) is plotted 
against a. This gives a straight line, whose slope, 
10~* cm, is a measure of the distance that the 
photons pass through GaP. The relatively rapid 
falloff of Q, at high energies indicates a carrier 
diffusion length very much shorter than 10~? cm." 
The spectral radiance of the very same GaP 
rectifier during electroluminescence is also seen 
on Fig. 2. The radiation peaks at photon energies 
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of 2.20+0.02 ev with a half-width of a few kT. 
This value is a few hundredths of an electron volt 
lower than the peak obtained from the van Roos- 
proeck-Shockley theory’ using the calculated a’s 
(derived from the Macfarlane-Roberts analysis), 
but agrees well with the peak obtained using the 
above-discussed additional low-energy absorp- 
tion. It is interesting to note that a comparison 

of the high-energy side of the spectral curve of 
electroluminous emittance to the van Roosbroeck- 
Shockley theory indicates a photon travel through 
GaP of about 10~* cm. This is in rather good 
agreement with the value obtained from the photo- 
current vs hv analysis. Some of the additional 
low-hv emission is most likely due to defect sites 
mentioned above. These and self-absorption cause 
the peak to shift somewhat to lower energies. 

We conclude that most of the electroluminescent 
phenomena on hand are quite distinct from those 
originating from recombinations via defect states 
as reported previously.’»? In all the latter cases 
the photon distribution peaks at energies lower 
than the energy gap value obtained from free 
carrier recombination electroluminescence re- 
ported in this Letter, or from the analysis of 
other optical data. Also in all defect center re- 
combination cases a linear dependence of the 
light on current has been found. The agreement 
between the above-discussed characteristics of 
electroluminescence, optical absorption, and 
photoelectric response leads us to identify the 
electroluminescent phenomenon as one primarily 
due to band-to-band bimolecular transitions. 

We have also seen bimolecular recombination 
kinetics for radiative transitions of carriers in- 
jected across a rectifying junction in indium 
phosphide. Braunstein® reported the presence of 


recombination radiation in InP at an energy close 
to the band gap. We have obtained the spectral 
distribution of this radiation. It is shown in 
Fig. 2. The position of the peak at 1.25 ev is in 
good agreement with the band gap values reported 
in the literature.** 

The authors wish to express thanks to R. Braun- 
stein, A. R. Moore, and J. I. Pankove for helpful 
discussions. 
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ELECTRON SPIN RESONANCE EXPERIMENTS ON SHALLOW DONORS IN GERMANIUM 


G. Feher, D. K. Wilson,* and E. A. Gere 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 4, 1959) 


Electron spin resonances on shallow donors in 
silicon have been studied by several authors.’ 
In germanium there are various factors which 
make the observation of electron spin resonance 
signals more difficult and, so far, no resonances 
have been reported on shallow impurity states. 
The recent conflicting experimental results per - 
taining to the donor ground state in germanium*»® 


led us to reinvestigate the electron spin resonance 
behavior in this material. In the present Letter 
we wish to report on the observation of electron 
spin resonances from both bound and nonlocalized 
electrons in germanium. The results indicate 
that the ground state of the arsenic and phosphorus 
donors in germanium is a singlet. 

The experiments were performed at temper- 
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FIG. 1. Electron spin resonance signal from arsenic 


and phosphorus donors in germanium with the magnetic 
field H parallel to the [100] direction. T= 1.3°K; 
ve = 9000 Mc/sec. 


atures between 1.3°K and 4°K and microwave fre- 
quencies of ~9000 Mc/sec and ~24000 Mc/sec. 
The details of the spectrometer have been de- 
scribed elsewhere.® 

At low donor concentrations the electrons are 
localized on the impurity and a resolved hyper - 


fine interaction with the donor nucleus is observed. 


This is illustrated for arsenic (/ = $) and phospho- 
rus (J= 4) doped germanium in Fig. 1. The exper- 
imental results on the line width AH (full width 
between inflection points of dx’’/dH), hyperfine 
splitting, and g values are summarized in Table I. 
The experiments on the arsenic-doped sample 
were repeated at 24000 Mc/sec. The g value, 

AH (with H parallel to [100]), and hf splitting were 
found to be the same as at the lower frequency. 
For comparison purposes the data on silicon’,® 
are also listed in Table I. 


———__ 


The g shift in germanium is about two orders 
of magnitude larger than in silicon. It is pro- 
portional to A/AE, where A is the spin-orbit 
coupling which in the valence band of germanium 
is an order of magnitude larger than in silicon 
and AE is the energy separation of the nearest 
bands that are admixed. The large observed 
g shift may have contributed to the fact that 
these lines have previously been missed. 

The line widths quoted in Table I were obtained 
with H parallel to the [100] direction. When H 
deviates from this direction an appreciable line 
broadening occurs. This broadening together 
with the necessary small donor concentration 
in germanium makes the observation of the lines 
more difficult than in silicon. The broadening 
may arise, either from the anisotropic part of 
the hf interaction with the germanium nuclei or 
from a spread in g values due to a local pertur- 
bation of the donor wave function’ (e.g., due to 
strains or overlaps). The observed Gaussian 
line shape is consistent with either picture. Ex- 
periments are under way to shed light on this 
question. 

The magnitude of the hf splitting is about a 
factor of two smaller than in silicon. By a sim- 
ple extrapolation of the theory of shallow donors 
in silicon due to Kohn and Luttinger,*»® one ex- 
pects the hf splitting in germanium to be approxi- 
mately an order of magnitude smaller than in 
silicon. This discrepancy with the experimental 
results is consistent with the large observed line 
width and indicates that the peaking of the con- 
duction band wave functions at the nuclei is about 








Table I. Summary of electron spin resonance data on phosphorus- and arsenic-doped germanium and silicon. 2 > 
Because of different passage conditions» the line width in silicon has to be multiplied by 0.85 before comparing it 


to germanium. 











Donor Line width Total hf 
concentration AH splitting | ¥(0)|? 
Donor (cm~) g value (oersteds) (oersteds) (cm~) 
Germanium 
Phosphorus 8 x 10'4 1.5631 + 0.0003 10 21 0.21 x 104 
Arsenic 5x 10% 1.5701 + 0.0003 11 107 0.87 x 10% 
Silicon 

Phosphorus 1.5 x10 1.99850 2.8 42 0.43 x 107 
Arsenic 1.810 1.99837 3.2 212 1.73 x 104 





“gee reference 1. 
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a factor of six larger in germanium than in sili- 
con. However, in deciding whether the singlet 

or the triplet is the lowest state only the pres- 
ence of an appreciable isotropic hf interaction 

is important. Kohn and Luttinger*»® have shown 
that the ground-state wave function can be made 
up of a linear combination of four wave functions 
corresponding to the four minima of the conduc- 
tion band. For three of these combinations (which 
form the degenerate triplet state) the value of 

the wave function at the origin is zero, and there- 
fore cannot give rise to the observed isotropic 
hyperfine interaction. The singlet state is formed 
from a linear combination that has a finite value 
of |~(0)|? and should give rise to a hf splitting. 
The experimental observation of a splitting indi- 


cates therefore that the ground state is the singlet. 





The absolute amplitude of the signal and its var- 
iation with temperature excludes the possibility 
that one is observing an excited state. 

At high donor concentrations” the electron is 
no longer localized on a particular impurity but 
moves throughout the crystal giving rise toa 
conduction process (“impurity band conduction”) 
even at the lowest temperatures. Under these 
conditions the hyperfine interactions are averaged 
out and a single resonance line is observed." 
The line width for such nonlocalized electrons 
in germanium is found to be anisotropic as illus- 
trated in Fig. 2 which shows the results ona 
phosphorus-doped sample (10'7 P/cm’, T =1.3°K, 
v, = 9000 Mc/sec). The width presumably arises 
from a spread in electronic g values which, in 
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ANGLE IN DEGREES BETWEEN H AND CRYSTAL AXIS [19] 


FIG. 2. Line width as a function of the orientation 
of the magnetic field with respect to the crystalline 
axes. T= 1.3°K; vp ~ 9000 Mc/sec; 10’ P/em’. In 
this impurity concentration range one observes a 
single resonance line due to nonlocalized electrons. 


contrast to silicon,’ are incompletely averaged 
out by motional narrowing.” An increase in tem- 
perature from 1.3°K to 4°K decreased the line 
width in the above phosphorus sample by ~ 10 %. 
Similarly, an increase in the concentration from 
3x10'* P/cm' to 10'7 P/cm® decreased the line 
width in the [100] direction from 6 oersteds to 

4 oersteds. These results and the observed 
Lorentzian line shape support the picture of a 
motionally narrowed line.” 

The electronic g value of the nonlocalized elec- 
trons was found to be isotropic and approximately 
equal to the g value of the bound electrons (e.g., 
for 10'’ P/cm’, g=1.5636+0.0003; and for 
10’? As/em*, g=1.5698+0.0003). This indicates 
that in this concentration of “impurity banding” 
the electron spends most of its time on the im- 
purity. 

Antimony-doped germanium was investigated 
in the concentration range between 10'* Sb/cm* 
and 3x10'* Sb/cm*. The line shape was found to 
be extremely asymmetric and the g value varied 
between ~1.6 with H parallel to the [100] direc- 
tion and 1.9 with H parallel to the [110] direction. 
This striking difference in behavior indicates 
that either the triplet-singlet levels overlap or 
that perhaps the triplet is lower in antimony. 
Strain experiments are in progress which should 
decide on this question. 

We profited greatly from discussion with many 
members of the Bell Telephone Laboratories. 

In particular we wish to acknowledge the stimu- 
lating discussions with W. Kohn, P. W. Anderson, 
and G. Weinreich. We also wish to thank G. Bem- 
ski and J. Eisinger for their help in the 24000- 
Mc/sec experiment. 
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"This would cause a mixing of the singlet and triplet 
states. However, in the [100] direction the four con- 
duction band minima are equivalent and the mixing 
would not cause a broadening. This picture is con- 
sistent with the result that the narrowest line is ob- 
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1% silicon this occurs at concentrations of the order 
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purity Bohr radii in germanium, the corresponding 
concentrations are lower by between one and two orders 
of magnitude. 

''In silicon this line was first observed by Portis, Kip, 
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ABSORPTION OF COMPRESSIONAL WAVES IN SOLIDS FROM 100 TO 1000 Mc/sec 


J. Lamb, M. Redwood, and Z. Shteinshleifer 
Electrical Engineering Department, Imperial College, London, England 
(Received May 25, 1959) 


We have determined the absorption of com- 
pressional waves in crystal quartz [X -direction], 
fused silica, germanium [100], and silicon [111]: 
the crystallographic directions specify the direc- 
tion of propagation. Measurements were made 
by the “pulse technique” using guided waves in 
cylindrical specimens.'~* The dependence of the 
amplitude absorption coefficient on frequency is 
shown in Figs. 1-4. 

In contrast to the findings of Granato and 
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FIG. 1. Log(absorption coefficient) versus log(fre- 
quency) in crystal quartz. Propagation is along the 
X-axis. o=20°C;e=-77°C. At 20°C, @ is proportional 
to f 1-82 
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Truell and their co-workers*~’ for germanium, 
we observed little difference between the absorp- 
tion values for five crystals of silicon from dif- 
ferent sources. These had dislocation densi- 
ties—revealed by etch pits—from a few hundred 
to over 10* per cm*. Two specimens were n- 
type with average resistivities of 0.1 and 14 
ohm-cm, and three were p-type with resistivi- 
ties of 3.20 and 30 ohm-cm. The axis of each 
cylinder was within 3° of the [111] direction. The 
absorption values in all five specimens agreed 
to within +5 %. 
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FIG. 2. Log(absorption coefficient) versus log(fre- 
quency) in fused silica (Thermal Syndicate, O.G. 
grade) at 18°C: a@ is proportional to f?-!2, 
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FIG. 3. Log(absorption coefficient) versus log(fre- 
quency) in germanium. Propagation is parallel to the 
[100] direction. o= 18°C; A=-158°C; the absorption is 
proportional to f!-** at 18°C. 


Granato and Truell® found a wide range of ab- 
sorption values in different specimens of ger- 
manium; their lowest values are close to, but 
slightly higher than, our values (Fig. 3). 

At 20°C the absorption coefficient, a, for sili- 
con and for germanium is approximately propor- 
tional to the square of the frequency. As a rough 
approximation, this is also true for crystal 
quartz and fused silica. 

The values for crystal quartz at 1000 Mc/sec 
give quantitative support to estimates which have 
been made using other techniques.*»* We are 
continuing this work. 

We wish to thank Mr. R. IL. Walker of the Ser- 
vices Electronics Research Laboratory and Dr. 
K. Hoselitz of the Mullard Research Laboratory 
who supplied the crystals of silicon, and Mr. S. 
E. Bradshaw of the General Electric Company’s 
Research Laboratory who provided crystals of 
germanium for this work. 





<Nw 
SS 


/ | 
? 





% (db/em ) 
Ne 





+ 
2 





0.6 fh i 


| 


0.2 fa 
































Ou! f (Mc/sec) 
100 200 400 600800 











FIG. 4. Log(absorption coefficient) versus log(fre- 
quency) in silicon. Propagation is parallel to the [111] 
direction. V=197°C;0=20°C;e=-77°C; A=-158°C. 
At 197°C, @ is proportional to f?** and at -158°C to 
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EFFECTS OF RARE EARTH IMPURITIES ON FERRIMAGNETIC RESONANCE 


IN YTTRIUM IRON GARNET 


J. F. Dillon, Jr., and J. W. Nielsen 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 11, 1959) 


The first published ferrimagnetic resonance 
work on single-crystal yttrium iron garnet’ (YIG) 
gave room-temperature line widths which were 
narrower than those observed for any other com- 
pound up to that time. However, the curve 4H(T) 
showed a maximum at 40 or 50°K, and the varia- 
tion of the field for resonance with crystal direc- 
tion,” i.e., the surface H,,,(hk/), indicated that 
the anisotropy surface at liquid helium tempera- 
tures was convoluted beyond anything reported 
previously. Since then one of the authors (JFD)* 
and others* have observed much narrower lines 
and, associated with them, AH(T7) curves® with 
much lower maxima which were, however, still 
near 50°K. These crystals did not show the wide 
excursions of Hyeg(hkl) at low temperatures. In 
view of the clear possibility that these effects 
were caused by the presence of rare earth im- 
purities, a series of YIG crystals doped with 
each of the rare earth ions were grown in order 
to identify experimentally the factors responsible 
for the line-width peak and the anomalous aniso- 
tropy. Recently Kittel® has considered theoreti- 
cally the effect of rare earth ions on ferrimag- 
netic resonance in the rare earth garnets and in 
rare-earth-doped YIG. This paper gives an ex- 
tract of the data obtained on YIG, containing 
small concentrations of rare earth ions. These 
show quite clearly that the line-width peak is at 
least predominantly a matter of rare earth im- 
purities, that it is not always at 40 or 50°K, and 
that the anomalous anisotropy is associated also 
with the same impurities. 

The crystals used were grown from flux. 
Highly polished spheres of 0.015- to 0.019-in. 
diameter made by the two-pipe technique” were 
mounted on sapphire rods with the [110] axes 
along the axis. The experiments were done near 
9.3 and 24.0 kMc/sec using TE,., mode rectan- 
gular cavities in transmission. The high-fre- 
quency cavity was situated within a double Dewar. 
Data were obtained over a wide temperature 
range by boiling away the cryogenic liquid and 
allowing the cavity to warm at a rate determined 
by the heat leak. Temperature was measured 
with a thermocouple or a carbon resistance 
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thermometer in contact with the wall of the cavity. 
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FIG. 1. Variation of ferrimagnetic resonance line 
width with temperature for several rare-earth-doped 
YIG crystals. Dopings are given in atomic percent on 
the figure. Line width is the field interval between 
points of half the maximum absorption. The data for 
YIG(0.6%Tb) were taken at 9.3 kMc/sec, the others 
at 24 kMc/sec. 


The line widths were determined from record- 
ings made of the absorption against magnetic 
field. All field measurements were made by 
nuclear magnetic resonance techniques. 

Figure 1 shows the 4H(T) curves for four dif- 
ferent samples. We shall discuss in this Letter 
only two impurity ions, terbium and europium. 
We have introduced the terbium in three different 
concentrations.* Comparing the three YIG(Tb) 
samples, we see that the line width increases 
monotonically with terbium concentration. There 
is a pronounced peak in the line width at about 
90°K for the YIG(0.04% Tb) and perhaps slightly 
lower for the YIG(0.1% Tb). The europium-doped 
sample serves as a demonstration of the very 
different effect encountered in introducing var- 
ious rare earth ions. Whereas the ground state 
of Tb*** is "F,, that of Eut++ with its less than 
half-filled shell has the singlet ground state "F,. 
Since J=0, we might expect to a first approxi- 
mation that the europium impurities will not con- 
tribute to the line width of the main iron reso- 
nance. 

Figure 2 shows 2... (hkl) in (110) for YIG(0.1% 
Tb) as well as for a nominally pure YIG made 
with a specially purified yttrium oxide. The for- 
mer shows a surface even more convoluted than 
that reported before for a sample with unspeci- 
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| FIG. 2. Field for resonance at 1.5° K, against angle 
in the (110) plane for YIG free of intentional rare earth 
doping (scale at right, f = 23192 Mc/sec) as contrasted 
with the corresponding plot for YIG doped with 0.1% 
terbium (scale at left, f=22990 Mc/sec). Line widths 
are given for several of the maxima and minima of the 

’ curve for the doped sample to show the line-width 
anisotropy. 


fied impurities.?, The curve for the YIG is 

grossly similar to that observed at room tem- 
t | perature, but there are a number extra peaks 

| several degrees wide and perhaps forty oersteds 

high superimposed on the expected simple curve. 
It appears now that a number of the rare earth 
ions will have distinctive anisotropy and H_,. eg “l) 
surfaces associated with them. A comparison of 
the two curves of Fig. 2 suggests that one of the 





impurities in our YIG is terbium present ina 
concentration of a few parts in a million. 

A. M. Clogston and L. R. Walker have joined 
one of the authors (JFD) in attacking the theoret- 
ical problem set by these plots and correspond- 
ing data on YIG doped with other rare earth ions. 

A large anisotropy in line width is observed to 
be associated with the anomalous H,,,(hk/) sur- 
face. In the case of the YIG(0.1% Tb) at 1.5°K, 
several of these values are given on Fig. 2. Thus 
within a single crystal we have available a range 
of damping parameters. 

These are being studied to gain an understand- 
ing of the physical processes involved in the 
damping associated with ferrimagnetic resonance. 

We intend to describe in detail the effects of 
small concentrations of the other rare earth ions 
in YIG on AH (hkl), H,pg(hkl), and g, allasa 
function of temperature, in subsequent papers. 

The authors acknowledge with thanks many dis- 
cussions with A. M. Clogston, E. G. Spencer, 
and R. C. LeCraw. H. E. Earl provided valuable 
technical assistance. 
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LOW-TEMPERATURE LINE-WIDTH MAXIMUM IN YTTRIUM IRON GARNET 


E. G. Spencer, R. C. LeCraw, and A. M. Clogston 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 11, 1959) 
































Measurements of the ferromagnetic resonance » 
line-width maximum in single-crystal yttrium 
iron garnet (YIG) were previously reported by 4F 
Dillon.’ His data show that the line width in- 
creases with decreasing temperature, reaches 12+ 
a maximum below liquid nitrogen temperature, 
and decreases again to approximately the room- 
temperature value at liquid helium temperatures. 
We have been investigating the nature of this 
line-width maximum at low temperatures to de- 
termine if the effect is intrinsic to ferromagnetic 
resonance and to determine its origin. The line 
widths reported in this paper are determined 
from the absolute absorption of the sample and 4 
the saturation magnetization vs T curve meas- 
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ured on single-crystal YIG. The large effects on 2h os menus 

the line width AH at room temperature due to in- . } 
homogeneous broadening by scattering from pits . £ I i 

on the surface of a YIG sphere were evaluated 0 50 100 50 200-250 300 
previously.” Consequently, we measured AH vs —— T°K 

T for a series of three spheres prepared using FIG. 1. AH vs T for single-crystal spheres of YIG 
polishing papers having mean grit sizes of 15, 5, 0.014 to 0.017 inch in diameter. Hy, oriented along 

and 0.3 microns, respectively. These measure- the [111] axis. The labels on the curves indicate the 
ments were made at 9300 Mc/sec. Figure 1 mean grit size of the final polishing paper used to pre- 
shows that the ratios of maximum line widths to pare the spheres. The lowest curve is Fig. 3 repeated 


room temperature line widths are 2.5, 4, and here for comparison. 


13, respectively, the larger ratios occurring for 
the better polished spheres. It is also clearly 
shown that the contribution to the line width due 
to surface preparation is essentially additive 
over the temperature range. We conclude there- 
fore that the low-temperature line-width maxi- 
mum is an intrinsic property of the material and 
not of the surface. Further, if perfect polishing 
were attained the line-width maximum would 
still be approximately 6 oersteds in this sample. 
Earlier investigations on the same samples* 
have shown a small but repeatable frequency de- 
pendence of line width at room temperature. The 
line width AH (full width between half maximum 
values) is 0.44 oersted at 3000 Mc/sec and 0.52 
oersted at 9300 Mc/sec along the [111] axis. , ra 
Taking these measurements over the tempera- F 3000 MC 


4H (OERSTEDS) 
ad oi ey a 











ture range, there are three essential features. 0 5 ah 1 "a 
See Fig. 2. The line width maximum at 3000 . 

Mc/sec is reduced by 40% in magnitude and is "s 

shifted from 40°K to 30°K. Also it is seen that FIG. 2. AH vs T measured at 9300 and 3000 Mc/sec. 
the room-temperature frequency dependence of Hao oriented along the [111] axis. 
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line width is related to the mechanism of the low- 
temperature line-width maximum. 

In considering the constituents that are used in 
growing* YIG, we recognized that the bulk yttrium 
oxide contains significant rare earth impurities. 
Further, the measured’ spin-lattice relaxation 
times of the rare earth ions indicated that these 
ions could contribute strongly to the observed 
temperature variation of line widths in YIG. A 
detailed theory of this effect has been prepared 
by Kittel.® 

The problem of purifying the yttrium oxide with 
respect to the rare earths was taken to Dr. Mark 
M. Woyski, of the Lindsay Chemical Division of 
the American Potash and Chemical Company. By 
using multiple exchange columns and by special 
processing against a small yield, he achieved at 
least a thousandfold reduction in the total rare 
earth content. The estimated total rare earth 
impurities are less than 0.1 part per million. 

The results for YIG grown from the special 
purity yttrium oxide are given in Fig. 3 and for 
comparison superimposed on Fig. 1. It is seen 
that the rare earth contribution’ to the AH maxi- 
mum has been reduced from approximately 6 
oersteds to slightly over 0.1 oersted, a ratio of 


0.45 
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FIG. 3. AH vs T on a 0.020-inch diameter YIG 
sphere grown from yttrium oxide containing total rare 


earth impurities estimated to be less than 0.1 ppm. 
Hj, oriented along the (111) axis. 


about 50:1. It is tempting to identify the fine 
structure with individual rare earth ions, espe- 
cially if we associate the peak at 100°K with 
Dillon and Nielsen’s terbium line. Additional 
measurements will be needed to eliminate the 
possibility of other mechanisms, such as weakly 
interacting magnetostatic modes. An over-all 
feature of the data is the general decrease in line 
width from room temperature to liquid helium 
temperatures. This is the first time, to our 
knowledge, that this trend has been observed in 
ferromagnetic resonance. It is probable that the 
line width in this sample is limited by inhomo- 
geneous broadening brought about by crystal im- 
perfections, surface effects, and the rare earth 
ions near 4.2°K; and by thermal®,® processes at 
higher temperatures. 

From the practical point of view it is interest- 
ing to note that »’’ (res) =14000 at 4.2°K. From 
the theoretical point of view it can be expected 
that measurements on single-crystal YIG grown 
from such purity will yield more precise evalua- 
tion of the fundamental mechanism involved in 
spin-wave scattering and spin-lattice relaxation 
phenomena. 

We should like to thank J. W. Nielsen for his 
advice through many discussions on the chemistry 
of garnets, and for his growing the garnets on 
which these data are based; J. F. Dillon for 
many helpful discussions; and H. E. D. Scovil 
for discussions on his measurements of rare 
earth relaxation. 
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3R. C. LeCraw and E. G. Spencer, J. Appl. Phys. 
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~ 43, W. Nielsen, J. Appl. Phys. 29, 390 (1958); J. W. 
Nielsen and E. F. Dearborn, J. Phys. Chem. Solids 
5, 202 (1958). 

5B. Bleaney and H. E. D. Scovil, Proc. Phys. Soc. 
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6c. Kittel, post-deadline paper 1959 Boston meeting 
of the American Physical Society; Phys. Rev. (to be 
published). 
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Letter [Phys. Rev. Letters 3, 30 (1959)] report on 
measurements in which deliberate amounts of specific 
rare earth ions have been added; this procedure allows 
them to measure effects due to particular ions. 
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OXIDES WHICH SHOW A METAL-TO-INSULATOR TRANSITION AT THE NEEL TEMPERATURE 


F. J. Morin 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 5, 1959) 


In a previous paper’ it was pointed out that the 
lower oxides of titanium and vanadium behave 
like metals. It was suggested that this was so 
because in the oxides at the beginning of the 3d 
series the nonbonding (f,g) orbitals extend out 
far enough to overlap and form a narrow conduc- 
tion band. It was also pointed out that polycrys- 
talline samples of Ti,O, and V,O, appear to show 
a transition from metal to insulator at the Néel 
temperature. Now, this transition has been con- 
firmed with single crystals of these oxides and 
also found to occur in the oxides VO and VO,,. 
This Letter reports electrical conductivity meas- 
urements made on polycrystalline TiO,” and on 
single crystals of Ti,O,,* VO, V,0,, and VO,, * 
as additional evidence for the existence of a con- 
duction band of f,g orbitals, and suggests a 
model to explain the conductivity transitions 
which are found. 

The single-crystal samples were 0.1 mm in 
size, too small for the usual four-point measure- 
ment, so two pressure contacts were used. 
Therefore, the absolute value of conductivity 
quoted is somewhat in doubt due to contact re- 
sistance and dimensional uncertainty. To obtain 
a Ti,O, sample of suitable uniformity, a small 
fragment was broken from a large boule. It was 
unoriented. The vanadium oxide crystals were 
easily oriented from their growth habit. The 
TiO sample was a small rectangular bar and was 
measured by the four-point method. 

Typical conductivity curves are shown in Fig. 1. 
In addition to these results, TiO was found to 
be a metal over the whole temperature range in- 
vestigated: 300°K to 1.5°K. The transition tem- 
peratures from conductivity measurements are 
summarized in Table I and compared with those 
from heat capacity (C,) and magnetic susceptibil- 
ity (x). In V,O, and VO the transitions were very 
abrupt, with no measurable time or temperature 
dependence. In VO, and Ti,O, the transitions 
were spread out over a range of temperature as 
shown in Table I, and some time dependence of 
conductivity was observed on cooling. Above the 
transition temperature, in the metallic state, all 
of the oxides showed a positive temperature coef- 
ficient of resistance typical of metals. If it is 
assumed that all the d electrons are available for 
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transport in this state, the electron mobility lies 
approximately in the range 10~* (TiO) to 10-7 
(VO) cm?/volt sec indicative of a very narrow 
conduction band. Finding a hysteresis about the 
Néel transition seems a little unusual. It should 
be recalled, however, that Cp and x measure- 
ments have been made most frequently only with 
increasing temperature. The Cy entries and the 
first two x entries in Table I are of this kind and 
compare well with the transition temperature 
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FIG. 1. Conductivity as a function of reciprocal 
temperature for the lower oxides of titanium and 
vanadium. Measurements were made along the [100] 
direction in VO, and along the C axis in V,O; and VO). 
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Table I. Conductivity transition temperatures in degrees Kelvin are compared with Néel temperature from heat 


capacity, Cy» and magnetic susceptibility, x. 











Conductivity 
Oxide Structure cooling heating Cp Xx 
TiO rock salt 
vo rock salt , 114 126 oe 
V0; corundum 153 165 169 168° 
vo, monoclinic® 340-325 335-350 3404 340° 
Ti,O, corundum 500-400 400-500 473! 390-4908 





8c. T. Anderson, J. Am. Chem. Soc. 58, 564 (1936). 


by, Wucher, J. recherches centre natl. recherche 
sci. Labs. Bellevue (Paris) 4, 254 (1952). 

Slightly distorted rutile structure; G. Andersson, 
Acta Chem. Scand. 10, 623 (1956). 

dy. Jaffray and A. Dumas, J. recherches centre 


measured by conductivity on heating. The last x 
entry was measured both with heating and cool- 
ing and shows a hysteresis.*® 

In TiO and VO, which are rock salt structures, 
each cation is connected to its twelve nearest 
neighbor cations by overlapping tg orbitals. In 
the metallic state these orbitals constitute a par- 
tially filled band. This band contains six states 
per atom and is one third full in TiO and one 
half full in VO. Let the crystal become an anti- 
ferromagnet with the spin alignment reversing 
between (111) planes of cations as found in NiO 
for example. Then electrons on one half of the 
nearest neighbors have spin up and on the other 
half have spin down. Thus the band is split into 
two sets of bands, one set for spin up and one 
set for spin down. In VO each set of bands con- 
tains three lower states per atom which are full 
and three higher states per atom which are 
empty. Thus VO is an insulator below the transi- 
tion temperature. In TiO the lower bands are 
only two thirds filled so TiO shows no transition 
to an insulator. This model is similar to one 
suggested by Slater® for MnO, in which each set 
of bands contains five lower states per atom 
which are filled and five higher states per atom 
which are empty. However, in MnO the tpg Or- 
bitals do not overlap, the band width is zero, and 
the crystal is an insulator both above and below 
the transition. 

The situation in the remaining oxides is much 
more complicated for reasons of structure and 
also because the transition in conductivity occurs 


natl. recherche sci. Labs. Bellevue (Paris) 5, 360 
(1953-54). 

©E. Hoschek and W. Klemm, Z. anorg. u. allgem. 
Chem. 242, 63 (1939). 

B. F. Naylor, J. Am. Chem. Soc. 68, 1077 (1948). 

&See reference 2. 


even though the bg orbitals are partially empty. 
A more general point of view has been suggested 
by Anderson’ and is based on the observation 
that antiferromagnets which involve indirect ex- 
change are insulators. It is suggested that there 
are two processes in competition: (1) the corre- 
lation effect which tends to localize the d elec- 
trons and which is aided by setting the cation 
spins in an antiferromagnetic alignment; (2) the 
tendency of electrons to delocalize themselves by 
spreading into a band and thereby gaining kinetic 
energy. That these two processes are fairly en- 
ergetic is indicated by the conductivity change of 
a factor 10° in VO and V,O, at the transition. A 
splitting of the conduction band by ~0.4 ev is re- 
quired to produce this amount of change in the 
carrier concentration. This may be compared 
with the energy of the Néel temperature which is 
only ~0.01 ev and which may represent the en- 
ergy difference between the two competing pro- 
cesses. This more general point of view also 
offers an alternative explanation for the differ- 
ence between TiO and VO. The electron mobility 
in TiO is ten times that in VO. This suggests 
that the conduction band is much wider in TiO so 
the electrons gain enough kinetic energy on 
spreading into it that the metallic state is always 
lower in energy than the antiferromagnetic state. 
The transition in all of the oxides can be des- 
cribed by a partially filled band which splits into 
an upper band which is empty and a set of lower 
occupied levels which are localized and have an 
antiferromagnetic spin alignment. 
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'F, J. Morin, Bell System Tech. J. 37, 1047 (1958). 

"Sample obtained from A. D. Pearson and described 
in J. Phys. Chem, Solids 5, 316 (1958). 

3Sample obtained from Linde Company. 

‘Crystals of the vanadium oxides were grown by 
H. Guggenheim by the hydrothermal process. 

5 Additional evidence was obtained that a hysteresis 
at the Néel transition is not unusual. With R. G. 
Shulman the nuclear magnetic resonance of Co®® in CoO 


was measured with slow warming and cooling about 
Ty. The transition was gradual, and on cooling ap- 
peared to occur several °K below the warming transi- 
tion. 

83. C. Slater, Handbuch der Physik (Springer-Verlag, 
Berlin, 1956), Vol. 19. 

™P. W. Anderson, Phys. Rev. (to be published) and 
private communication. 








MULTILEVEL PULSED-FIELD MASER FOR GENERATION OF HIGH FREQUENCIES* 


S. Foner, L. R. Momo, and A. Mayert 
Lincoln Laboratory, Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received June 1, 1959) 


A proposed’? multilevel pulsed-field ruby 
maser has been successfully operated at 4.2°K 
as an oscillator at both 12.61 kMc/sec and 19.15 
kMc/sec with a pumping frequency of 12.61 
kMc/sec. The results demonstrate the feasibil- 
ity of generating or amplifying very high fre- 
quencies on a pulsed basis by conversion of mag- 
netic field energy, supplied to an inverted spin 
distribution, into coherent radiation. The prin- 
ciples of operation and characteristics of the 
maser, and field dependence of spin-lattice re- 
laxation time are discussed. 

The pertinent transitions and the energy level 
diagram of Cr** in Al,O, for applied fields per- 
pendicular to the c axis® are shown in Fig. 1. 
For simplicity, the AM=2 transition between 
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FIG. 1. Energy level diagram for Cr** in Al,O; with 
applied field perpendicular to C axis. Operating points 


for pumping at v3=1,+v,, and emission at v; and », are 
indicated. 
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levels 1 and 3 is saturated* at a pump frequency 
Vs =V,+¥,, Corresponding to a fixed dc bias field, 
H,, above any field where cross relaxation ef- 
fects’ may occur, such as at »,. The cw three- 
level maser® would operate at y,.<<v,, when H 
=H,, with an inverted spin population between 
levels 2 and1. If H is now increased to H; ina 
time short compared to the spin-lattice relaxa- 
tion time, oscillation or amplification should be 
obtainable at any frequency v;>v, for which suit- 
able maser cavity resonances occur. A gain of 
energy in the ratio of v;/ v, over the continuous 
three-level maser is obtained at the expense of 
energy extracted from the pulsed magnetic field. 
Wavelengths of about 1 mm would be generated 
at a field of 100 kilo-oersteds. 

The maser reflection cavity, a solid rectangular 
parallelpiped of ruby,’ was designed with the two 
lowest modes (TE,,, and TE,,,) occurring at v, 
=12.61 kMc/sec and y,=19.15 kMc/sec, respect- 
ively, with undercoupled loaded cavity Q’s of 
about 2000 and 4800, respectively. The pump 
power and radiated energy were transmitted to 
and from the single cavity coupling iris by a tap- 
ered dielectric transition section attached to K «, 
band waveguide. Both the dielectric transition 
and cavity were coated with a 0.001-in. silver 
plating in order to assure uniform penetration of 
the pulsed magnetic field throughout the cavity 
volume. The pulsed field, H,, was generated in 
a solenoid (13 cm length, 7 cm i.d.) by discharge 
of a 2000-:f variable voltage capacitor bank.® A 
peak field of 9.4 kilo-oersteds was obtained at 


1000 volts with a half-period of about 3 millisec- 
onds [see Figs. 2(A), 2(B)]. The maser cavity 
assembly and surrounding glass Dewars were in- 
serted in the solenoid, and the entire assembly 
was fixed between the pole faces of an electro- 
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FIG. 2. Oscillograms of pulsed-field maser operation 
at 4.2°K. (A) Pulsed magnetic field versus time-—hori- 
zontal time scale 500usec/cm. (B) Pulsed magnetic 
field versus time—horizontal time scale 100usec/cm. 
(C) Envelope detector output showing 19. 15-kMc/sec 
emission of pulsed-field maser—horizontal scale 100 
psec/cm. Emission is obtained during both rise and 
decay of field pulse. (D) 19.15-kMc/sec emission— 
horizontal time scale expanded to 20usec/cm. An 
intense burst of radiation occurs at the onset of emis- 
sion, accounting for the short blank portion of the oscil- 
logram. 





magnet. The field orientations were: H,¢ per- 
pendicular to the c axis, and H, perpendicular to 
both the c axis and to H,. The mutually orthog- 
onal orientation of H,, H,, and the c axis assured 
that the resultant H was always greater than H,, 
and furthermore, that H always remained per- 
pendicular to the c axis. The field dependence of 
the energy levels therefore remained as indicated 
in Fig. 1. 

The pumping power was gated off within 50 
usec of the start of the pulsed field rise, and 
oscillations were observed at either v, or 
during the traversal of H, through the appropriate 
resonance. Within the accuracy of the pulsed 
field measurement, emission occurred as pre- 
dicted from Fig. 1. When the pump power was 
not turned off, population redistribution between 





levels 4-3 and 3-2 was observed (see Fig. 1). 
Typical oscillation characteristics at 4.2°K for 
the y, mode are shown in Figs. 2(C) and 2(D). 

The expanded scale of Fig. 2(D) shows the pul- 
sating character® of the oscillations, and the 
sensitive dependence of the emission pulses on 
the particular magnetic field sweep rate. The 
number of pulses and total energy output could 

be controlled by varying the pulsed-field sweep 
rate—the more rapid traversals resulted in 

lower peak output, and for a sufficiently fast rise 
time, no oscillations were observed. In this way 
it was possible to avoid emission at », completely, 
and yet obtain v, emission at the peak of the 
pulsed field. When only a small portion of the 
inverted spin system radiated at y,, further 
emission at y, was observed during the field de- 
cay. The peak power output at vy, was greater 
than 300 uw, and the integrated output energy for 
a typical case, e.g., Fig. 2(D), was greater than 
0.04 erg. In both cases measurements were made 
of energy arriving at the crystal detector. Cal- 
culations show that a total energy of about 2 ergs 
would be available in the cavity if compiete sa- 
turation between levels 1-3 were obtained, and 

if cavity losses and the spin-lattice relaxation 
were neglected. A considerable increase in ef- 
ficiency is expected from an optimized micro- 
wave system. 

Operation of this particular pulsed field maser 
would be limited to low magnetic fields if the 
spin-lattice time varied inversely as some large 
power of the magnetic field.*° Measurements of 
spin-lattice relaxation time in ruby" at 4.2°K 
combined with the present results at v, indicate 
that such a dependence is not yet dominant. Re- 
sults of experiments as a function of temperature 
and field will be presented elsewhere. 

An alternate scheme of multiple pumping’»” 
would circumvent most of the spin-lattice time 
field dependency, and would minimize require- 
ments on the pulsed field system. Spin inversion 
with 4M =1 transitions between levels 1-2, 2-3, 
and 3-4 in that order may be attained with a re- 
latively small pulsed magnetic field at a reason- 
able pump frequency v;. With appropriate orien- 
tations of H, and H,, and with a suitable maser 
material of large zero-field splitting, emission 
at a AM =2 or 4M =3 transition with v;~3v; would 
then permit generation and/or detection and am- 
plification at millimeter wavelengths on a pulsed 
basis with a rather simple arrangement. For 
instance, 10° spins radiating at 300 kMc/sec 
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with a repetition rate of only 10 per second would 
yield an average power of 20 milliwatts—a power 
somewhat above that presently available (con- 
siderably less than 1 pw). If multiple pumping 
could be employed, pulsed fields of the order of 
a kilo-oersted would be sufficient. 

We are indebted to B. Feldman and J. J. Kelley 
for valuable assistance, and.to Dr. R. H. Kings- 
ton, Dr. H. J. Zeiger, and Dr. S. H. Autler for 
informative discussions. 
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Lincoln Laboratory, a center for research operated by 
Massachusetts Institute of Technology with the joint 
support of the U. S. Army, Navy, and Air Force. 

tPresent address: Convair, San Diego, California. 

1g, Foner, Colloque International de Magnétisme de 
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'The energy level diagram is shown for D<0 as indi- 
cated by recent measurements of Geusic, Peter, and 
Schulz-du Bois, Bell System Tech. J. 38, 291 (1959). 

‘In principle AM=3 pumping may be feasible, and/or * 
spin inversion by a 180 degree pulse or adiabatic fast 
passage may be employed for higher efficiency. 

5Bloembergen, Shapiro, Pershan, and Artman, Phys, 
Rev. 113, 445 (1959). 

‘N. Bloembergen, Phys. Rev. 104, 324 (1956). 

"Pink ruby (nominally 0.1% cr) supplied by Linde 
Air Products Company. 

8s, Foner and H. H. Kolm, Rev. Sci. Instr. 28, 799 
(1957). | 

*Makhov, Kikuchi, Lambe, and Terhune, Phys. Rev. 
109, 1399 (1958). 
“TWA, H. Cooke, Reports on Progress in Physics (The 
Physical Society, London, 1950), Vol. 13, p. 276. 
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SMEKAL-RAMAN TYPE MODIFIED X-RAY SCATTERING 


K. Das Gupta* 
Department of Physics and Astronomy, Ohio State University, Columbus, Ohio 
(Received June 4, 1959) 


The spectroscopic analysis of x-ray beams, 
scattered from some low atomic number ele- 
ments, viz., carbon, beryllium, and aluminium, 
made with the help of the two-crystal spectrom- 
eter, multicrystal spectrograph, and focusing 
spectrograph by DuMond,’ Bearden,” Gingrich, * 
Ehrenberg,‘ and a few others, verified the exist- 
ence of Rayleigh (unchanged frequency) and Comp- 
ton scattering but could not detect the existence 
of the phenomenon of the partial absorption of 
photons of Smekal- Raman type. 

With scatterer carbon and a monochromatic 
incident x-ray beam of energy hv, the frequency 
v’ of the modified line of Smekal- Raman type can 
be calculated from the equation hv’ =hv - 5E, 
where 65E in a typical process will correspond to 


the energy difference between the initial (1s72s?2p") 


and the final (1s'2s?2p*) state of the carbon atom. 
The value of 5E for the element carbon as ob- 
tained from the x-ray K-absorption data is 282.2 
ev and the wavelengths of the Smekal- Raman type 
modified lines due to incident copper and molyb- 
denum radiations, viz. (1) Cu Ka, =1537.4, (2) 
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Cu Ka, =1541.2, (3) Cu Kp, =1389.4, (4) Mo Ka, 
=707.8, (5) Mo Ka,=712.1, and (6) Mo KB, = 630.98 
xu (Siegbahn x units), will be (t) 1593.3, (2) 1597.4, 
(3) 1434.9, (4) 719.5, (5) 723.9, and (6) 640.2 xu, 
respectively. 

With the help of a bent mica (d=2.554 A) Cau- ' 
chois spectrograph (diameter of the focal circle 
=9 inches) and also a bent quartz (d=1.178 A) 
spectrograph (diameter of the focal circle = 5.5 
inches) and using Ka and Kf radiations of copper 
and molybdenum from the sealed-off x-ray tubes, 
with graphite (reactor quality), polystyrene (for 
carbon), and beryllium as scattering specimens, 
the modified lines of Smekal- Raman type have 
been observed in calculated positions using both 
the spectrographs for each of the scattering spe- 
cimens. Figures 1 and 2 are the tracings of the 
original microphotometer records, clearly show- 
ing (1) Rayleigh lines, (2) Compton bands, and (3) 
Smekal- Raman type of modified lines in x-ray 
scattering. 

The sealed-off x-ray tubes (Matchlett) were run 
at 12 to 15 ma at 40 kv and the period of exposure 
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FIG. 1. Spectroscopic analysis of Cu target radia- 
tion scattered at ~90° by polystyrene with the help of 
the bent mica Cauchois spectrograph. The photometer 
record clearly shows (1) Rayleigh lines, viz., CuKa, » 
and KB,; (2) Compton bands due to Ka, » and KB ;; (3) 
modified Smekal-Raman lines due to Cu Ka, , and Kp,. 
Dispersion on the film =21.6 xu/mm. The Ka, 2 lines 
are not resolved. 


varied between 100 and 150 hours. The target ra- 
diations from both copper and molybdenum were 
spectroscopically analyzed by giving prolonged 
exposures, and, only in the case of the copper 
target tube, Ka lines of nickel and La of tungsten 
appeared in the spectrum. The molybdenum tar- 
get was, however, free from any such surface 
contamination of nickel and tungsten. The scat- 
tered radiation from graphite and polystyrene 

has been studied at scattering angles of approxi- 
mately 90° (Fig. 1) and 97° (Fig. 2). The effec- 
tive width of the bent crystal was about 6 mm and 
the slit in between the scatterer and the bent 
crystal is made of pure copper which allowed a 
beam of 6 mmx4 mm. The approximate distance 
between the scatterer and the focal spot of the 
target is about 1.5 inches. The effective irra- 
diated area of the piece of graphite of 5 mm 
thickness is 1 cmx6 mm. The piece of graphite 
was kept at 45° to the direction of the incident 
beam. The scattering experiments with carbon 
and with beryllium have been done by keeping the 
Spectrograph setting exactly the same, only the 
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FIG. 2. Spectroscopic analysis of Mo target radia- 
tion scattered at ~97° by graphite, with the bent quartz 
Cauchois spectrograph. The photometer record shows 
clearly (1) Rayleigh lines, viz., Mo Ka, and Ka,, (2) 
strong Compton band, and (3) the modified Smekal- 
Raman line due to Mo Kay; that due to Ka, overlaps 
with the short-wavelength end of the Compton band and 
is just visible in the original negative. The dispersion 
as measured in the original negative = 16.1 xu/mm. 


piece of graphite is replaced by a similar piece 
of beryllium. The observed modified lines due 
to carbon appeared at 282-ev gap and that due to 
beryllium at 116-ev gap on the low-energy side 
of the primary Ka of copper or of molybdenum. 
The preliminary experiments made with the sam- 
ples of extra pure magnesium and aluminium as 
scatterers, at different setting angles of the 
spectrograph, confirmed beyond doubt that the 
modified lines due to graphite, polystyrene, and 
beryllium appearing at the calculated positions 
cannot be due to any impurities either in the tar- 
get or in the sample. 

Although modified lines have been obtained with 
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different scatterers, viz., Be, C (graphite, poly- 
styrene, apiezon oil), Mg, Al, Cu, and Fe at 
calculated positions, most of them are, however, 
extremely weak in intensity and it has not been 
possible to take their microphotometer records. 
It is hoped that by using the high light gathering 
power. spectrograph, high contrast films, strong 
X-ray sources, and selecting suitable angles of 
scattering so as to have the minimum background 
due to the tail of the Compton band or the higher 
order Bragg reflections due to the white part of 


the target radiation, it will be possible to study, 
from this newly observed modified scattering, 
the term values of the elements and also the en- 
ergy levels of the solid state. 





* 

Permanent address: Department of Physics, Uni- 
versity of Calcutta, Calcutta, India. 
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GROUND STATE A-DOUBLING TRANSITIONS OF OH RADICAL* 


G. Ehrenstein and C. H. Townes 
Department of Physics, Columbia University, New York, New York 


and 


M. J. Stevenson 
Research Laboratory, International Business Machines Corporation, Poughkeepsie, New York 
(Received June 4, 1959) 


One of the molecules whose presence in inter- 
stellar space may be detectable by means of its 
radio-frequency spectrum is the OH radical.’ 
Attempts at observing it with radio telescopes 
have been unsuccessful thus far. To make future 
searches more fruitful, frequencies of the appro- 
priate absorption lines were measured in the 
laboratory with a Zeeman-modulated spectrome- 
ter. 

For OH, at temperatures prevailing in inter- 
stellar space, only the lowest rotational state 
(J = 3) of the *Iy, electronic state will be appre- 
ciably populated. This lowest energy level is 
split by A doubling into two levels 1666.4 Mc/sec 
apart. Each of the two A-doublet levels is split 
further by hyperfine interactions with the hydro- 
gen nucleus, so that the spectrum consists of 
four absorption lines. The two stronger transi- 
tions, having no change in the total angular momen- 
tum F(AF=0), were measured and their frequen- 
cies are given in Table I together with frequencies 
calculated from expressions given by Dousmanis, 
Sanders, and Townes.* According to theory, the 
intensity ratio between the F =2— F =2 and 
F=1-+F=1 transitions should be 9:5, and this 
was verified experimentally. Substantially the 
same ratio should prevail in interstellar space. 
The transitions with AF =+1 are expected to be 
about ten times weaker and have not been observed 
as yet. 


40 


Table I. Observed and calculated frequencies for A- 
doubling transitions in J =} rotational state of Mls,. 
electronic state of OH radical. 











Experimental Calculated 
Hyperfine frequency frequency 
transition (Mc/sec) (Mc/sec) 
F=2-F=2 1667.34 + 0.03 1666.5 
F=1-F=1 1665. 46 + 0.10 1664.6 








OH radicals were produced outside the absorp- 
tion cell by an electrodeless radio-frequency 
discharge.* The microwave spectrometer was 
conventional in design except for the Zeeman- 
modulation technique. The absorption cell was a 
glass-lined, coaxial waveguide about 23 feet long 
and 14 inches in diameter. This coaxial cell was 
used both to propagate low-frequency microwave 
power and to provide Zeeman modulation. The 
modulating magnetic field was obtained by pass- 
ing a large current through the center rod and back 
along the outer conductor. This produces a mag- 
netic field between the two conductors which ideally 
does not extend outside the waveguide. The very 
small leakage fields which were obtained produced 
only very low pickup in the detecting circuits, 
whereas the pickup was quite troublesome when 
a conventional solenoid was wound around the 
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waveguide and used for Zeeman modulation. The 
modulating current consists of a 100-kc/sec sine 
wave current superimposed on a direct current, 
in such ratio that the minimum of the sin wave 
corresponds to zero current. With a modulator 
whose maximum output was 25 peak-to-peak 
amperes, adequate modulation fields were ob- 
tained for both ?f},,. and *I,,. states of OH. This 
technique for Zeeman modulation can also be 
applied to a standard Stark cell with the Stark 
septum acting as the center conductor. 

The authors wish to thank the Sperry Gyroscope 
Company and Professor Yardley Beers for lend- 


ing us klystron tubes used in these measurements. 
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EVIDENCE FOR SHELL-MODEL STATES IN Sc*T 


C. M. Class, R. H. Davis,* and J. H. Johnson 
The Rice Institute, Houston, Texas 
(Received April 6, 1959) 


The excited states of the mirror nuclei Ca** and 
Sc“ are of considerable interest because of the 
possibilities they present for testing the predic- 
tions of the shell model under quite favorable 
circumstances. According to this model, the 
ground states of these nuclei are described by 
single particles in the f.,. shell outside a Ca® 
core described by closed (1s)*(1p)'*(1d, 2s)** 
shells. Single-particle excited states are pre- 
dicted approximately in the order py, Pyo, fs, 
8yo, Igo, and gx. In addition, states of more 
complex character are expected due to the ex- 
citation of nucleons from the Ca® core.’ These 
states should be bound in the case of Ca* for 
which the separation energy of a neutron is 8.37 
Mev, while for Sc*! they are expected to be un- 
bound since the separation energy for a proton 
is only 1.63 Mev.” Whether or not the difference 
in the binding of the states has an important effect 
on their ordering is unclear. 

The experimental evidence supporting the ex- 
istence of the single-particle states is rather 
scant and comes for the most part from the 
Ca*(d, p)Ca* reaction.* By application of strip- 
ping theory to the angular distributions of the 
resolved proton groups, it is possible to assign 
the orbital angular momenta of the captured neu- 
trons and by a comparison of the relative cross 
sections at the maxima of the angular distribu- 
tions, a distinction can be made between states of 
mainly single-particle character and those of 
more complex character. On this basis, the p,. 


and p,. levels have been tentatively assigned at 
energies of 1.947 and 3.950 Mev, respectively.* 
Seven other levels between 2 and 4 Mev have been 
assigned / values (J <2), but none of these are 
formed with sufficient intensity to suggest possi- 
ble single-particle character. In addition, thir- 
teen states in this interval have been located 
which show no stripping maxima. The data avail- 
able on the states of Sc** are less complete. 
States formed by /=1 protons in the Ca®(d, n)Sc* 
reaction have been reported at excitation energies 
of 1.86 and 2.07 Mev,° one of which by analogy 
with Ca“ is probably the p,,. single-particle state. 
Levels at 2.25 and 3.44 Mev are known from (9, y) 
measurements,° and at energies above 6 Mev a 
number of states have been located by the study 
of the Ca*(p, p’e*+)Ca® reaction.’ 

In the work reported here, additions to the body 
of information available on Sc* have been made 
by studying the elastic scattering of protons from 
Ca*.® As is well known, when elastic scattering 
measurements are made with good resolution, 
one is able not only to locate excited states of the 
compound system but also to make spin and pari- 
ty assignments for the states by an analysis of 
the resonance shapes. This analysis is much 
simplified when, as in the case at hand, the tar- 
get nucleus has zero spin.® Another feature of 
interest in this case results from the low excita- 
tion energy with which Sc*! can be formed, which 
enables one to examine much of the energy in- 
terval relevant to shell-model considerations. 
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The details of the experimental method and 
the complete analysis of the results will be pub- 
lished elsewhere.” To summarize briefly the 
essential points here, the measurements were 
carried out at bombarding energies between 1.2 
and 5.5 Mev, corresponding to excitation ener - 
gies in Sc** between 2.8 and 7 Mev. The region 
below an excitation energy of about 2.8 Mev was 
inaccessible because of the weak penetration of 
the Coulomb barrier by low-energy protons. The 
protons, scattered from thin (~2 kev) Ca targets 
evaporated onto 100-yg/cm? thick carbon foils, 
were detected at 90, 125, and 150 degrees with 
respect to the beam direction by a CsI crystal 
spectrometer which satisfactorily resolved the 
elastic groups from calcium and carbon. The 
shapes of the resonances were analyzed by appli- 
cation of dispersion theory, aided by the use of 
an electronic computer." Essentially unambig- 
uous assignments of the 7 values of the protons 
forming forty-two states were made as a result 
of this analysis, and it was possible to make j- 
value assignments to about twenty states with 
considerable certainty. In addition to /(j) and 

nm, the widths of the levels were also obtained 
from the analysis after suitable corrections were 
made for the effects of finite target thickness and 
energy spread of the beam. The I°’s so obtained 
needed no correction for effects of competing 
reactions in this case since the (p, p) process is 
the only particle reaction which occurs until a 
bombarding energy of about 5 Mev is reached, 
when the (p, p’) reaction develops. The (p, y) 
partial width is assumed negligible. From know- 
ledge of the TI and the 7 value of a state, the re- 
duced width 7*, given by )“=I'/2kRP), is readily 
obtained, where P, is the penetrability of a proton 
having angular momentum /7, k is its wave num- 
ber, and RA is the nuclear radius which in this 
case was taken to be 5.13 fermis. The magnitude 
of this quantity relative to the Wigner limit for 
this quantity gives a measure of the single -parti- 
cle character of the state involved. In Fig. 1 the 
values of }*/y.? for states of given / (or j) are 
plotted as a function of the excitation energy. 

The results shown in Fig. 1 display several 
noteworthy features: 

(1) The marked peaking of 7*/y,,” for different 
1(j) values suggests the influence at these ener- 
gies of the single-particle states expected from 
the shell model although these states appear in 
a somewhat different order than predicted.’ The 
discrepancy of order is not surprising and calls 
for the choice of a more suitable potential for 
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FIG. 1. The reduced widths of the states observed 
in elastic scattering of protons, expressed as a frac- 
tion of the Wigner single-particle limit, are shown 
as a function of the excitation energy in Sc’. The 
lines connect the data points for states formed by pro- 
tons of a fixed 7 value. 


determining the positions of these states. It 
might also be remarked that s states, for which 
the efficiency of detection is greatest due to the 
favorable penetrability factors, are relatively 
few and that their reduced widths are significantly 
smaller than for those states expected from the 
shell model. This would seem to imply that the 

s states in this region of excitation involve a com- 
plex mode of excitation of the Ca® core. 

(2) While there are too few states contributing 
to make a strong argument, it appears that the 
envelope of the reduced widths of states formed 
by protons of a given /(j) value shows a character 
which is illustrative at relatively low energies of 
excitation of the giant resonances discussed in 
the complex potential model” and in a more de- 
tailed fashion in the intermediate-coupling model 
of Lane, Thomas, and Wigner.’* The occurrence 
of these maxima may also be viewed as supporting 
the interpretation of the strong maxima found in 
























ee ee ee 


~~ SS Shhh 




















VoLUME 3, NUMBER | 


PHYSICAL REVIEW LETTERS 








Jury 1, 1959 

















Table I. Experimental single-particle level spacings. 
AE (Mev) AE (Mev) 

Levels (present work) (reference 3) 

Ifsj2- Ure 6.4 =4 

2P1/2- 2P3/2 2.2 ~1.8 

Psa- Us 4.5 "s 

189/2- 2P1/2 1.5 act 

2ds /2- 2P1/2 1.9 





the spectra of protons from (d, p) reactions for 
this region of nuclei as due to the influence of 
single-particle states on the level structure.** 

(3) If one assumes that the groups of states 
according to /(j) value are in fact identifiable 
with the appropriate single-particle states, then 
the mean energy spacings of these states are of 
interest in connection with theories of the origin 
and magnitude of the spin-orbit force. The 
spacings observed in our measurements are sum- 
marized in Table I, and for comparison previous 
estimates of these spacings are included when 
available. 
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ENERGY SPECTRUM OF PHOTONEUTRONS FROM OXYGEN 


Carmelo Milone 
Istituto di Fisica dell’ Universita, Centro Siciliano di Fisica Nucleare, Catania, Italia 
(Received March 30, 1959; revised manuscript received June 2, 1959) 


The energy spectrum of the photoneutrons from 
oxygen has been studied by irradiation of a water 
target with a 31-Mev collimated bremsstrahlung 
beam from the Brown Boveri betatron of Torino 
University. 

Photoneutrons emitted at about 90° with the y- 
ray beam have been recorded by means of the 
proton recoil tracks in Ilford L, plates 400 yp 
thick. A water wall 60 cm thick screened the 
plates against the spurious neutrons coming di- 
rectly from the betatron, as in a previous work. 

Plates were scanned and proton recoil tracks 
were selected following the method used previ- 
ously.? 


ee 


The neutron spectrum is deduced from the pro- 
ton recoil spectrum taking into account the cross 
section for (m,p) collisions in the hydrogen of the 
emulsion. The relation E, =E,/cos*@ was intro- 
duced only for 6>10°. Small corrections are due 
to the probability of escape of the tracks from 
the emulsion’ and to the absorption and scatter - 
ing of neutrons in the target. The experimental 
energy spectrum is given in Fig. 1(a). 

Under the restrictive assumption that only the 
direct O(y,) process is effective and that the 
residual O** nucleus is left in the ground state, 
the expected neutron spectrum may easily be in- 
ferred from the O(y,m) cross section,® taking 
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FIG. 1. (a) Experimental photoneutron spectrum. 
N=neutrons per AE, =0.5 Mev (arbitrary units). (b) 
and (c) Spectra calculated from the O(y,”) cross sec- 
tions. 


account of the bremsstrahlung spectrum. The 
behavior of the neutron spectrum thus inferred is 
given in Fig. 1(b). 

Although the above assumption is very restric- 
tive and the O(y,) cross section’ is available 
only in large energy steps, a rough agreement 
may be observed in the high-energy region be- 
tween the experimental and calculated spectra. 
This gives evidence of a large contribution of di- 
rect process in the O(y,m) reaction. On the other 
hand, for all elements previously investigated,* 





evaporation accounts for the greater part of the 
emitted neutrons. A more detailed comparison 
between the experimental and the calculated 
spectra in the high-energy region may be made 
up to Ey =25.1 Mev taking into account Spicer’s 
cross section’ calculated in 0.5-Mev steps. The 
spectrum (c) in Fig. 1 inferred from Spicer’s 
cross section in the giant resonance region shows 
a fine structure as the experimental spectrum. 

In the low-energy region a contribution of the 
evaporative process may account for the differ- 
ence between the experimental and calculated 
spectra. 

The neutron peak around 2.5-3 Mev may be due 
to neutrons leaving O** in excited states. The 
two neutrons peaks A and B are to be attributed 
to photons of energy E., around 22 and 24 Mev, 
respectively. Indeed, taking into account (i) the 
O**(y,) threshold energy (E;=15.6 Mev), (éi) the 
energy of the first excited level in O*(E, =5.3 
Mev), and (iii) the behavior of the O(y,) cross 
section beyond 27 Mev,° it appears that after the 
emission of the neutron groups A and B the O** 
must be left in its ground state. In this particular 
case E. =15.6+ (16/15)E,. 

It gives this writer great pleasure to express 
his thanks to Professor R. Ricamo and Professor 
G. Wataghin for their generous help and to Pro- 
fessor G. Cortini for valuable discussion of the 
program of the research. Grateful thanks are 
due to the Catania and Torino laboratory col- 
leagues for their friendly collaboration. 
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SHELL EFFECTS AT HIGH EXCITATION ENERGIES* 


H. W. Newson and M. M. Duncan 
Duke University, Durham, North Carolina 
(Received May 25, 1959) 


It has been known for some time that a plot of 
the level spacings of rather highly excited com- 
pound nuclei (as measured by s-wave neutron 
resonances) fluctuates rapidly as a function of 
mass number.’»? The fluctuations are evidently 
connected in some manner with the shell struc- 
ture of nuclei since prominent maxima of the 
level spacings appear at neutron numbers of 50, 
82, and 126. Rosenzweig*® has calculated relative 
shell-model level densities. We are presenting 
similar calculations based on an even simpler 
model. 

We shall attempt to account for level spacings 
of even-A compound nuclei by considering the 
most significant mode of excitation of compound 
nuclei to be excitations of 0* pairs of like nu- 
cleons among the levels of a spherical shell- 
model scheme. For our calculations we assume 
that there are a minimum number of unpaired 
particles (i.e., two in odd-odd or even-even com- 
pound nuclei with nonzero angular momentum). 
One can make a simple count of the number of 
possible different ways of arranging the two un- 
paired nucleons in a shell scheme and still ob- 
tain a J” consistent with the spin and parity of 
the target nucleus and an incident s-wave neu- 
tron. This count, which has been carried out 
only for odd-odd compound nuclei between N = 50 
and 82, is roughly proportional to 2J+1 (as long 
as J is not too large), in qualitative agreement 
with the predictions of a simple Fermi gas 
model.* We therefore define a level spacing 
Dy =2(21+1)Dop5=1/p. which should be approxi- 
mately independent of all but pair excitations of 
nucleons; J is the spin of the target nucleus and 
Dp g is the observed spacing including all 
S-wave channels. 

If one assumes, for the moment, that a major 
shell (only partially filled) is isolated, it is easy 
to calculate the number of possible ways of ar- 
ranging indistinguishable pairs in the shell as a 
function of an excitation energy which is meas- 
ured in units of S, an average spacing for pairs 
between subshells of a major shell. These sim- 
ple pair excitations alone are not enough to ac- 
count for the observed level densities. However, 
there are two (complex pair) modes of excitation 
available: Because of the large separation energy 


(about 2 Mev per nucleon’) between major shells, 
one and only one pair can be excited into a higher 
major shell. As an alternative, one or possibly 
more pairs may be brought up from the filled 
subshells below the unfilled major shell, as will 
be discussed later. These complex pair excita- 
tions can occur in a number of different ways 
(index i) which leaves an energy E’-E ; for the 
simple pair excitations. It is clearly possible to 
calculate the number of possible configurations 
at a given nuclear excitation energy, E’ above 
the lowest energy state with J", by sharing the 
available energy in all possible ways between 
simple and complex pair excitations. If the dis- 
tribution of the number of simple pair excitations, 
r(E’-E;), is calculated from a shell model, and 
ry is assumed independent of i, then the total con- 
figuration density at E’ is the sum over all pos- 
sible values of E;. If, for simplicity, one as- 
sumes E;=nS’, where n=0, 1, ..., then 





, 


ne 
1/Dy = po(E’) = r(E’ - nS’) = r(E’ -E)dE/S’. 
n=0 0 

(1) 


The maximum energy, M, which can be absorbed 
by simple pair excitation in the partially filled 
shells, is roughly equal to or less than the ex- 
citation energy, B, of the compound nuclei when 
35<N<90 and between 120 and 130 if the pair 
spacing, S’, is about 0.5 Mev. For these cases 
Eq. (1) becomes 


Po(B2M) =1/Dy=n7n,/S’. (2) 





nz and my are the number of arrangements of the 
indistinguishable proton pairs and neutron pairs, 
respectively, within their major unfilled shells. 
The results of these calculations based on Eq. 
(2) are shown in Fig. 1(a) for odd-odd compound 
nuclei. The predicted spacings near the rare 
earths are much too low (dashed curve) due to 
the fact that M>> B for these compound nuclei so 
that Eq. (2) breaks down very badly. The higher 
calculated points in this region are based on an 
approximate correction which assumes a trian- 
gular distribution for y, and that B=6 Mev, S’ 
=50 kev, andS=500 kev. In Fig. 1(b) the predic- 
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tions of Eq. (2) are compared to some of the 
even-even compound nuclei. It is seen that the 
general trend of the experimental results and in 
some cases their detailed behavior are correctly 
predicted. 

Sample counts have shown that the general trend 
of the predicted spacings is unaffected by as- 
































10” ; - -——— - - 
10°F | 4 
2228 Odd-Odd #Even- Even 3 

fo 
—_ di 

ay ft 4 
ro To 4 

1 of a [° #50 | 
VE lo! - 4 9: 
E “id t| « 7 

‘4 Vv Y + | = 
Ta ‘ ° | \ — 3 
— 3 . 4 re) = 
= F ° ¢ “oR % 4 
~— a 3° ! 4 
o + ' re) e- 
val BY ; 
/ ’ . 

a 4 

Ok rl H = te 
3 ’ ; > 

5 ; q 

: ‘ bs + _ 

f ‘ : ba 

‘ io d 

1E V/ 3 3 

: i" 3 

. 50 82 126+ 82 - 

: (a) T (b) - 

0 1 i i i i j 
30 50 70 90 0 50 70 #490 


FIG. 1. Calculated and measured values of D, are 
plotted against the number of neutrons N in the target 
nucleus for even-even and odd-odd compound nuclei. 
Experimental points are shown as circles. Calculated 
points are shown as small solid dots connected by 
straight solid lines except that the predictions of Eq. 
(2), where it is invalid, are connected by dashed 
lines. Experimental errors are not indicated. The 
best known spacings (in general the smallest) are un- 
certain by + 20% or more; the most uncertain points 
are those near N=50 where the errors are of the 
order of a factor of two or more in either direction. 
The same value of S’ (50 kev) has been used in all 
calculations. 
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suming four rather than two unpaired nucleons. 
In any case, this number must be small since 
the pairing energies® [about 23(2j+1)A~* Mev] 
limit it. The effect of this restriction on the 
number of possible configurations is now being 
investigated in detail by Rosenzweig.” However, 
sample calculations (based on the simplest cases) 
indicate that a count of all the energetically pos- 
sible “broken pair” configurations predicts too 
low a spacing and, hence, we must assume that 
the stability of most of these configurations does 
not permit them to be observed as sharp reso- 
nances. Consequently, we are unable (in general) 
to estimate the value of S’ from the shell model, 
since some of the observed resonances may be 
due to an unknown fraction of the broken pair con- 
figurations. In addition, we may also observe 
configurations arising from the excitation of 
pairs out of the filled subshells below the shell 
of interest. Since the conventional shell-model 
diagram is only a semiempirical summary of in- 
formation about excitations above the Fermi 
level, we cannot predict the energetics of this 
last mode of excitation. It is rather surprising 
that a single value of S’ (50 kev) yields as good 
agreement between calculated and measured 
values as that in Fig. 1. 

We have profited greatly from discussion with 
Dr. M. F. L. Pryce and Dr. E. Merzbacher. We 
are also indebted to Dr. J. A. Harvey for his 
valuable estimates of spacings from the experi- 
mental data® and to Miss E. Winton for help with 
the calculations. 


*Work supported in part by the U. S. Atomic Energy 
Commission. 
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TRIPLE ISOMERISM IN Ir'*t 


Gertrude Scharff-Goldhaber and Michael McKeown 
Brookhaven National Laboratory, Upton, New York 
(Received June 12, 1959) 


About eight years ago we started a systematic 
program to search for long-lived isomers pro- 
duced by pile neutron irradiation. Samples of 
most elements, weighing ~50 mg, were placed 
for one month inside the Brookhaven reactor 
and were then allowed to “cool.” An analysis of 
the radiations from these samples carried out 
two years after the irradiation gave evidence only 
of well-known activities. A renewed study under- 
taken in the beginning of this year revealed a new 
activity induced in the Ir sample: Although the 
well-known 74.37-day ground-state activity’ 

Ir'*& should have decayed to ~1 disintegration 
per sec, according to our computation, we ob- 
served a gamma spectrum, identical with that 

of Ir'*& (see Fig. 1), but corresponding to a 
source strength of 2820 + 140 disintegrations 

per sec. As Ir'™& is known to have a 1.45-min 
isomer,’ decaying mainly by a 58-kev isomeric 
(E3) transition to the 74-day ground state, we 
are forced to conclude that the activity observed 
stems from a third, longer lived isomer of Ir'™. 
This then constitutes a case of a nuclear isomeric 
triplet. Only one other case of triple isomerism 
is known, again in an odd-odd nucleus: Sb!*,? 
No definite spin assignments were possible so 
far for the metastable states in Sb'** because of 
the very low energies of the isomeric transi- 
tions. Hence, it seemed of particular interest 
to study the mode of decay of the new long-lived 
isomer of Ir!®, 

We shall distinguish the two metastable states 
of Ir'* by naming the 1.45-min isomer Ir**"""1 
and the new long-lived isomer Ir'”2. Measure- 
ments of the source intensity of Ir'””’2 through 
several weeks showed no decay within the statisti- 
cal limits of error, allowing us to state a lower 
limit T,,>5 yr. 

The identity of the gamma-ray spectrum 
accompanying the decay of Ir'*”"2 with that of 
Ir'*& indicates that Ir'*”"2 decays by an isomeric 
transition. In order to determine the multipole 
order of this transition, we first studied the L 
X-rays emitted from a thin source of Ir’?”2, 
using a 3-mm_-thick Nal(T1) crystal as detector. 
We found that the ratio J L x-rays/sie key from 
this sample was ~16 times larger than that from 
an equally thin Ir’ source, studied in the same 
geometry. There was no measurable difference, 


however, between the Ix x-rays‘/s.5 kev Tatios 

of the two sources. This implies that the iso- 
meric transition is mainly converted in the 

L shell and is therefore a high electric multipole. 
Next we examined the conversion electrons from 
this transition by using a thin Ir'*”"2 source 
placed on a 1-cm-thick anthracene scintillator, 
covered with a 0.2-mg/cm? Al foil. We found 
unresolved L and M conversion electron lines 
from a 161+5 kev transition. (See Fig. 2, curve 
A.) The electron spectrum was compared with 
that of Ir**8 (see Fig. 2, curve B), which was 
quite flat in the same energy region. The energy 
calibration was carried out by means of the con- 
version lines of the isomeric transitions in Te!*>” 
(See Fig. 2, curve C.) 
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FIG. 1. Comparison of the scintillation spectra of 
the gamma rays from Ir'**8(e) and Ir'#?2 (>5 yr) (). 
The slightly higher gamma-ray intensity in the low- 
energy region is due to increased backscattering in 


the thicker Ir””2 source. The peak between the 
K x-ray peak and the 201-205 kev photopeak is due to 


backscattering in the lead shield. 
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No unconverted 161-kev gamma rays were found. 


However, because of the weak Ir'””"2 — at 
our disposal we can only say that J. poy*lsi6 kev 
<0.1. From the energy, the low K/L WL ratio, and 
the lower limit for the half-life of the isomeric 
transition, the most probable multipole assign- 
ment is E5, as follows from the values for the 
computed half-lives and activation cross sections 





10 000 r ~~ : 
7 1i6it 5 kev 
ee 


“LM CONVERSION ° 
4} LINES 









Ir 922 (>Syr) 
1 000 f 








4000 T T T T T 





Ir'92® (74 day) 






| OOO F 


COUNTS 


a 











10000 T T Tr T T 


r 30,9 kev 105.4 kev (L LINE) 
r ! 





5 

. G2 kev (K LINE) 
4 
4 
4 


1 OOO fT 











40 80 120 60 200 


o————® £( hev) 


FIG. 2. Conversion electron spectra of Ir'®?””2 
(>5 yr) (curve A) and Ir'**& (74 day) (curve B). The 
energy calibration was carried out by means of the 
K and L conversion electron lines of the 110-kev transi- 
tion and of the L conversion electron line of the 35-kev 
transition in Te'™ (curve C). The positions of the 
L and Mconversion lines of the isomeric transition in 
Ir? were determined by taking the theoretical con- 
version coefficients (see Table I) into account. 
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given in Table I. For an £5 transition of 161 key 
for Z=77, ax =6.3, ay =568 (Sliv),® a,,=238 
(Rose).* The theoretical half-life for a single 
proton £5 transition of this energy is 19.6 yr. 
Assuming that this is the actual half-life, we 
compute from the observed source strength an 
activation cross section o=5.2 x10?" cm? for 
pile neutrons. Since natural Ir was used as the 
target, we estimate that a fraction of this cross 
section, <3x10~?" cm’, is due to an (n, 2n) 
process in Ir'®. 

As Fig. 2, curve A, shows, no conversion 
electrons from a 58-kev transition were ob- 
served in the electron spectrum of Ir*®””2, 
fore, the E5 transition does not end in the 1.45- 
min state Ir'””1. Unless it is followed by a low- 
energy (<30 kev) transition, which would have 


Table I. Relation between electric multipole order, 
half-life, and activation cross section computed for a 
161-kev transition. 








E4 Es E6* 

aK 1.5 6.3 11.6 
a7? 122.2 568 5730 

ay ord 21.4 238 3760 
Bot 145.1 912.3 9501.6 
T,, (sec) © —1.55x*10° «5.65104 8.8210" 
T (yr) 3.38 x 107° 19.6 2.79 x 10° 
Gact (em) f 5.2107 5.7 x10 





“The values for the E6 conversion coefficients were 
obtained by extrapolating loga versus electric multi- 
pole order. 

see reference 3. 

“See reference 4. 

din order to take screening effects roughly into 
account, the M-shell conversion coefficients for 
Zeff=70 were used. 

©The gamma half-lives were computed using the 
formulas given by S. A. Moszkowski, in Beta- and 
Gamma-Ray Spectroscopy, edited by K. Siegbahn 
(Interscience Publishers, Inc., New York, 1955), 

p. 391. 

f The activation cross sections 7, were computed 
using the observed number of disintegrations/sec, 
the total neutron flux, and the computed half-life 
given in the line above. The following relation be- 
tween half-life T (in years) and activation cross sec- 
tion (in cm?) holds: 





o y(t M2) = 9.06% 6 -28,7,5. .06/T 


There- 








> = snr Fe 6 hlUrMhlhlUC rele lUrhUCUCUr 


Qa fteo fF BPD & =e HH oe 


Ye oO = kt eo 


ano ot wo wf ee 








ere- 








VoLUME 3, NUMBER 1 


PHYSICAL REVIEW LETTERS 


Jury 1, 1959 








escaped detection, it may be assumed to go 
directly to the 74-day ground state which is 
known to have a character 4+ or 5-.’ From this 

it follows that if ,umg=4+, then /jyis2m, = 9F 

and [j-192m, = 1+. On the other hand, if ,isg=5-, 
then Itpi92m, =10+ or 0+, and Tiyisem, =2+ or 8+, 
respectively. Since a weak beta continuum had 
been reported for Ir**71,1 of unknown end-point 
energy, we decided to study which states in Pt’? 
are populated by beta decay of this isomer. Thus 
we hoped to be able to make a spin assignment 

for Ir'*?”"1, from which the spins of the other two 
isomers would follow. The Ir”: sources con- 
sisted of enriched Ir™ (85.9%), irradiated for 

10 sec by neutrons from our reactor. We meas- 
ured the beta-spectrum of the 1.45-min activity 
and the beta-gamma coincidences, as well as 

the gamma spectrum and the gamma-gamma 
coincidences. From our results it follows that 

the Pt'* ground state (0+) and also the two low- 
est excited states at 316 kev (2+) and 611 kev 

(2+) are populated in the beta decay of Ir*”"1. 

The branching ratios found were 7x10-°, 8x10~°, 
and 2.5x10~° per disintegration, respectively.® 
Figure 3 shows the proposed disintegration scheme 
of the isomeric triplet. It contains the unique spin 
assignments for the three isomeric states derived 
from this work. (The assignment 2- for Ir!?"1, 
although compatible with the beta branching, is 
excluded, because Ir’**& would then be 5+.) Also, 
the previously established’ spin assignments of 

the levels of Pt'®*, which are fed in the beta 

decay of Ir'*”": and Ir’*8, are given. From 


Ir 92 mQ 9F 





FIG. 3. Disintegration scheme 


of the triple isomers of Ir'®?, The — IF 
decay of Ir'®?”2 and the analysis os a2 
of the beta decay of Ir'#?1, as ~— 


well as the spin assignments for 
the three isomers, result from 
the present work. Because the 
logft values for the beta transi- 
tions from Ir'® and from Ir'9?""1 
are comparable, it is not possible 
to decide which set of the beta 
transitions is allowed and which 

is first forbidden. Hence the 
ambiguity in the parities (upper 

or lower set) still persists. Only 
that part of the decay scheme of 
Ir'®°8 which is relevant for the 
discussion is shown. 






these it follows that either all three beta branches 
of Ir'*?”, are allowed and the three high-energy 
beta branches from Ir'8 (representing ~99% 

of all beta decays) first forbidden, or vice versa, 
depending on which parities (upper or lower) are 
chosen. The logft values derived from our stud- | 
ies of the beta decay of Ir'*”: are almost the 
same as those for the beta branches of Ir!” 
reported previously. It is therefore not possible 
to make definite parity assignments for the 

three isomeric states. 

By studying the 1.45-min growth of the inten- 
sity of the 468-kev gamma ray from the 4+ state 
of Pt'* which is not populated in the decay of 
Ir'*?1 we obtained the ratio for the initial acti- 
vation cross sections o;(Ir'”"1): 9;(Ir’8) = 0.84. 
Taking the total activation cross section from the 
literature,® o;(Ir'®8) =(7.0+2.0)x10- cm’, we 
obtain o;(Ir'™8) =(3.8+1.1) x 10-*? cm? and 
o;(Ir??""1) =(3.2+0.9) x10? cm*. The latter 
value is to be compared with the value of 
(2.6+1.0) x10~* given in the literature.® 

The small cross section found for formation 
of Ir'*?”"2 is in agreement with expectations 
from the rule’ that o, decreases with increas- 
ing |!-I,|, where J, denotes the spin of the 
compound state. For Ir'™+n, I, =$+ 3. 

We wish to thank E. der Mateosian and G. 
Harbottle for their kind collaboration in var- 
ious phases of this work. 
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OXYGEN-20* 


Nelson Jarmie and Myron G. Silbert 
Los Alamos Scientific Laboratory, University of California, 
Los Alamos, New Mexico 
(Received June 10, 1959) 


This is a preliminary report of the observation 
of the nuclide O”*, heretofore unknown. This 
isotope of oxygen is of interest because it is new, 
because it is a member of the group of nuclides 
with a closed shell of eight protons, and par- 
ticularly because it is the lightest nucleus known 
with isotopic spin T=2. 

We observed O° by detecting the protons from 
the O'*(t, »)O”° reaction. The experimental ar- 
rangement and equipment have been described in 
detail previously.’ Briefly, 2.6-Mev tritons 
from an electrostatic accelerator bombard a gas 
target enriched in O'*. Reaction fragments are 
analyzed in a double-focussing magnetic spec- 
trometer and are detected in a CsI crystal scin- 
tillation spectrometer. The varying sensitivity 
of the CsI to different particles and the momentum 
analysis of the fragments determine the mass and 
energy of these fragments. As an additional check 
that the fragments in the O° reaction were pro- 
tons, their behavior was observed when a thin 
aluminum foil was placed over the CsI crystal. 
Protons from targets of normal oxygen, 25% 
O'* ? and 96% O** were observed. Background 
runs were also taken with nitrogen and methane 
in the target to check on possible contaminant 
reactions. All target gases were analyzed with a 
mass spectrometer. To eliminate the possibility 
that the proton groups could have come from 
(He*, p) reactions caused by the He** component 
of the beam, runs were taken using a HT* beam. 
It was shown that the results below were defi- 
nitely from the O**(t, p)O”° reaction. 

Two proton groups associated with O° were 
observed. These are assumed to be due to the 
ground state and first excited state of O°. The 


50 


@ for the reaction O**(¢, p)O”° is 3.12 + 0.04 Mev. 
This then gives a preliminary value for the mass 
of O”° of 20.01036 + 0.00004 amu or a mass ex- 
cess of 9.65+0.04 Mev. The beta disintegration 
energy’ for the O”°-F*®° decay is calculated to be 
3.75 Mev. The first excited state of O”° is found 
to be at an energy of 1.70+0.05 Mev above the 
ground state. No other energy level was seen up 
to about 4.05 Mev. The errors stated are stand- 
ard deviations. 

This mass for O”° indicates a more stable nu- 
cleus than has been expected. For instance, 
Talmi and Thieberger* with a shell-model cal- 
culation predicted a mass excess of 11.4 Mev. 
The measured mass of O”° corresponds to a 6.5- 
Mev excited level in F*° and to a 16.7-Mev level 
in Ne”°. These levels should be the positions of 
the first T=2 state in these nuclei. 

More precise values on O*°, as well as infor- 
mation on the energy levels of O'*, N**, and N", 
will be reported in a complete paper. 

Our thanks go to Eugene Haddad for some help- 
ful advice on this experiment. 
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Work performed under the auspices of the U. S. 
Atomic Energy Commission. 

‘Nelson Jarmie and Robert C. Allen, Phys. Rev. 
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~ the 25% O'* was kindly furnished by A. O. Nier, 


University of Minnesota. 

3See S. Katcoff and J. Hudis, J. Inorg. Nuclear 
Chem. 3, 253 (1956), for an attempt to measure the 
oO” half-life. 

‘I, Talmi and R. Thieberger, Phys. Rev. 103, 718 
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Although the neutral decay mode 6,°~7°+ 7° is 
expected to occur in addition to the well-known 
charged mode 6,°~7~+7+, there has been no 
complete experimental proof of its existence be- 
cause of the difficulty of detecting the neutral de- 
cay products with high efficiency. Ridgway, 
Berley, and Collins’ and Osher, Moyer, and 
Parker* have reported counter experiments in 
which they observed gamma rays of roughly the 
correct energy originating several centimeters 
downstream of internal targets in the Cosmotron 
and Bevatron. The latter authors also report 
that the decrease in the numbers of such gamma 
rays with distance downstream from the target is 
consistent with the assumption that the gamma 
rays originate from =*, 6°, and A° decays whose 
lifetime is roughly that measured for the ordinary 
charged decays of these particles. Eisler et al.* 
have reported some gamma rays associated with 
strange particle production in a propane bubble 
chamber; Crawford et al.* have reported similar 
results from a hydrogen bubble chamber; and 
Boldt et al.* have observed showers occurring in 
conjunction with charged A° and 6° decays ina 
multiplate cloud chamber. In none of these ex- 
periments has it been possible to observe indivi- 
dual events in which all four gamma rays arising 
from the neutral 6° events were detected. 

In order to study this process in detail, as well 
as others involving gamma rays, we have con- 
structed a liquid xenon bubble chamber 30 cm in 
diameter and 25 cm in depth. The chamber oper- 
ates at -21°C at a vapor pressure of 370 psi. The 
liquid xenon has a density of 2.18 g/cm® and a 
radiation length of 3.9 cm. We have taken about 
160000 photographs of this chamber exposed to 
a1.0-Bev 1” beam at the Berkeley Bevatron. 

Figure 1 shows the production of a A° followed 
by the decay A°~p+2. Also visible in the pic- 
ture are four electron pairs whose parent gamma 
rays appear to have originated at a common point 
about 1.2 cm from the origin of the A°. If the 
measured directions of the pairs are used to ex- 
tend the lines of flight of the gammas back toward 
their origin, it is found that these lines all pass 
through a cube less than 5 mm ona side. The 
volume of this cube is fully accounted for by 


DIRECT PROOF OF 6,° NEUTRAL DECAY* 


J. L. Brown, H. C. Bryant, R. A. Burnstein, R. W. Hartung, D. A. Glaser, J. A. Kadyk, 
D. Sinclair, G. H. Trilling, J. C. Vander Velde, and J. D. van Putten 
Harrison M. Randall Laboratory of Physics, University of Michigan, Ann Arbor, Michigan 
(Received June 9, 1959) 
















































FIG. 1. Example of the decay 0,°~21"~4y. The 6° 
is produced at point 1 in association with a A° which 
decays into a proton and 1 (as indicated) both of which 
stop in the chamber. The lines of flight of the four 
gamma rays, given by the measured direction of the 
electron pairs A, B, C, and D, intersect near 2 where 
the 6,° decay took place. 


measurement errors, and therefore it can be 
concluded that within these errors the four gam- 
mas do originate at a common point. These four 
gamma rays are interpreted most naturally as 
arising from the decays 6,°~1°+7°, 1°~2y; the 

6° having been produced in association with the 
A°. We have found nine other four-gamma events 
of this same type in analyzing about one-half of 
the pictures. 

In support of the interpretation of these events 
as neutral decay modes of 6,° particles the follow- 
ing points should be noted: (1) The background of 
such four-electron-pair events arising either 
from prongless neutron stars in which two n° are 
created, or from the chance association of unre- 
lated gamma rays is completely negligible. (2) 
Although the photographs were scanned for all 
charged V° decays, with or without associated 
gamma rays, every charged V° decay associated 
with four electron pairs was a A°~p+7~, unac- 
companied by any other visible strange-particle 
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decay. This fact strongly indicates that the four 
gamma rays are indeed to be -.':ributed to a par- 
ticle of strangeness +1. (3) On the basis of the 
known 6,° lifetime, ® less than one of our ten 
events can be attributed to a neutral decay mode 
of the @,°. Although these considerations leave 
little doubt as to the fact that the four-gamma 
events do arise from neutral 6,° decay modes, 
some question may remain as to whether these 
decay modes are indeed of the type 0,°—12°+ n° 
followed by 1°~2y or whether perhaps some of 
the gammas are produced directly in the decays. 
If one knows how to pair the gammas belonging 
to the same 7°, it is possible to compute the 0° 
mass from the observed directions of the 6° and 
the four gammas, thus making a direct check 
that the decay scheme 0,°~7°+7° has really been 
observed. Unfortunately the results of this cal- 
culation are quite sensitive even to small meas- 
urement errors. Making use of the lower limit 
for the energy of each gamma ray imposed by the 
visible ionization loss of its associated electron 
shower, however, we were able to pick assign- 
ments of gamma rays which made the kinematics 
for the decays consistent with the mode 0°—7°+1° 
followed by 1°~2y for both 7s. We therefore 
believe that the decay scheme 6,°~27° is the most 





likely interpretation of the events, although on 
the basis of our data such schemes as 6,°~7°+ 2 
or 6,°~4y cannot be ruled out. 

Work is in progress on other aspects of the ex- 
periment, including branching ratios for neutral 
decay of A° and 6°. 

We wish to express our gratitude to the staff 
of the Lawrence Radiation Laboratory for making 
this work possible, and especially to our many 
friends at the Bevatron who helped us very much 
during this experiment. 
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EXPERIMENT ON CHARGE INDEPENDENCE IN PION INTERACTIONS 


D. Harting, J. C. Kluyver, A. Kusumegi, R. Rigopoulos, A. M. Sachs,* G. Tibell, 
G. Vanderhaeghe, and G. Weber 
CERN, Geneva, Switzerland 
(Received May 28, 1959) 


Although the existing experimental results in 
pion physics are consistent with the assumption 
of charge independence of pion-nucleon forces, 
the statistical accuracy of the data leads to an 
uncertainty of + 15% in quantitative tests of the 
concept.' Direct experimental tests? offer the 
possibility of substantially reducing this uncer- 
tainty before being limited by the accuracy of 
theoretical calculations of the Coulomb and mass 


difference effects in the specific reactions studied. 


The present paper describes an experiment at 
600 -Mev proton energy on the reactions 
p+d—H* +n", 
p+d—He*+n°, 


for which (apart from Coulomb and mass correc- 
tions) the charge independence hypothesis predicts 
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the ratio of the differential cross sections to be 
2 at any c.m. angle. The c.m. energy of the 
emerging pions is 220 Mev, well above the (3,3) 
resonance. The cross sections are compared by 
detecting only the recoiling nuclei, which are 
identified by momentum and time-of-flight selec- 
tion, with pulse-height analysis as an additional 
check. 

A sketch of the experimental arrangement is 
given in Fig. 1. The 600-Mev proton beam trav- 
erses a 2-cm thick liquid deuterium target with 
0.1-mm thick Mylar windows and then a secondary 
emission chamber monitor. The solid angle for 
the scattering is defined by a 2x3 cm? tungsten 
collimator, 4 cm thick, which is placed 500 cm 
from the target at an angle of 11.3° to the incident 
proton beam. Of the H* and He* nuclei transmitted, 
the higher energy group, corresponding to 52° in 
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FIG. 1. Experimental arrangement. The counter 
dimensions are given in cm. 


the c.m. system, was measured. A few proper- 
ties of these particles are given in Table I. The 
particles passed by the collimator are analyzed 
by a magnet, the current in which is adjusted to 
bend the beam of either the H® or the He® nuclei 
over 26° into a telescope of three scintillation 
counters. The separation between counter 1, 

just after the magnet, and the pair of counters 
2and 3 is 400 cm. The signal from counter 1 is 
delayed by a cable of a length corresponding to 
the time of flight of the H® and He® nuclei over 
this 400-cm path and triple coincidences are 
recorded simultaneously with two coincidence 
circuits, one with 5.5 mysec and the other with 
4.0 mysec resolution. The selectivity of this 
momentum—time-of-flight system ensures that 
no particles with a mass number different from 
three can be counted if they come directly from 
the target, i.e., are not scattered in or after the 
collimator. The identity of the counted particles 
is confirmed by recording, from one of the count- 
ers, the amplitude of the pulses that give rise to 
triple coincidences. These pulses are split and 
fed through two fast gates in parallel, each trig- 
gered by a different coincidence circuit. The out- 
put of each gate is recorded on a separate multi- 
channel pulse-height analyzer. By taking the 


Table I. Properties of the H® and He® nuclei counted 
in the present experiment. The range and energy loss 
are for plastic scintillator, the time of flight for 400- 
cm pathlength. 








Time of 
T pe/Z flight Range Energy loss 
(Mev) (Mev) (musec) (g cm™*) (Mev cm? g~) 
# 300 1332 31.1 23.2 7.2 
He? 300 666 31.1 5.8 28.9 








data from the well-defined H® and He* peaks of 
the pulse-height spectra (shown for He® in Fig. 2), 
particles of the wrong charge and accidentals 
producing small pulse heights are eliminated. 

Measurements on H° are alternated with He* 
runs. The thickness of counter 1 is reduced 
from 6 to 1.5 mm for the He* measurements to 
equalize the multiple scattering for both kinds of 
particles. Other effects of the larger multiple 
scattering of He® are negligibly small, as the 
trajectory of the heavy particles goes through 
vacuum from the target to the analyzing magnet 
and from there on through hydrogen. 

The main difficulty of the experiment is the 
very small cross section of the reactions which, 
at the angle considered, is 14 ub/sterad for H°. 
Three H® particles are counted for every 10" 
protons in our geometry, whereas the number of 
protons and deuterons scattered into the same 
solid angle is a factor of 10* higher. Most of 
these particles are removed by the analyzing 
magnet, but Fig. 3, which gives the relative 
number of triple coincidences as a function of 
the delay in counter 1, shows that the time-of- 
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FIG. 2. Pulse-height spectrum in counter 2, gated 
by triple coincidences, for He® settings. The number 
of counts is given per channel for 8 x 10"° incident pro- 
tons. 
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FIG. 3. Time-of-flight curve for H® settings. 


flight selection still has to discriminate against 
rather a large number of deuterons. 

The direct proton beam crosses the experi- 
mental room in a vacuum pipe but nevertheless 
contributes most to the background in the count- 
ers, approximately one thousand singles for 
each heavy particle recorded, notwithstanding 
the heavy shielding. To avoid corrections for 
accidental coincidences, the measurements have 
been made with a proton beam of 2x10” protons 
per second, one-fifth of the maximum intensity. 

The results of the counting, with their statisti- 
cal errors, are listed in Table Il. The value for 
the ratio of the measured cross sections follow- 
ing from these numbers is 2.29+0.05. 

In addition to the statistical uncertainty, possi- 
ble sources of systematic errors must be con- 
sidered. The reliability of the result depends 
mainly on two assumptions: first, that all the 
H® and He® nuclei, passing through the collimator, 
are counted and secondly, that all the recorded 
counts with the correct pulse height originate in 
H® and He® nuclei going through the three counters. 
To test the validity of these assumptions all meas- 
urements were repeated under marginal condi- 


Table Il. Results of the experiment. The number of 
counts is given for 10" incident protons. The errors 
include only counting statistics. 








Ho He? 
Full target 279.8 +2.6 136.2 +1.6 
Empty target 16.1+1.1 21.2+1.2 
Full - empty 263.7+2.8 115.0+2.0 
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tions. When the delay in counter 1 was set off 
from its correct value by an amount equal to the 
resolving time of the coincidence circuits, less 
than 1% of the counting rate was left. In the ratio 
this effect can be neglected. When the counters 
2 and 3 were shifted together by their own width 
in the horizontal direction, the counting rate 
dropped to less than 1% of the central position 
value. The corresponding figure after vertical 
movement was at most 2% for both H® and He’. 
Depending on their origin these marginal counts 
would have to be added to or subtracted from the 
measured numbers, leading to a negative or posi- 
tive correction of (2 +2.8)% in the ratio. As this 
origin is unknown (the counts might be due to 
multiply scattered H* and He® nuclei or to other 
scattered particles with the correct time of flight) 
no corrections were applied, but the uncertainty 
was expressed in an additional error to the ratio 
of 4%. Other sources of uncertainty are the drift 
in the monitor and possible deviations from 100% 
efficiency of the counting system, each contri- 
buting a 1% error. 

After transformation to the c.m. system and 
combining the errors in quadrature, the final 
result for the ratio of the cross sections at a c.m. 
angle of 52° is 








R =2.26+0.11. 


The difference between this result and the value 
2 predicted by isotopic spin considerations alone 
is somewhat larger than the 4% which, according 
to a preliminary calculation by Kohler,* should 
be added for Coulomb and mass corrections. 
Measurement of the ratio at a different angle is 
now planned. 

The authors wish to thank Professor G. Ber- 
nardini for helpful discussions, Dr. G. Fidecaro 
who kindly lent us the target, and Dr. B. Hedin 
for his work on the extraction of the proton beam. 
We are indebted to Mr. V. Beck and Mr. P. Actis 
for technical assistance and to the cyclotron staff 
for their co-operation. ic 
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SEARCH FOR RARE DECAY MODES OF THE ».* MESON* 


Juliet Lee and N. P. Samios 
Columbia University, New York, New York 
(Received June 4, 1959) 


With the use of a hydrogen bubble chamber we 
have endeavored to detect the following two decay 
modes of the ».* meson: 


ut—et+e~+et, 


pt—et+v+D+et+e-. 


The decay p+ ~e*++e~ +e* has not been observed 
although there are no known selection rules which 
would forbid it. Its nonexistence, however, has 
been well enough accepted to be incorporated into 
theoretical proposals such as the universal V-A 
theory’»? which excludes neutral currents. Re- 
cently Kinoshita and Sirlin® and also Eckstein and 
Pratt* have calculated the photon energy spectrum 
and decay rate for the process p+~et+v+V+y. 


As a result, the decay mode p+~et+t+v+D+et+e-, 


where the photon is internally converted, should 
occur a small fraction of the time. This last 
process would appear similar to the direct decay 
of the p»+~e*++e~ +e* in that three prongs would 
emerge from the end of a »* stopping.® A sepa- 
ration between these two processes is possible 
since the energy of the three electrons for the 
non-neutrino decay would be equal to the » mass 
(i.e., 106 Mev) and their vector momentum would 
be zero. On the other hand, in the neutrino de- 
cay the energy of the secondaries would be less 
than 106 Mev and momentum would not balance. 
It is of interest to investigate both of these re- 
actions in order to determine the validity of the 
above assumptions and calculations. 

Slow 7* mesons were stopped in a liquid hydro- 
gen bubble chamber 12 in. in diameter, 6 in. in 
depth, placed in a magnetic field of 8800 gauss. 
On the average, 10 ns stopped per picture, de- 
caying into e+ +v+D. The 7* meson served to 
identify the 1 * meson and 2.2 x10° such p*’s 
were observed. The experiment consisted in 
examining the ».+ endings for those events which 
had 3 prongs and then measuring the vector 
momenta of these secondaries. 

The possible sources of background were the 
following: 

(1) u+~e++v+D+y, where the ray was ex- 
ternally converted within 2 mm of the vertex. 
The probability for this process can be calcu- 
lated using the 7-ray energy distribution given 
by Kinoshita and Sirlin* and the known conversion 


probability in liquid hydrogen. Using a length of: 
2 mm, the number of these events expected is 
0.03. 

(2) A y ray produced by an interaction in the 
walls accidentally materializing on a p* ending. 

In all the pictures, 170 such random pairs were 
observed. A y* ending can be located to within 

a cube whose side is 2mm. Therefore, the prob- 
ability that one of the 170 pairs could have occurred 
at the end of one of the ».*+ endings is =170 x10 

x (0.2)°/10* =1.4x10~. 

Both of these effects are negligible. The scan- 
ning efficiency was determined by looking for 
events in which the positron from a yp* decay 
made a visible knock-on electron within 2 mm of 
the vertex. The comparison of two independent 
scannings of the film gave 90 % efficiency for 
finding these two-prong events. The detection of 
the three-prong events should be much easier, 
corresponding to an efficiency of the order of 95%. 

From a total sample of 2.2x10° y*’s, three 
events were found where three charged prongs 
(2 positrons and 1 electron) emerged from the p* 
ending. The total energy of the three secondaries 
was in each case 51+3.5 Mev, 5541.5 Mev, and 
33 +2 Mev and their resultant momenta 4743.5 
Mev/c, 5441.5 Mev/c, 2942 Mev/c, respectively. 
Figure 1 is a picture of one of the above events. 


ye aes 
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FIG. 1. Example of an event where 3 secondary 
prongs emerge from the p-meson ending. It is identi- 
fied as 


nt—pt+7 


et+et+e-+v4+Dd, 
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Of the three, none can be attributed to the direct 
decay p*~e*++e~ +e* since the total energy of 
the secondaries was less than 106 Mev and there 
was a net momentum unbalance. If the events 
from this decay mode were distributed according 
to a Poisson distribution, this would correspond 
to a branching ratio 


R=(ut—et+e~ +et)/(ut—et+v+d)=1x10% 


for 90% confidence. 

According to the previous criteria, the three- 
prong events are all examples of the decay 
ut—et+v+i+et+e. This gives a ratio 


R’ =(ut—et+v+d+et+e )/(pt—et+v+D) 
=(1.5+1.0)x107S. 


This ratio R’ can be estimated using the internal 
conversion probability, py, as given by Kroll 
and Wada,° 


2a 2E 11 5 E, " 
onaz na) aaa le -E )]: 
a 


where E_ =photon energy, m=electron mass, 


™m,, =meson mass [terms of order (E, /m,) 


mM 


have been neglected], and using also the transition 


probability for the radiative decay as given by 
references 3 and 4, Using a lower energy cut- 
off of 5 Mev for the y ray, this gives for R’a 





value of approximately 5x10~° which is essen- 
tially in agreement with the experimental result. 

An accurate determination of the branching 
ratio is not feasible by means of a bubble-chamber 
experiment of this type, the purpose of the pre- 
sent endeavor being to show the existence of the 
decay mode and an order of magnitude of the 
decay rate. 

We wish to thank our scanners, Mrs. Dina 
Goursky, Mr. John Impeduglia, Miss Mary John- 
son, and Mr. Alex Rytow for their painstaking 
work and the Columbia Bubble Chamber Group 
for taking the pictures. 
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PION PRODUCTION BY PIONS* 


Walton A. Perkins, John C. Caris, Robert W. Kenney, Edward A. Knapp, 
and Victor Perez-Mendez 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received June 4, 1959) 


The existence of a pion-pion interaction has 
been postulated to explain the nucleon structure’ 
and the peak’ in the 7~-proton total cross sec- 
tion near 1 Bev. The direct interaction of the in- 
cident pion and a virtual pion in the meson cloud 
surrounding the nucleon could contribute signif - 
icantly to the production of an extra pion in pion- 
nucleon collisions. In order to study this inter- 
action, we have investigated the reaction 7 +p 
—n*+ +7 +n in the energy region from 260 to 430 
Mev. The z* meson in the final state provides a 
unique signature for this reaction, since the 
highest energies are just barely above the thres- 
hold for the production of two secondary pions. 

The experimental arrangement is shown in Fig. 
1. Negative pion beams produced by the 730-Mev 
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proton beam of the Berkeley synchrocyclotron 
striking an internal target were magnetically se- 
lected and collimated as shown in the figure. The 
resulting pion beams, with kinetic energies of 
260, 317, 371, and 427 Mev and energy spread 
of +2.5%, impinged on a 4-inch-thick liquid hy- 
drogen target. 

The positive pion was unambiguously identified 
by its characteristic » decay, by use of the elec- 
tronic gating equipment described elsewhere.* 

In this experiment the coincidence 12345 triggered 
the delayed gate which made a coincidence with 
the pv pulse; the energy band of the detected pions 
was determined by the thickness of the copper 
absorber placed between counters No. 3 and No. 
4. The detection efficiency of the counter tele- 
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. Experimental arrangement (not to scale). 


scope was previously determined by using a mag- 
netically analyzed beam of positive pions. 

Measurements were taken for 1* mesons emit- 
ted at 60°, 90°, 125°, and 160° in the barycentric 
system for incident 7~ kinetic energies of 317, 
371, and 427 Mev; for a kinetic energy of 260 
Mev only a single measurement was made at 90° 
(lab). To obtain the angular differential cross 
section an integration was performed over the 
energy of the 7* meson at each angle. The angu- 
lar differential cross sections are nearly iso- 
tropic at 317 and 371 but are peaked forward at 
427 Mev. The differential cross sections were 
integrated over the angle of the 7+ meson to give 
the total cross sections shown in Fig. 2. The 
errors shown are the combined errors of count- 
ing statistics and of integration over energy and 
angle. 

Recent measurements have been made by Zinov 
and Korenchenko® at incident energies of 307, 
333, and $370 Mev for the combined reactions 
1 +p-n*++°+n and 7°+p—1 +7°+p. These meas- 
urements, taken for charged pions emitted at 
106° (c.m.), are in reasonable agreement with 
the results presented here. 

The dashed line in Fig. 2 represents the theo- 
retical prediction according to Rodberg” and 
Kazes* based on the static model which assumes 
only a direct pion-nucleon interaction. We inter- 
pret the large systematic discrepancy between 
this theory and experiment as evidence for a 
pion-pion interaction (i.e., here the 7~ interacts 
directly with the 7+ in the meson cloud surround- 
ing the nucleon). 

The following Letter® discusses the effect of a 
pion-pion interaction on the total cross section 
for this reaction. The solid line in Fig. 2 shows 
the result of including this interaction in the 
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FIG. 2. Total cross section for 7-~+p—-1++n-+n as 
a function of the incident 7~ kinetic energy. The dashed 
curve represents the static-model prediction by Rod- 


berg’ and Kazes®; the solid curve is the pion-pion 
interaction result by Rodberg. ® 





theoretical prediction. 

We should like to acknowledge the continued 
interest and support of Professor A. C. Helm- 
holz. We also wish to thank Mr. James Vale and 
the cyclotron crew for their assistance and 
cooperation during the course of the experiment. 
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PION PRODUCTION AND THE PION-PION INTERACTION* 


——_. 


Leonard S. Rodberg 


Department of Physics, University of California, Berkeley, California 


(Received June 4, 1959) 


In the preceding Letter’ the cross sections for 
1” +p—-1~ +7*+n at incident energies between 
260 and 430 Mev were presented. In this note we 
wish to discuss the theoretical interpretation of 
these results. 

One may first attempt to understand this process 
using the well-known pion-nucleon interaction. 
This has been done previously’ within the frame- 
work of the static model, which should remain 
applicable for the low energies involved here. In 
this model the incident pion interacts directly 
with the nucleon, and the nucleon “shakes off” the 
final pions. This approach yields a cross section 
for 1" +p-—1~ +2*+n which is less than 0.1 mb 
for Ejn- £400 Mev. This small cross section re- 
sults from (a) the relatively small, nonresonant, 
cross section for scattering of the initial and 
final pions on the nucleon (the incident energy is 
more than 100 Mev above the position of the 3-3 
resonance, and the mean kinetic energy of the 
final-state pions is more than 50 Mev below the 
resonance energy); (b) the weak pseudovector 
coupling (f? ~0.08) which produces the extra 
meson; and (c) the small amount of phase space 
available to the outgoing P-wave pions. Thus 
the direct pion-nucleon interaction makes a neg- 
ligible contribution to pion production in the 250- 
450 Mev region. 

We have therefore attempted to fit the results 
of the preceding Letter using a direct pion-pion 
interaction. The contribution from the inter- 
action of the incident pion with a pion in the 
nucleon cloud has been computed using the im- 
pulse approximation. The sizable cross section 
at low energies suggests an S-wave n-7 inter- 
action, and the rapid rise of the cross section 
with energy requires a P-wave interaction (all 
angular momenta are in the 1-7 center-of-mass 
system). The S-wave interaction has been des- 
cribed by phase shifts satisfying [v/(v+1)}” 
xcot 6,/=1/a, with the isotopic spin /=0 and 2, 

y the square of the 7-7 c.m. momentum, and 

a, =5A, @,=24.° The P-wave cross section has 
been computed from the low-energy form 
[A/(v+1)}“cot 6,2 =1/a,. A fit to the observed 
total cross section was obtained with |A/ =1/10 
and |a,| =1/5 (units 1, =1); the result is shown 
in Fig. 2 of the preceding Letter. The P-wave 
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phase shift obtained here is consistent with that | 


found by Frazer and Fulco* in their investigation 


of nucleon structure. | 


The calculation was performed assuming that 
only the one-meson exchange term contributes. 
The meson-nucleon rescattering corrections can 
be computed by incorporating a 1-a interaction 


into the formalism of reference 2. In the one- } 


meson approximation which was used in reference 
2, this adds an additional inhomogeneous term 
to the integral equations obtained there. These 
new equations can be solved using the same 
methods as in the earlier work. If the produc- 
tion amplitude is represented by D,, the correc- 
tion due to the rescattering of meson k is 
~[ie#5ss(*) 5ind,,(k) D,, assuming that only pion- | 
nucleon scattering in the T =J = - state is im- 


portant. There should then not be large correc- 


tions in the present energy region. As in the 


case of photomeson production,® this is because ~* 


the pion production takes place far from the 


nucleon. At higher energies, since the momen- 


tum transfer is of the order of the incident c.m. 
momentum, the mean radius of interaction be- 
comes smaller and the rescattering corrections 
become correspondingly larger. 


Over the energy range considered here the ; 


ratio of S-wave to P-wave contributions varies 
from 15 to 1/3; they become equal at Ej, = 400 
Mev. The present experiment does indicate a 
sharp break within this region. For a given 











partial wave in the 7-7 system, the present model 
predicts a pion angular distribution in the over-_ , 
all center-of-mass system favoring forward 
angles, and a pion energy distribution tending 
slightly toward higher energies than predicted by 
phase space. Both the angular distribution and 

the energy spectrum’ of the 7* are consistent 

with this at 317 and 427 Mev, but differ markedly 
at 371 Mev. This abrupt c e can be inter- 
preted as a manifestation of the rapid change in 
the S/P ratio and the resulting the S-P interfer- 
ence. It does not seem possible to make quanti- 
tative estimates of these effects since the inter- 
mediate pion is far off the mass shell. Its “mass” , 
varies with the angle and energy of the emitted 
pions, and the dependence of the cross section on ) 
this mass is difficult to estimate. The “free” 
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n-7 scattering cross section can be determined 
only by the extrapolation procedures recently 
suggested by Goebel and Chew and Low.’ 

The angular- and energy-correlations in the 
n-7 c.m. system should also be sensitive to the 
s/P ratio: a fore-aft asymmetry should appear 
within this energy range, and higher c.m. ener- 
gies should be favored as the P-wave interaction 
becomes dominant. An observation of both final 
pions is required to determine the relevant angle 
and energy. As discussed in the previous para- 
graph, it is difficult to make a reliable estimate 
of these effects using the present phenomenologi- 
cal approach. 

According to the model proposed here, the P- 
wave, T=1, a-7z interaction should dominate for 
Eine > 400 Mev. It follows that the charged-to- 
neutral production ratio (17 +p—-1~ +1* +n)/ 

(1 +p-~1~ +n° +p) is predicted to be 2, the square 
of the ratio of the pion-nucleon coupling constants 
for charged and neutral meson emission. This 
interaction also yields no production of 27°. For 
1*+p, only the reaction 1*+p—1++7° +p is per- 
mitted by the 7 =1 interaction. Asa result, the 
inelastic cross section for 1*+p should be 1/3 
that for 7 +p; experimentally,®»® for Ejnc = 500 
Mev, this ratio is about 1/4. Since the probability 
of finding two pions in the nucleon cloud is small, 
this model also predicts that double-pion produc - 





tion should be unlikely compared to single-pion 
production. 
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CONSEQUENCES OF ATOMIC CONVERSION FOR THE INTERPRETATION OF EXPERIMENTS 
ON THE SPIN-DEPENDENCE OF MUON ABSORPTION* 


V. L. Telegdi 
The Enrico Fermi Institute for Nuclear Studies, The University of Chicago, Chicago, Illinois 
(Received June 1, 1959) 


The ground state of a u-mesonic atom of non- 
zero nuclear spin J is split by the hyperfine in- 
teraction into two states, of total angular mo- 
menta F =J+4. It was pointed out by Bernstein, 
Lee, Yang, and Primakoff' (hereinafter referred 
to as BLYP) that (a) the lifetimes of the members 
of this doublet may, as a consequence of the in- 
coherence of the two states, be distinct, and that 
(b) these lifetimes will actually differ if the muon 
absorption rate depends on F, i.e., if the rele- 
vant weak interaction is at least partly spin-de- 
pendent. The total muon disappearance rate is 
the sum of the absorption and the decay rates, 
and the latter is indeed to an excellent approxi- 
mation independent of F. BLYP estimated the 


fractional difference 5 between the two lifetimes 
on the basis of a simple nuclear model,” and 
suggested an experimental test for the spin-de- 
pendence of muon absorption: Demonstrate that 
the rate of appearance of decay electrons (orig- 
inating from y~ bound in a monoisotopic target 
with J#0) does not, as a function of time, follow 
a simple exponential. The logarithm of the de- 
cay curve should indeed exhibit (in the sense de- 
fined below) a positive curvature proportional to 
&*. This quadratic dependence would of course 
prevent one, as BLYP emphasized, from estab- 
lishing which member of the doublet absorbs 
faster, i.e., from settling a question of particular 
theoretical significance. 
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In the discussion of BLYP it was assumed im- 
plicitly that atomic processes do not contribute 
to the lifetime of either state, and in particular 
that the energetically higher member of the 
doublet cannot be converted into the lower one at 
an appreciable rate. Such a conversion could 
take place in principle by the ejection of Auger 
electrons and, to a lesser extent, by M1 radia- 
tion. It is the purpose of this note to explore the 
experimental consequences of such atomic con- 
version processes which, in our opinion, may 
not be negligible in situations of actual physical 
interest. 

We refer to the higher state of the doublet by 
1, to the lower one by 2, and indicate the muon 
disappearance rates by \, and ),, respectively; 
let R stand for the conversion rate from 1 to 2. 
The populations m, and m, of 1 and 2 are governed 
by 


dn, /dt = -(,,+R)m, 
dn,/dt = - 4%, +Rn,, (1) 


as a function of time, ¢. The rate of appearance 
of decay electrons is for all ¢ proportional to 
Ao(m, +), indicating with », the muon decay rate 
(by hypothesis the same for 1 and 2). In the ex- 
periment proposed by BLYP, one measures the 
curvature K of f(#) =log(n,+m,). We adopt the 
standard definition 


K =f"(1 of*)@, 


where the primes indicate differentiation with 
respect to ¢. The sign of K is determined by the 
sign of f’’. From (1) one has directly 


Ff? = [Oy -rg)?mymy + (Ay -Ag) Rm, (n,+0,))/(m, +m), (3) 


(2) 


which shows that K> 0 when R=0, i.e., in the 
case considered by BLYP. The curvature will, 
however, necessarily be negative (K <0) when the 
inequalities 


Ay AQ, <9, (4a) 
R> ld, -A_| /(L+n,/n,), (4b) 


obtain under the conditions of the experiment. 
(4b) is of course a function of time; as K goes 
monotonically to zero with increasing ?, 


R> ld, - Ag! /[1 +, (0)/n,(0)], (4b’) 


is sufficient to ensure, with (4a), K <0 for all ¢. 
Odd Z - odd A nuclei, the most suitable targets 

for the BLYP experiment, have in general mag- 

netic moments py in agreement with the picture 
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of anJ; ,3 unpaired proton. The sign of yy de- 
termines which of the states F =/ +4 lies lowest, 
so that (4a) can be written as 


A_>Ay, When pwy>O, V=L+3) (5a) 


Ay>A. when ywy<0, Y=L-3) (5b) 
indicating with \ +(-) the disappearance rate in 
the F =1+(-)4 state.[Note that N’° is the sole re- 
levant nuclide for which (5b) applies.] The ex- 
treme case of a nucleus with .)>0 is the proton; 
condition (5a) requires hence that muon absorp- 
tion proceed faster from the singlet than from the 
triplet state in the basic absorption act. This 
will be the case when the ratio of the Gamow- 
Teller and of the Fermi coupling constants is 
negative,® e.g., for the V-A interaction now well 
established in many weak reactions. Condition 
(5b) may also be fulfilled by the same type of 
coupling, because here the muon and the proton 
are in a relative singlet state when F =] +3. 

Assuming that the members of the doublet are 
statistically populated at t=0, condition (4b’) 
becomes 


R/\d, -r_|>d/(@I+1) when Hy? (6a) 


R/Id, -r_|>@+1)/(27 +1) when <0. 


For Al (Z =13, J = §) a detailed calculation,” based 
on a V -xA type interaction and the measured‘ 
muon absorption rate, gives \_ -A, =1.2x10° 
sec™’. According to a recent estimate,* the con- 
version by Auger effect on the conduction elec- 
trons in metallic Al should occur at a rate R 

= 1(21+1)"*x1.5x10° sec’. Thus, even if this 






(60) | 





estimate were too high by a factor two, (6a) 
should be amply fulfilled for Al. As |, “ae | 
~1/Z, while R~Z® (with a>1), elements with 
Z>13 should a fortiori yield K <0. 

We conclude that in the BLYP experiment 





(a) a negative curvature can arise only if the 
basic muon absorption process is faster for anti- 
parallel muon and proton spins, e.g. for a V-A 
type interaction; 

(b) a zero curvature does not necessarily imply 
the absence of spin dependence; 

(c) the conversion rates from the higher to the | 
lower member of the p-mesonic hyperfine doublet | 
are for Z 213 sufficiently large to lead to a nega- | 
tive curvature once the condition (a) holds. 

These conclusions, in particular (c), have an 
obvious bearing on the neutron asymmetry ef- 
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fects in p-capture recently discussed by Uberall® 
The author was prompted to the considerations 

presented here some time ago by preliminary ex- 

perimental results” which suggest that K<0 for 

Al; he is grateful to R. Winston for assistance 

with the analytic formulation of the curvature 

problem. He takes pleasure in thanking H. Prima- 

koff and J. Bernstein, who had independently 

drawn similar conclusions, for much friendly 

advice. 
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Research supported by a joint program of the Office 
of Naval Research and the U. S. Atomic Energy Com- 
mission. 

‘Bernstein, Lee, Yang, and Primakoff, Phys. Rev. 
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tions, see H. Primakoff, Revs. Modern Phys. (to be 
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this at the 1958 Midwest Theoretical Physics Conference 
at St. Louis (p. 189 of the mimeographed ween 

‘J. C. Sens, Phys. Rev. 113, 679 (1959). 

This estimate was made by H. Primakoff and the 
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reference 2. R depends primarily on the probabilities 
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y =2.2 was inferred from Knight shift and optical hfs 
data. One finds, for an F=]+4— F=]-} transition, 


R=(2I+1)-\(64n/9)Z "ay (me /m,)*m,cM, 


and 10<Z*<12 for Al. 
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SUPPRESSION OF P-STATE CAPTURE IN (K,p) ATOMS* 


T. B. Day, G. A. Snow, and J. Sucher 
Department of Physics, University of Maryland, College Park, Maryland 
(Received June 3, 1959) 


It is of obvious importance, in studying the re- 
actions of K” mesons at rest in liquid hydrogen 
or deuterium, to have some idea of the atomic 
orbit from which the K mesons are captured. 
Since the low-energy scattering data’ seem to 
indicate a large S-wave interaction, it is sim- 
plest to ignore P-wave contributions and assume 
that the reactions take place from the 1S orbit. 
(In the case of deuterium, even the neglect of 
P-wave interactions does not prevent an appre- 
ciable capture from the 2P orbit.?) 

However, this neglect, in the case of hydrogen, 
of the 2P orbit contributions has been very hard 
to justify. Since transition rates due to K -p 
interactions are so very much greater than those 
of radiative transitions, even a very small P- 
wave interaction can compete favorably with the 
radiative 2P~1S transition rate, as long as the 
2P state is occupied at all. Moreover, it is 
usually argued that the cascading mesons will 
predominantly pass through the 2P state (except 
for loss due to the direct capture from higher 
nS orbits) since the nature of the selection rules 
for dipole radiation favors populating the cir- 
cular orbits. 

We would like to point out in this note that suc- 
cessive Stark effect collisions of a highly ex- 
cited (K~,p) atom with nearby protons in the 


liquid hydrogen will drastically reduce the pro- 
bability that the K” meson ever reaches the 2P 
orbit. 

Consider the following picture. The (K ,p) 
atom is moving with some velocity v(~ 10° cm/sec 
from considerations of thermal motion, recoil 
from previous radiative transitions, and possible 
pickup of the K” in flight) through the liquid hy- 
drogen. The average radius of a unit sphere 
(hereinafter called a unit cell) in liquid hydro- 
gen (density 0.07 g/cc)® is~1.8 A. Hence, on the 
average about one in ten traversals of a unit cell 
will result in the mesonic atom’s passing within 
a distance of one electron Bohr radius from a 
proton, and thus being subjected to the intense 
electric field of the proton.‘ To estimate the 
transition rates due to this Stark effect, let us 
compute it in the m=6 level, between the 6P and 
6S states. 

A measure of this rate is given by 


¥g(6P~6S) = |{6S |E-F|6P)| 
=E , (aS |rcosé|6P), (1) 
where E,,(a,) is some average electric field due 


to the proton at one Bohr radius, and 7 cosé@ 
refers to the K  meson’s coordinate. One ob- 
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tains® 


(6S |r cos@ |6P) = 30.8 ay-) (2) 


where a;- is the K-meson Bohr radius, and 
— Foe 
E 4 (ao) e* /a,*. (3) 


Then, we find 
¥g(6P-6S) =2.0x10"* sec™. (4) 


This is comparable to the direct capture transi- 
tion rate from the 6S level,* given by 


7, (ns) = (4.710'")/n® sec™, (5) 
7, (6S) =2.2 x10" sec”. (6) 


The radiative 6P-1S transition rate (obtained 
from Table XV of reference 5 upon multiplica- 
tion by the K~-p reduced mass = 648m,) is 


Y¥p=1-26x10% sec™*. (7) 


Stark transition rates for other possible (n, 1 
~n, 1-1) transitions will be comparable in mag- 
nitude to that of Eq. (4), and decrease as n de- 
creases, so that at n=2, 


¥g(2P-25) 0.2 x10"* sec™. (8) 


Note that y, is much larger than the reciprocal 
of the transit time (~10** sec™'), which, in turn, 
is much larger than YR Hence, during the time 
that a highly excited mesonic atom is within a 
Bohr radius of a hydrogen atom,’ radiative tran- 
sitions can be ignored, while many Stark transi- 
tions can take place, so that some steady-state 
distribution is approached. Ignoring, for the 
moment, the capture from the nS state, we ap- 
proximate the equilibrium distribution due to the 
Stark transitions with a simple (2/+1) distribu- 
tion.* Now, when the capture from the nS state 
is included, the mP (and, of course, the nS) states 
are immediately depopulated.* Hence, we expect 
that in each close collision, an appreciable frac- 
tion of the K” mesons in the n level are lost due 
to capture from the 2S state. 

We are now in a position to make a quantitative 
estimate of the total effect. Since the Bethe and 
Salpeter tables only go up to m=6, we use that 
level as an example of a starting level and pro- 
ceed as follows. Take 3600 K mesons distributed 
in the 6S to 6H level with a (2/+1) weighting at 
time equal to zero. Let the mesonic atoms move 
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and decay via radiative transitions for the time 
necessary to cross 10 unit cells.’® After that 
time consider the atoms as passing within a Bohr 
radius of a neighboring proton. This has the ef- 
fect, as noted above, of completely depopulating 
the 6P level and reshuffling the remaining par- 
ticles into a (2/+1) distribution again. Then 
start again with the new (2/+1) distribution on 
another 10-cell journey, etc. 

A lower limit for the total loss in the first 10- 
cell journey is™ 300 particles from 6P plus 16.9 
from 6D, 6F, 6G, and 6H combined. When the 
3600 K” mesons are gone, the loss due to 6P-€§ 
~capture will be (3/3.169) x3600 (=3408). The 
remaining K mesons are lost by direct radiative 
transitions and are distributed as in Table L 
Thus, Table I shows that only ~45 K” mesons of 
the original 3600 get to the 2P level by direct 
6D-2P radiative transitions. On following the 
cascading mesons through the lower levels, 
which are also subject to the Stark effect, one 
finds that the population of the 2P level is raised 
to ~52 out of the original 3600. 

So <1.4% of the K mesons starting in the 
n=6 level distributed according to a (21+1) dis- 
tribution reach the 2P level. Furthermore, 
similar considerations apply at even higher n 
levels. To estimate the percentage of those 
starting in a higher level which reach the 2P 
level, we use the fact obtained from the n=6 
level calculation, that most of the particles 
reaching 2P came via direct 6D~2P radiative 
transitions. This percentage was directly pro- 
portional to the radiative transition rate for 
6D-2P. For a given initial 1, the radiative tran- 
sition rate goes as n~* (see reference 4). 

For example, at n=15, the S-state nuclear 
capture rate is still much faster than all radia- 
tive rates and the inverse of the transit time. 
Hence, our conditions are still fulfilled, and 


Table I. Distribution by radiation from n=6 levels, 
starting with 3600 particles distributed in the 6S to 
6H states with a (2/+ 1) weighting. (All levels shown 
are taken as stable for purposes of the diagram.) 





5P/3.71 5D/8.9 5F/17.5 5G/31 
4P/7.98 4D/16.0 4F/21.5 

3P/17.36 3D/26.1 

2P/44. 59 
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now, of those mesons which start in an n=15 
level, the percentage that reaches the 2P state 
pecomes ~0.1%. Thus, our calculations indicate 
that of the K mesons captured by protons in 
liquid hydrogen in highly excited states, <1% 

will ever reach the 2P level.” 

The above result implies that when a K meson 
comes to rest in liquid hydrogen, the nuclear 
capture reaction essentially always occurs with 
the relative orbital angular momentum of the 
(K",p) system equal to zero. This has many im- 
portant consequences: 

1. The isotropic angular distributions of the 
decays of £*, = , and A° hyperons'® produced 
from (K ,p) captures at rest shows, according 
to the arguments of Treiman, ‘* that the spins of 
these three hyperons are each equal to one-half. 

2. The relation between the data for (K , p) 
reactions in flight and those for (K~,p) captures 
at rest, assumed by many authors,** is now 
justified. In addition, the (A°,2z) production 
data from (K ,p) atoms can be used to determine 
the parity of the K ~ meson relative to the (A°, p) 
system.’® 

3. The predictions of Amati and Vitale*’ con- 
cerning the branching ratios from (K ,p) capture 
reactions at rest, based on the assumption of 
global symmetry, may now more significantly be 
compared with existing data. Subject to the re- 
servations pointed out by Dalitz,**® the data vio- 
late the predicted inequality. 

4. It is expected that the Stark effect transi- 
tions will also dominate in the capture at rest of 
p, =, and = hyperons in hydrogen, so that it 
is justified to assume that the capture reactions 
take place from S states. 

5. Similar considerations are expected to hold 
for captures at rest of K , =, and = particles 
in deuterium. 

The results of this calculation indicate that 
essentially no mesonic x-rays should be produced 
when aK meson comes to rest in liquid hydro- 
gen. This result can be checked experimentally 
with a search for such x-rays. 

We are indebted to Professor L. Madansky for 
Suggesting the possible importance of the Stark 
effect in (K ,p) mesonic atoms, and to Dr. R. G. 
Glasser for directing our attention to the higher 
n levels. 
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transit time. However, for m=2, IT ~ 2x10'? sec™! 
compared with yp(2P—1S) = 4x 10'! sec™!. Thus, in 
integrating over the path of the (K ,p) system in a 2P 
state, the Stark effect is rather small, especially when 
the 2S level shift due to nuclear capture is included 
{not done in Eq. (b) above }. 

‘°This average time for crossing one unit sphere is 
#(2Ro/v) where A, is the radius of the unit sphere 
~1.8A, and v is the velocity. The same factor $ ap- 
plies in crossing, on the average, a sphere of radius 
equal to one Bohr radius. 

'!'This estimate of the loss underestimates the true 
effect since successive repopulation and depletion of 
the nP level many times within one collision is ignored. 

12 Clearly the results depend on the exact value chosen 
for the velocity. The loss due to radiation, and hence 
the percentage reaching the 2P state, is inversely 
proportional to this velocity. However, the velocity 
would have to be changed by a large factor to change 
the absolute percentages appreciably. 

137,. Alvarez et al., Nuovo cimento 5, 1026 (1957). 

4S. B. Treiman, Phys. Rev. 101, 1216 (1956). 
SR. H. Dalitz and S. F. Tuan, Ann. Phys. (to be 


63 















VOLUME 3, NUMBER 1 


PHYSICAL REVIEW LETTERS 


JuLy 1, 1959 





published); Jackson, Ravenhall, and Wyld, Nuovo ci- 
mento 9, 834 (1958). 

‘61. B. Okun’ and I. Ia. Pomeranchuk, J. Exptl. 
Theoret. Phys. (U.S.S.R.) 34, 997 (1958) [transla- 
tion: Soviet Phys. JETP 34, 688 (1958)]. 


‘TD, Amati and B. Vitale, Nuovo cimento 9, 895 
(1958). 

18R. H. Dalitz, 1958 Annual International Conference 
on High-Energy Physics at CERN, edited by B. Ferrettj 
(CERN, Geneva, 1958), p. 197. 











S-MATRIX POLES IN THE DETERMINATION OF PARITIES* 


G. Feldman and T. Fulton 
John Hopkins University, Baltimore, Maryland 
(Received May 18, 1959) 


Recently it has been conjectured’ that there 
exist dispersion relations in momentum transfer 
variables. A number of authors”»* have suggested 
using an approximation to these relations in order 
to determine such physical quantities as parities, 
coupling constants, etc. This approximation con- 
sists of extracting the contribution of certain 
single-particle intermediate states (“the pole 
terms”) to S-matrix elements, and analyzing the 
experimental data on the assumption that these 
contributions predominate. In this note we wish 
to point out some potential difficulties inherent 
in such an analysis. 

One of the processes which has been examined 
in reference 2 is the production of A’s in 1-proton 
collisions, that is 


1 +p~A+K°. (1) 


If we assume that the Mandelstam conjecture’ is 
justified, there is a term in the matrix element 
of process (1) from the K-meson intermediate 
state illustrated in Fig. 1(a). This term will con- 
tribute provided the vertex parts 


(K°lj_4!K*), (Alj,-+|P) (2) 





FIG. 1. Possible pole contributions for the processes 
n~+p—-A+K* and 1-+p—-n+r°, 
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do not vanish. In this case the matrix element 
will have a pole in cosé@ (@ being the angle between 
the incoming 7 meson and outgoing K in the center- 
of-mass system) at a point a,. This point is in 
the nonphysical region, but can be near cosé=1. 
There are, of course, contributions to the ma- 
trix element arising from other intermediate 
states, for example the ones in which there is 
both a K and am meson. These terms lead toa 
branch point in cosé@ at the point a,. Specifically, 
for a 7-meson lab energy of 1.9 Bev, a, =1.24 
and a, =1.40. 

If one plots F(cosé@) using the experimental data, 
where 


F(cosé) =(a, - cos6)*da/dQ, (3) 


it is claimed that one can obtain direct evidence 
for or against the presence of the pole term. The 
single experiment, namely the angular distribu- 
tion for process (1) at a fixed energy, could be 
used to prove that the vertex parts (2) do not 
vanish. Hence it would follow that the relative 
K* -K° parity is odd. 

We wish to examine reactions analogous to 
process (1). Our aim is to show that use of the 
type of analysis discussed above may lead to 
ambiguous, and possibly incorrect, conclusions. 
Accordingly, we shall first consider the experi- 
mental angular distribution* for the process 


1 +p—7° +n (4) 


at a fixed 7-meson kinetic energy of 220 Mev in 
the laboratory system. In strict analogy to the 
method outlined for process (1), we consider 
whether the experimental information sheds any 
light on the existence of the 7+-meson inter- 
mediate state illustrated in Fig. 1(b). The matrix 
element M(cos@) (@ is the center-of-mass scatter- 
ing angle) for process (4) is, apart from kinema- 
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tic factors, of the form 


fgu Tu 
M(cosé@) . eh T(cosé), (5) 
where 
f=(m*1§0) gla"), (6) 
gu Tu, =(p19(0)_ _In) , (7) 


and I’ =1(y,) depending on whether the 7* has even 
(odd) parity relative to the p-m system. The con- 
stant c, is given by 


Cy =1+ 3p? /k?, (8) 


where k is the 7-meson momentum in the center- 
of-mass system. The next singularity, a branch 
point arising from intermediate states of two 

m7 mesons, occurs atc,, where 


C, =1+2p?/k?. (9) 
At 220-Mev lab energy, 
C =1.16 and c, =1.64. (10) 
We define 
G(cosé@) =(c, - cosé)*da/dQ. (11) 


G(cos@) is a function regular at cos@=c,. Since 
the nearest singularity is at cos@=c,, one may 
expand G about c, in a power series whose radius 
of convergence is c,-c,. Chew* suggests the use 
of experimental data within this circle of conver - 
gence in extrapolation to the pole.® If we further 
assume that the pole is dominant in this region, 
we obtain the result 


G(cos6) =(c, - cosé)?2> nates | M(cosé) |? 


~f*g?(1 -cosé), (12) 


where we have taken ['=y,. (The result using 
I'=1 is markedly different.) 

In Fig. 2 we plot the experimental points for 
G(cosé) in the region 1 >cos@2=0. We observe 
that these points can be fitted to a straight line 
passing through cos@=1, i.e., the distribution 
(12). One might consider this as evidence for 
the presence of the 7+-meson intermediate state. 
This implies that the vertex part (6) is nonzero 
and forces the untenable conclusion that the 7° 
meson is a scalar (or that parity is not conserved 
in strong interactions). Also plotted in Fig. 2 is 
the result obtained using the best fit for the ex- 
perimental do/d& as determined in reference 4. 





G (cos 8) 











Cos © 


FIG. 2. Plot of G(cos@) =(Cy-cos0)*da/dQ (cy =1. 16) 
for 1~-p charge-exchange scattering at 220 Mev. The 
experimental data are from reference 4. The straight 
line is obtained on the assumption that the pole term 
is dominant. The dashed curve is based on the best 
fit to do/dQ given in reference 4. Errors in this curve 
are not shown. 


If we apply the Chew criterion® to this case, we 
would be restricted to using experimental data 
in the region 1 2 cos@ 20.68. The best fit to the 
experimental data extrapolated to this region no 
longer lies on the straight line through cosé =1. 
However, one need not have assumed that only 
the pole contributes. To prove that the pole term 
is not present, it is necessary to show unambig- 
uously that G(cos@) has a second order zero at 
cosé =C,. 

One can carry out a similar procedure for the 
reaction 


1 +p—-m +p, (13) 


at 220 Mev. Again, in the absence of other in- 
formation, we can ask whether the experimental 
angular distribution indicates the existence of 
the 7° intermediate state. The analysis proceeds 
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as before. Once more, one cannot rule out the 
presence of this pole. Its existence would imply, 
in this case, that parity is not conserved in 
strong interactions. 

The illustrations we have employed indicate 
that attempts to establish the presence or ab- 
sence of certain vertex parts (and thus determine 
parities) by the above means are, at best, ques- 
tionable. Of course, to make our point we have 
been selective in the choice of experiments. 
Nevertheless, this is in keeping with the spirit 
of the pole approximation as presently advocated. 

The conjecture that poles in cosé@ appear in the 
S matrix is no doubt reasonable. However, an 
analysis of experiments which assumes that the 
poles predominate may require accurate meas- 
urements sufficiently close to the pole, yet far 
from other singularities. What is meant by 
“accurate,” “sufficiently close,” and “far” is 
somewhat nebulous. In fact, in the above illus- 
trations [and in process (1)], in order to establish 
the presence or absence of a pole, it is necessary 


to distinguish between a possible zero near cosé@ =1 


and a double zero at cos@=c, (or a,). At higher 


energies the poles approach the physical region, 
that is the points cos@=+1. Analysis of experi- 
ments at these higher energies probably would 
not eliminate the ambiguities, since the other 
singularities also approach these limit points. 





* 
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‘The corresponding circle of convergence for process 
(1) lies entirely outside of the physical region. Never- 
theless, Taylor? has suggested using data in the region 
12 cosé = 0. 
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INTEGRATION OF SECONDARY CONSTRAINTS 
IN QUANTIZED GENERAL RELATIVITY. P. W. 
Higgs [Phys. Rev. Letters 1, 373 (1958)]. 


In this Letter’ it was stated that three of the 
secondary constraints’ in the quantized theory of 
gravitation are equivalent to the statement that 
the wave functional ¥{g,.(x)} is invariant under 
general transformations of the three coordinates 
x. From this the conclusion was drawn that ¥ 
must depend only on the three eigenvalues of the 
Ricci tensor Ry, (x). (As before, Latin indices 
run from 1 to 3; Greek indices from 0 to 3.) It 


appears on more careful examination of the trans- 


formation properties of © that the former state- 
ment is not quite correct, so the conclusion is 
erroneous. The purpose of the present note is 
to clarify the invariance properties of ¥ and to 
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present the correct solution of the secondary 
constraints 


KH W=0, (1) 
where 


Ys rs 
x = -2 
u Bysu" @,7 ) s 


a 

with 7 = ~46/0g 
Let us consider just the localized transforma- 

tions of the three spatial coordinates, x’” =x” 

+a’(x), where the infinitesimal functions a’ (x) 

vanish at infinity. Then 


u u u 
og, ,(2)=-(a Bysu*? Pus” enw” 


and so, the primary constraints having been 
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satisfied already, 
oy =(d°x dg, (e)6¥/dg, (x) 


oa ifeex a (x), (x) - 2 dSa"(x)g x” 


The first term vanishes if the constraint (1) is 


satisfied; the second is an integral over the sur-' 


face at infinity which vanishes only if a“(x) is 
localized. We therefore have to construct func- 
tionals ¥{g,.(x)} which are invariant only under 
localized transformations. 

The problem is solved by transforming Ers 
into a form which satisfies certain transversality 
conditions: 


ae 8 FT 
8%) =9 ? Ff ab?” 


If y4(x) are three independent functions such that 
the differences y*-x@ are localized and a" ab 
satisfies suitable conditions, then we know that 
¥ is independent of y*. By suitable conditions 
we mean that a" ob has only three independent 
components and the functions y? are uniquely 
determined by 85%) together with these condi- 
tions. It turns out that the well-known harmonic 
coordinate conditions fulfill these criteria. 

It is convenient to introduce the contravariant 


densities h”S, related to g,. by h”§g.4=(- g)”75"). 


Then we define the “transverse gravitational 
potentials” 


M)=(dety® yy? RG), @) 


where y@(x) is the regular solution of the gener- 
alized Laplace equation 


Ys 
w"%y |) =0 (3) 


which satisfies the asymptotic condition (locali- 
zation) 


lim (y? -x") = 0. 
|x | 





The functions y? and consequently also the trans- 

verse potentials are uniquely determined by 
rs) The constraints (1) now tell us that 

v! 2 y{nTao (y)}. The transversality condition, 

which follows from (2) and (3), is 


h =0, (4) 


so there are only three independent transverse 
components. They may be exhibited explicitly 
by performing a Fourier transformation and 
resolving the transform of nTab parallel and 
perpendicular to the wave vector k,: we define 


Ye) =[ety exp(-ik-y)e! 27) +0}, (5) 


where ef (k) is an orthonormal triad of vectors 
_ scalar product here is Euclidean) such that 

=\|kle*g. Then, by (4) and (5), the only non- 
maha components of y4) are y"', , Ya. 
[In addition, the reality condition 74J(-k) =y*#(k) 
is to be imposed.] Finally, the wave functional 
may be written as ¥{y4J(k)}. 

The transverse potentials which have been 
defined here reduce in the linearized theory to 
those which were obtained by Arnowitt and Deser,* 
except that their 7“/ was traceless on account of 
the eighth constraint. It remains to be seen 
whether Dirac’s Kr condition implies a similar 
restriction in the full theory. 


‘Equation (11) of this paper should read 


u.@ 
(Hx), H po= 1H 1 (x)8/8x 5 (x-y). 
2p, A. M. es Proc. Roy. Soc. (London) A246, 
333 (1958). 
5R. Arnowitt and S. Deser, Phys. Rev. 113, 745 
(1959). 





EFFECT OF NUCLEAR ELECTRIC DIPOLE 
MOMENTS ON NUCLEAR SPIN RELAXATION 
INGASES. P. A. Franken and H. S. Boyne [Phys. 
Rev, Letters 2, 422 (1959)]. 


In this Letter it was estimated that careful 
measurements of nuclear spin relaxation times 
in noble gases such as He* and Xe!”® could re- 


veal the existence of nuclear electric dipole 
moments as small as 10~° or 10~* nuclear mag- 
neton. (One nuclear magneton =eh/2Mc =e x107™ 
cm.) Indeed, we interpreted the already available 
measurements on these gases to indicate upper 
limits of order 10~* to 10~* nuclear magneton. 
Professor E. M. Purcell has brought to our 
attention that we have made a serious error in 
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our implicit assumption that the velocity changes 
due to collisions in the gas are uncorrelated. 
Whereas an atom does experience a random walk 
in its displacements it is not true that the changes 
in velocity due to collisions are random. For 
example, if an atom is moving in the x direction 
it is more than 50 % probable that the next colli- 
sion will decrease the x component of velocity. 
Were this not true the atoms in the gas would 
rapidly depart from the Maxwell-Boltzmann 
velocity distribution. 

Since we implicitly assumed that the changes 
of velocity are uncorrelated, it becomes apparent 
that our calculation must be modified. Purcell 
has found, and we entirely agree, that our calcu- 
lation is valid only when w7, 21, where w is the 
Larmor frequency of the nuclear magnetic moment 
in the applied magnetic field and 7, is the mean 
time between collisions. 

In the event that w7-<1, Purcell finds, and we 
entirely agree, that the relaxation time we com- 
puted [Eq. (2) of the Letter] must be multiplied 
by the approximate factor [1/w7,.' Therefore 
the estimated maximum electric dipole moment 
must be multiplied by the approximate factor 
[1/w7,], which is of order 10° for the experiments 


discussed in the Letter. 

It seems entirely possible that experiments 
could be done with He and Xe under conditions 
where w7,=1, in which case the formulas of 
the Letter are applicable. For magnetic fields 
of ~ 20000 gauss and room temperature we esti- 
mate that gas pressures of ~1/20 atmosphere 
would be satisfactory. We see no objection to 
compressing or even liquifying the gases for 
short periods of time in order to perform the 
necessary spin polarization measurements. We 
estimate that upper limits of ~10~* nuclear mag- 
neton could still be established, but the experi- 
ments are difficult. Such upper limits would still 
be an improvement over what can be obtained by 
previously discussed methods.? 

We wish to express our gratitude to Professor 
Purcell for bringing this correction to our atten- 
tion, and for several valuable discussions. 





‘Professor Purcell has performed an exact calcula- 
tion of this factor which will be submitted to the Physi- 
cal Review. 

*E. M. Purcell and N. F. Ramsey, Phys. Rev. 78, 
807 (1950). 
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In this section are printed the abstracts of Articles that 
have been forwarded to The American Institute of Physics 
for publication in THE PHYSICAL REVIEW. In quoting 
information obtained from this section before the appear- 
ance of the corresponding Article, reference should be 
made to “Physical Review (to be published)” rather than 
to this Journal. 


SINGLE-PARTICLE ENERGIES IN THE MANY- 
FERMION SYSTEM. David J. Thouless, Lawrence 
Radiation Laboratory, University of California, 
Berkeley, California (Received January 19, 1959). 


We have examined the ground state of the many- 
fermion system, using perturbation theory, with 
particular attention to Brueckner’s method. A 
study of the single-particle propagators in the 
ground state of the system sheds some light on 
the relation between the single-particle energies 
used, for example, in the optical model of nuc- 
lear reactions, and Brueckner’s self-consistent 
energies. This leads naturally to the use of a 
renormalization factor to allow for propagation 
off the energy shell. It is possible to formulate 
the perturbation theory entirely in terms of graphs 
with no fermion self-energy parts, and this method 
is used to examine the approximations used for 
this problem, and to relate more closely the cal- 
culations of the ground-state energy and of other 
properties of the ground state. 


LOW-TEMPERATURE TRANSPORT PROPER- 
TIES OF COPPER AND ITS DILUTE ALLOYS: 
PURE COPPER, ANNEALED AND COLD-DRAWN. 
Robert L. Powell, Hans M. Roder, and William 

J. Hall, Cryogenic Engineering Laboratory, 
National Bureau of Standards, Boulder, Colorado 
(Received January 22, 1959). 


Experimental results between 4° and 300°K are 
given for (1) the thermal conductivity, electrical 
resistivity, and thermoelectric force and power 
of two high-purity coppers, one annealed and one 
cold-drawn 26%; and (2) the electrical resistivity 
of a series of seven samples of high-purity cop- 
per cold-drawn between 0 % and 20 % elongation. 
The total electronic thermal resistivities each 
consist of three terms: the intrinsic resistivity, 
W;; the imperfection resistivity, W,; and a devia- 
tion term, Wj,, indicating the departure from 
strict additivity of W; and W,. The intrinsic 
thermal resistivity and intrinsic electrical re- 





sistivity vary as T*-* and T**®, respectively, 
contrary to the predictions of the usual transport 
theory using Bloch approximations and assump- 
tions. The resistivity of pure copper is 1.545 
ohm cm at 0°C. The increase in imperfection 
electrical resistivity is approximately linear 
with increase in cold-drawn elongation. How- 
ever, the added resistivity is not independent 

of temperature (Matthiessen’s rule), but is about 
twice as great at the ice point as it is at 4°K. 
The change in thermoelectric power with draw- 
ing is positive at the lower temperature, but 
negative above 38°K. The Lorenz number does 
not approach the Sommerfeld value at the lowest 
temperatures, but flattens out to a value con- 
siderably smaller. A qualitative discussion for 
each of the various effects is given. 


RECOVERY OF ELECTRON-IRRADIATED 
COPPER. I. CLOSE PAIR RECOVERY. J. W. 
Corbett, R. B. Smith,* and R. M. Walker, Gen- 
eral Electric Research Laboratory, Schenectady, 
New York (Received January 28, 1959). 


Stage I recovery (14-65°K) of electron-irradi- 
ated, high-purity copper is shown to be com- 
posed of at least five substages: I,, 14-24°K; 
Ip, 24-28.5°K; Ic, 28.5-33°K; Ip, 33-48°K, and 
Ip, 48-65°K. 1,, Ip, and I, recovery are dis- 
cussed in further detail. (Discussion of I,, and 
Ir recovery is deferred to the following paper.) 
I B and I; are shown to have the characteristics 
of close pair recovery. It is inferred that I, also 
has these characteristics. The values of the 
activation energy for recovery are E 4 =0.05 
+0.01 ev, Ep =0.085+0.01 ev, and FE, =0.095 
+0.01 ev. An attempt to relate the observed 
close pair spectrum to specific interstitial- 
vacancy configurations is described. 


"Now at Hanford Laboratory Operations, Hanford 
Atomic Products Operations, Richland, Washington. 


RECOVERY OF ELECTRON-IRRADIATED COP- 
PER. II. INTERSTITIAL MIGRATION. J. W. 
Corbett, R. B. Smith,* and R. M. Walker, Gen- 
eral Electric Research Laboratory, Schenectady, 
New York (Received January 28, 1959). 


Several experiments are presented dealing 
with substages Ip and Ig of stage I (10°K - 80°K) 
recovery in electron-irradiated pure copper. It 
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is concluded that both Ip and If represent the 
recovery due to the free diffusion of a defect, 
presumably an interstitial atom, with an activa- 
tion energy for motion E,,, =0.12+0.005 ev. Ip 
is the result of correlated recovery (the inter- 
stitial returns to the vacancy from which it 
came) while I, is the result of uncorrelated re- 
covery (the interstitial travels to a distant sink). 
The migrating interstitial is shown to interact 
with vacancies, traps, and other interstitials 
leading, respectively, to annihilation, trapping, 
and cluster formation. The interaction between 
two interstitials leading to cluster formation is 
found to be comparable with the annihilation in- 
teraction between an interstitial and a vacancy. 
Clusters containing more than two interstitials 
are also formed. Previous work in copper is 
discussed, and it is concluded that this previous 
work is consistent with the above description of 
stage I recovery. The relationship of the pre- 
sent work to the general problem of assigning 
defect processes to the higher temperature re- 
covery states in copper is also discussed. 


*Now at Hanford Laboratory Operations, Hanford 
Atomic Products Operations, Richland, Washington. 


SPIN-WAVE SPECTRA FOR CANTED ANTIFER- 
ROMAGNETS AND FERROMAGNETS. R. Orbach, 
Department of Physics, University of California, 


Berkeley, California (Received January 16, 1959). 


The magnetic resonance conditions and the 
spin-wave spectra are found for canted antifer- 
romagnetic and ferromagnetic lattices, where 
the cant is produced by magnetocrystalline ani- 
sotropy fields which are noncollinear. The sub- 
lattices are thereby caused to cant towards each 
other in antiferromagnets and away from each 
other in ferromagnets. The cant of the antifer- 
romagnetic sublattices may produce a net moment 
(weak ferromagnetism) but does not alter appre- 
ciably the usual antiferromagnetic spin-wave 
spectrum in the presence of anisotropy. The 
static susceptibility parallel ( i) and at right 
angles (x,) to the vector difference of the ani- 
sotropy fields is calculated, and is shown to be 
altered from the noncanted result. It is shown 
that in antiferromagnets with no apparent weak 
ferromagnetism, canted sublattices may still be 
present, and can be detected by a nonzero ratio 
of x, to x, at 0°K. The ferromagnetic spin-wave 
spectrum shows a sudden change from the normal 
spectrum as soon as one introduces the noncol- 
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linear anisotropy fields. An optical branch is 
formed and a high-frequency k =0 magnetic 
resonance is expected. This resonance is a con- 
sequence of the two-sublattice character of the 
canted ferromagnet and may be termed an ex- 
change resonance. 


PRESSURE EFFECTS IN LUMINESCENCE. I. 
ISOBARIC EXPERIMENTS ON Nal(T1). L. Reiffel, 
Physics Research Division, Armour Research 
Foundation, Chicago, Illinois (Received January 
5, 1959). 


The effect of steady-state hydrostatic pressure 
on the properties of luminescent systems, with 
particular reference to the phosphorescence of 
Nal(Tl), is examined. It is shown that readily 
observable changes in decay mean life may be 
produced by very modest pressures. An inter- 
pretation of these effects is developed from a 
configurational coordinate model and also from 
thermodynamic considerations. Experimental 
results are in good agreement with these theo- 
retical expectations. For a 0.66-ev trap in 
Nal(T1), the change in center volume at the point 
of transition from metastable to emitting state is 
found to be +17.5x10™"* cm*®. Assuming a spe- 
cific model for the center allows computation of 
the sign and magnitude of the critical displace- 
ment of ions required for trap collapse—in this 
case +0.19 A. The force constant for the meta- 
stable state is computed to be 5.8 x10° dynes/cm; 
and the vibrational frequency, assuming an ef- 
fective mass equal to six I ions, is 3.410" 
sec * corresponding to a zero-point energy for 
this trap of 0.007 ev and crossover at a vibra- 
tional quantum number of about 47. The trans- 
mission coefficient for crossover to the emitting 
state appears to be of the order of 0.03%. Other 
possible experiments involving steady-state or 
transient pressure effects are outlined. 


TEMPERATURE DEPENDENCE OF THE WORK 
FUNCTION OF SILVER, SODIUM, AND POTAS- 
SIUM. C. R. Crowell and R.A. Armstrong,* 
Department of Physics, McGill University, 
Montreal, Quebec, Canada (Received October 
$1, 1958). 


Measurements have been made of the tempera- 
ture dependence of the work functions of freshly 
evaporated films of silver, sodium, and potas- 
sium deposited at pressures of the order of 10° 
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mm of mercury. The results for silver (-13.4 
x1075 volt/°K) were obtained by a capacitive 
contact potential method, those for sodium and 
potassium (-51x10-* and -26x10~° volt/°K, 
respectively) by the diode retarding-field method. 

The theoretical outline by Herring and Nichols 
has been extended to produce quantitative pre- 
dictions in fair agreement with the experimental 
results. Quantum defect calculations of the in- 
ternal energy have been used to determine the 
effect of thermal expansion. Some allowance for 
anisotropy and dispersion has been included in 
the calculations of the electrostatic effect of 
thermal vibrations at constant volume. This ef- 
fect is the largest contribution to the tempera- 
ture dependence of the work functions of mono- 
valent metals. The calculations point out the 
importance of the predominantly transverse 
modes of vibration. The effect of thermal vibra- 
tion at constant volume on the chemical potential 
has also been calculated allowing for anisotropy 
and dispersion, and is shown to be small (~ -2 
x10°5 volt/°K for monovalent metals). The theo- 
retical predictions (representing the assumed 
model within an estimated 25%) are -12x10°°, 
-45x10°°, and - 28x10~° volt/°K for silver, 
sodium, and potassium, respectively. 


"Now at Radio and E. E. Division, National Research 
Council, Ottawa, Canada. 


NEW APPROACH TO THE THEORY OF SUPER- 
EXCHANGE INTERACTIONS. P. W. Anderson, 
Bell Telephone Laboratories, Murray Hill, New 
Jersey (Received February 4, 1959). 


The theory of indirect exchange in poor con- 
ductors is examined from a new viewpoint in 
which the d (or f) shell electrons are placed in 
wave functions assumed to be exact solutions of 
the problem of a single d electron in the presence 
of the full diamagnetic lattice. Inclusion of d- 
electron interactions leads to three spin-depend- 
ent effects which, in the usual order of their 
sizes, we call: superexchange per se, which is 
always antiferromagnetic; direct exchange, 
always ferromagnetic; and an indirect polariza- 


tion effect analogous to nuclear indirect exchange. 


Superexchange itself is shown to be closely re- 
lated to the poor conductivity, in agreement with 
experiment. By means of crystal field theory 
the parameters determining superexchange can 
be estimated, and in favorable cases (NiO, 


LaFeO,) the exchange integrals can be evaluated 
with an accuracy of several tens of percent. 
Qualitative understanding of the whole picture of 
exchange in iron group oxides and fluorides fol- 
lows from these ideas. 


VIBRATION FREQUENCY SPECTRA OF DIS- 
ORDERED LATTICES. I. MOMENTS OF THE 
SPECTRA FOR DISORDERED LINEAR CHAINS. 
C. Domb* and A. A. Maradudin, Physics Depart- 
ment, University of Maryland, College Park, 
Maryland, and E. W. Montroll and G. H. Weiss, 
Institute for Fluid Dynamics and Applied Mathe- 
matics, University of Maryland, College Park, 
Maryland (Received January 19, 1959). 


Using the theory of random walks on lattices, a 
combinatorial expression has been obtained for 
the even moments of the vibrational frequency 
spectrum of a randomly disordered, two-com- 
ponent, isotopic linear chain as functions of the 
concentrations of the two kinds of particles and 
of their mass ratio. Expressions for the even 
moments up to y,, are presented. 


“On leave from Physics Department, King’s College, 
University of London, London, England. 


VIBRATION FREQUENCY SPECTRA OF DIS- 
ORDERED LATTICES. 0. SPECTRA OF DIS- 
ORDERED ONE-DIMENSIONAL LATTICES. C. 
Domb* and A. A. Maradudin, Physics Depart- 
ment, University of Maryland, College Park, 
Maryland, and E. W. Montroll and G. H. Weiss, 
Institute for Fluid Dynamics and Applied Mathe- 
matics, University of Maryland, College Park, 
Maryland (Received January 19, 1959). 


Using a combination of the moment-trace meth- 
od and a new method, the “delta-function” meth- 
od, the vibrational frequency spectrum of a ran- 
domly disordered, two-component, isotopic 
linear chain has been computed for a wide range 
of the concentrations of the two kinds of particles 
and of their mass ratios. In addition the parti- 
cular case of a chain in which the mass of one of 
the isotopic constituents becomes infinite can be 
treated exactly, and the results of this analysis 
shed light on the form of the spectra for lattices 
with large but finite mass ratios for the two con- 
stituents. The spectra are characterized by the 
disappearance of the square-root singularity at 
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the maximum frequency which is found in or- 
dered one-dimensional lattices, and by the ap- 
pearance of impurity bands, the nature of which 
is discussed. Finally, the zero-point energy of 
a randomly disordered lattice is calculated and 
compared with the zero-point energy of an or- 
dered lattice and of the separated phases. 


+ 
On leave from Physics Department, King’s College, 


University of London, London, England. 


ELECTROLUMINESCENCE OF AIN. G. A. 
Wolff, I. Adams, and J. W. Mellichamp, United 
States Army Signal Research and Development 
Laboratory, Fort Monmouth, New Jersey (Re- 
ceived January 12, 1959). 


Electroluminescence in the visible region is 
observed when AIN is excited by either ac or dc. 
The emitted spectra consist of narrow bands in 
the 400-500 my region and broad continuous 
bands in the 500-700 my region. The broad 
bands vary according to the activator agents. 
With ac, the luminescence is in phase with the 
applied voltage. With dc, light is emitted in the 
vicinity of the cathode and decays with time. 
This decay can be explained by polarization. The 
presence of a barrier at the cathode is indicated. 
Impact excitation is proposed. 


STUDIES OF THE SEMICONDUCTING PROPER- 

TIES OF THE COMPOUND CsAu. W. E. Spicer, 

A. H. Sommer, and J. G. White, RCA Laborato- 

ries, Princeton, New Jersey (Received February 
5, 1959). 


Data are presented showing CsAu to be a semi- 
conductor. The optical absorption curve is typi- 
cal of a semiconductor with a threshold for direct 
optical excitation across the gap between 2.6 and 
3.3 ev and exciton peaks at about 2.6, 3.0, and 
3.3 ev. A relatively high optical absorption at 
smaller photon energies may indicate an indirect 
absorption process. Additional information has 
been obtained by measuring the spectral response 
curve of the photoemission from CsAu. The re- 
sistivity of the material is approximately 0.01 
ohm-cm over the range from 410°K to 4.2°K, in- 
dicating either degeneracy or impurity banding. 
X-ray studies show that CsAu has a CsCl struc- 
ture and a lattice constant of 4.263+0.001 A. 
Finally, some preliminary results obtained with 
other alkali metal-gold compounds are described. 
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ELECTROMAGNETIC PROPERTIES OF SUPER. 
CONDUCTORS. K. K. Gupta* and V. S. Mathur, 
Enrico Fermi Institute for Nuclear Studies and 
the Department of Physics, University of Chi- 
cago, Chicago, Illinois (Peceived January 2, 
1959). 





An integral relationship between the super - 
current and the vector potential of the electro- 
magnetic field, similar to Pippard’s phenome- 
nological equation, is derived from the gauge- 
invariant theory of Wentzel. Expressions for 
coherence distance and penetration depth have 
been calculated. 


*on leave of absence from the Tata Institute for 
Fundamental Research, Bombay, India. 


D BANDS IN THE BODY-CENTERED CUBIC 
LATTICE. Joseph Callaway,* Department of 
Mathematics, Queen Mary College, University 
of London, London, England (Received February 
2, 1959). 


A model of a crystal, consisting of positive 
point charges neutralized by a uniform distribu- 
tion of negative charge, is employed to study the 
form of the d-bands in a body-centered cubic 
lattice as a function of the lattice spacing. The 
wave functions are expressed as linear combina- 
tions of plane waves and the potential treated as 
a perturbation. It is shown that the perturbation 
series for the energy is a power series in Za, 
where Z is the atomic number and a is the lat- 
tice parameter. The leading term in the series 
is of the order (1/a)’, and the coefficients of 
successive terms in the series decrease rapidly. 
The first three terms are evaluated for the states 
of predominantly d symmetry at the center of the 
Brillouin zone, and the corner H. 


"Permanent address: Department of Physics, Uni- 
versity of Miami, Coral Gables, Florida. 


FREQUENCY SHIFTS IN HYPERFINE SPLITTING 
OF ALKALIS CAUSED BY FOREIGN GASES. 

H. Margenau, P. Fontana, and L. Klein, Yale 
University, New Haven, Connecticut (Received 
February 2, 1959). 


The difference in the dispersion force between 
an alkali atom in a particular hyperfine level be- 
longing to the ground state and a perturbing mole- 
cule is computed. These asymptotic forces, if 
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active alone, generally produce red shifts and 
suffice to account for the results obtained for 
the heavier buffer gases. Experimental data 
exhibit blue shifts for the lighter gases and there- 
fore indicate that the net frequency shifts are the 
result of exchange as well as dispersion forces. 
The former cannot be determined theoretically. 
Therefore, a simple model is devised in which 
the difference in the forces is given a positive 
trend at distances of separation smaller than d, 
while beyond d it is given the calculated form. 

The experimental data for all alkali-rare gas 
interactions can then be fitted by values of d 
which, for the different foreign gases, are of 
order of their gas-kinetic diameters. 


ABSORPTION OF GASEOUS HELIUM IN THE 
EXTREME ULTRAVIOLET. Norman N. Axelrod 
and M. Parker Givens, Institute of Optics, Uni- 
versity of Rochester, Rochester, New York 
(Received February 2, 1959). 


The absorption of gaseous helium at room tem- 
perature has been studied in the spectral region 
550 A to 150 A. The measurements of Lee and 
Weissler have been confirmed and extended. The 
data are in substantial agreement with the theo- 
retical predictions of Huang. 


NEUTRON ELASTIC AND NONELASTIC CROSS 
SECTIONS FOR BERYLLIUM. J. B. Marion, 

Los Alamos Scientific Laboratory, Los Alamos, 
New Mexico, and University of Maryland, College 
Park, Maryland, and J. S. Levin and L. Cranberg, 
Los Alamos Scientific Laboratory, Los Alamos, 
New Mexico (Received January 29, 1959). 


A time of-flight technique has been used to 
measure angular distributions for the elastic 
scattering of neutrons from beryllium in the 
energy range from 2.6 to 6.0 Mev. Neutron total 
cross sections were also measured and the fol- 
lowing nonelastic (total minus elastic) cross 
sections obtained: 2.60 Mev, 0.27+0.13 barn; 
3.50 Mev, 0.43+0.10b; 4.10 Mev, 0.5140.08b; 
5.00 Mev, 0.60+0.08b; and 6.00 Mev, 0.73 +0.07b. 
Differential cross sections were also obtained 
for the inelastically scattered neutrons leaving 
Be’ in the 2.43-Mev excited state. Integration of 
these angular distributions yielded the following 
inelastic scattering cross sections: 3.50 Mev, 
0.23+0.04 barn; 4.10 Mev, 0.25+0.04b; and 











5.00 Mev, 0.34+0.05b. A continuous spectrum 
of neutrons was observed which corresponds to 
the direct (m, 2m) reaction and/or the excitation 
of the Be® “level” at 1.7 Mev. This continuum 
was found to have an angular distribution roughly 
symmetric about 90°. The 2.43-Mev state of Be® 
was found to decay by neutron emission to the 
ground state of Be® with a branching ratio of 
12+5%. 


ELASTIC SCATTERING OF 3.7-Mev NEUTRONS 
FROM S, Fe, Co, Ni, Cu, AND Zn. M. K. 
Machwe,* D. W. Kent, Jr., and S. C. Snowdon, 
Bartol Research Foundation of the Franklin Insti- 
tute, Swarthmore, Pennsylvania (Received De- 
cember 15, 1958). 


The differential cross section for the elastic 
scattering of 3.7-Mev neutrons from sulfur, 
iron, cobalt, nickel, copper, ‘and zinc has been 
measured in a ring geometry using a modified 
Bonner-type scintillation detector. The meas- 
urements were taken over an angular range be- 
tween 10° and 160°. Angular resolution effects 
have been removed by a suitable iterative pro- 
cedure. Multiple-scattering effects have been 
removed by an approximate analytical treatment. 


* 
Now at Université de Montreal, Montreal, Canada. 


POLARIZATION OF 6- AND 7-Mev PROTONS 
ELASTICALLY SCATTERED BY NUCLEI. 
Robert E. Warner and W. Parker Alford, De- 
partment of Physics and Astronomy, University 
of Rochester, Rochester, New York (Received 
December 31, 1958). 


Polarization measurements have been made by 
first scattering from the target of interest and 
second scattering in a carbon polarimeter for 
which the mean second scattering angle is 48° 
and for which the mean value of the polarization 
is greater than +70 % for 5 Mev <E, <6 Mev. 
The polarization as a function of angle at 15° 
intervals from 45° to 135° is reported near 6 Mev 
for Al, Ti, V, Cr, Fe, Co, Ni, Cu, and Zn, and 
near 7 Mev for Al, Fe, Ni, and Cu. P(@) usually 
shows an oscillatory behavior, with the magni- 
tude of the maxima and minima never greater 
than about 40%. The polarization as a function 
of energy is reported for C at 45° from 5.7 to 
6.8 Mev, and for Fe and Cu at 60° and 120° 
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from 5.8 to 6.4 Mev. For C, P(45°) decreases 
from +92 % at 5.7 Mev to +36% at 6.8 Mev. For 
Cu, P(60°) and P(120°) vary slowly, if at all, 
with energy, but a rapid energy variation for Fe 
is observed, possibly due to resonance effects 
in compound elastic scattering. 


CROSS SECTIONS FOR ELASTIC SCATTERING 
AND REACTIONS DUE TO PROTONS ON N"*. 

S. Bashkin, R. R. Carlson, and R. A. Douglas,* 
State University of Iowa, Iowa City, Iowa (Re- 
ceived January 28, 1959). 


Absolute differential cross sections were meas- 


ured for the reactions N**(>,p)N**, N**(o, a,)C™ 
(ground state), and N“°(p, a,)C**(first excited 
state). Thin gas targets, enriched up to 98% in 
N'*, and a 6-mil CsI(T1) crystal detector were 
used. Numerous angles were studied for bom- 
barding energies between 1 Mev and 3.6 Mev. 
Spins and parities are discussed where reso- 
nances suggest the existence of excited states 
in O'*. Results are as follows: 





E p (kev) Resonant mode Spin, parity 
1028 payayy 17 
1210 paya, 5 
1640 p a, sd 
1890 ay 
1985 ay 2- 
3000 papa, 4* 
3300 p a, 2" 
3350 ay 2* 
3520 pagar, ee eg 


At the 1210-kev resonance, the incident energy 
for which the a, yield is maximum depends on 
the angle of observation. This is probably the 
result of interference and has bearing on the 
spin and parity assignment of this level. The 
anomaly at 1890 kev cannot be positively identi- 
fied as due to resonance. 


"Now at Universidade de Sa Paulo, Sao Paulo, Brasil. 


CROSS SECTION FOR ELASTIC SCATTERING 


OF PROTONS BY N“. S. Bashkin, R. R. Carlson, 


and R. A. Douglas,* State University of Iowa, 
Iowa City, Iowa (Received January 28, 1959). 


The elastic scattering cross section for protons 


on N'* was measured for bombarding energies 
from 900 kev to 4000 kev and for angles of 90°, 
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125.3°, 149.4°, and 160.9° in the center-of-mass 
system. A gas target chamber and a 6-mil thick 
CsI(T1) crystal detector were used. The cross 
section was measured to an absolute accuracy 
of +3.5%. An anomaly was observed at 3880 + 40 
kev corresponding to a state in O'* at 11.00 Mey 
previously seen in the yield of inelastic gamma 
rays. 





‘Now at Universidade de Sao Paulo, Sao Paulo, Brasil. 


TOTAL PHOTONUCLEAR ABSORPTION IN Al. 
M. V. Mihailovic, G. Pregl, G. Kernel, and M. 
Kregar, J. Stefan Institute, Ljubljana, Yugosla- 
via (Received January 15, 1959). 


A direct measurement of the total photonuclear 
absorption in Al was carried out with improved 
resolution using a Compton spectrometer for 
analysis of the bremsstrahlung and absorption 
spectrum. The maximum value for the cross 
section is found to be 100+10 mb; the integrated 
cross section is estimated to about 800 mb Mev. 


SCATTERING OF 29.2-Mev He* NUCLEI BY 
DEUTERONS. J. S. C. McKee and D. R. Sweet- 
man,* Department of Physics, University of 
Birmingham, Birmingham, England, P. V. 
March and W. T. Toner, Department of Natural 
Philosophy, University of Glasgow, Glasgow, 
Scotland, and W. M. Gibson, CERN, Geneva, 
Switzerland (Received February 3, 1959). 


Differential cross sections for He* particles 
elastically scattered from deuterium are pre- 
sented for angles from 10° to 40° (lab). Cross 
sections obtained previously for the same pro- 
cess are also shown, and there is good general 
agreement between the two sets of results. 


*Now at Atomic Weapons Research Establishment, 
Aldermaston, Berks, England. 


C** -N** MASS DIFFERENCE AND MASS EX- 
CESSES OF SOME LIGHT NUCLIDES. L. Fried- 
man, W. Henkes, * and D. Christman, Chemistry 
Department, Brookhaven National Laboratory, 
Upton, New York (Received February 5, 1959). 


Mass excesses of the nuclides H', D?, C!?, C™, 
N**, O'*, and A® were determined by measure- 
ment of 18 parent molecule ion doublets with the 
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mass synchrometer. The C’* -N** difference, 
156.44+0.29 kev, was obtained from measure- 
ments on the doublets [C’*O** -O,], [N**O"* -O, ], 


(c“H,-H,O], [C’*CH,-CO], and [N,H, -O,] and on 


other doublets which give some of those above. 

This is in excellent agreement with the value of 

156+ 1 kev obtained from determination of the 
' (¢ beta end point. The consistency of these data 
| supports the 1-kev upper limit for the neutrino 
rest mass derived from the H* -He* mass dif- 
ference and the H® beta end point, and also sup- 
ports the argument that mass synchrometer data 
are essentially free of systematic errors of the 
magnitude required to account for discrepancies 
between nuclear reaction and mass-spectro- 
graphic mass excesses of C’” and O"*. Mass ex- 
cesses of the light nuclides determined were 
slightly higher than recently published mass syn- 
chrometer data, but agreement was generally 
within 1MU- 


7. 
Present address: Robert Kochstrasse 8, Marburg/ 
Lahn, Germany. 


POSITIVE-PION PRODUCTION IN PROTON- 
PROTON COLLISIONS AT 450 Mev. Lee G. 
Pondrom,* Enrico Fermi Institute for Nuclear 
Studies, University of Chicago, Chicago, Illinois 
(Received January 28, 1959). 


The external proton beam of the Chicago cyclo- 
tron has been used to study the reaction p+p 
-1*+n+p at 450 Mev. The reaction p+p—1*t+d 
was also measured as a check on the equipment. 
A double-focusing wedge magnet was used to ob- 
serve the momentum distribution of pions at 
angles of 14°, 20.5°, and 30° in the laboratory 
system. The over-all resolution of the magnet 
system varied from 1.5% to 2.2%. A least- 
squares fit of the p+p-~-17*+d data in the center- 
of-mass system gives do/dQ « (0.38 + 0.20 +cos70) 
with o=1.47+0.12 millibarns. The spectra for 
the unbound reaction were fitted to Ss, Sp, Ps, 
and Pp spectrum shapes predicted by the phenom- 
enological theory of pion production. This fit 
gives do/dQ « (0.40 + 0.08 +cos*6) with o=1.80 
+0.7 millibarns. The results of this fit have 
been compared with the cross section for neutral 
pion production at 450 Mev. This comparison 
| Was found satisfactory for the amounts of the 
Ss, Sp, and Pp spectrum shapes found in this 
experiment, but an excessively large amount of 
the Ps shape, 0.38+0.10 millibarn, was required 











to fit the data. The results show that o,,(Pp) is 
much smaller than has been previously reported. 
An excitation function for the unbound cross sec- 
tion is predicted by this fit in terms of the pheno- 
menological model, and is found to increase too 
rapidly at high bombarding energies. 


oe 
Now at the Aeronautical Research Laboratory, 
Wright-Patterson Air Force Base, Ohio. 


PROTON COMPONENT OF THE PRIMARY 
COSMIC RADIATION. Frank B. McDonald, 
State University of Iowa, Iowa City, Iowa, and 
William R. Webber, * Imperial College, London, 
England (Received January 14, 1959). 


The proton component has been extensively 
studied at high altitudes on a series of six Sky- 
hook balloon flights at various latitudes using 
the Cerenkov-scintillation counter technique. 
The intensity of primary protons has been meas- 
ured at 4°, 41°, 53°, and 55° geomagnetic lati- 
tude. The proton differential energy spectrum 
has been measured directly in the region 250- 
750 Mev. It is observed that primary alphas 
and protons have the same form of rigidity spec- 
trum from 1 Bv to 17 Bv. A series of cutoff 
rigidities are measured in the vicinity of ) = 53° 
and 55°. The intensity and composition of fast 
splash albedo is determined at ,=4°, 53°, and 
55°. Significant time variations are observed 
between the three high-latitude flights. Differ- 
ences observed in the intensity and energy spec- 
tra of these three flights are discussed and 
strong restraints are placed on the possible mo- 
dulating mechanism. 


*Now at the Physics Department, University of Mary- 
land, College Park, Maryland. 


PRODUCTION OF ©* HYPERONS BY 990-Mev 
POSITIVE PIONS IN LIQUID HYDROGEN. A. R. 
Erwin, Jr.,* and J. K. Kopp, Brookhaven National 
Laboratory, Upton, New York, and Harvard Uni- 
versity, Cambridge, Massachusetts, and A. M. 
Shapiro, Brown University, Providence, Rhode 
Island, and Harvard University, Cambridge, 
Massachusetts (Received February 2, 1959). 


The production and subsequent decay of six =* 
hyperons was observed in the central region of a 
hydrogen bubble chamber exposed to 7* mesons 
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of kinetic energy 900 +30 Mev at the Cosmotron. 
The corresponding total production cross section, 
corrected for scanning inefficiencies, is 0.1849-09 
mb. This value, combined with weighted aver- 
ages of £° and = production cross sections in 
hydrogen obtained by Brookhaven, Columbia, and 
Berkeley groups, is used to compute the three 
triangular inequalities imposed by charge inde- 
pendence. The inequality [o(=*)]”*+[o0(Z=~)]” 
-[20(2°)] =0.1420°2 (mb)”? >0, is the only one 
that might not be satisfied. The near-zero value 
of the left-hand side indicates the difficulty to 

be encountered in demonstrating a failure of 
charge independence at this energy. 


. 


Now at the Department of Physics, University of 
Wisconsin, Madison, Wisconsin. 


MODEL OF A LINEAR HARMONIC OSCILLATOR 
IN THE GENERAL THEORY OF RELATIVITY. 
Henry Zatzkis, Mathematics Department, Newark 
College of Engineering, Newark, New Jersey 
(Received January 16, 1959). 


A model of a linear harmonic oscillator in the 
general theory of relativity is examined. It is 
based on a classical model of a point mass vi- 
brating harmonically about the center of an ideal 
liquid sphere. The motion is no longer harmonic. 
However, to a first approximation, it is still har- 
monic, but the period of vibration depends both 
on the amplitude and on the curvature of the 
space. Both quantities tend to slow down the mo- 
tion. In the limiting case of flat space the classi- 
cal frequency is again restored. 


THREE-PION CONTRIBUTION TO THE ELEC- 
TROMAGNETIC STRUCTURE OF THE NUCLEON. 
B. Bosco and V. De Alfaro, Istituto di Fisica 
dell’Universita, Torino, Italia, and Istituto 


es 


Nazionale di Fisica Nucleare, Sezione di Torino 
Italia (Received February 2, 1959). 


In this paper the three-pion contribution to the 
scalar electromagnetic form factors of the nuc- 
leon is computed by using the fixed-source meg, 
theory without rescattering corrections. The 
(y,37) interaction is taken phenomenologically as , 
a point interaction. It is shown that for those 
values of the (7,37) coupling constant compatible | 
with photoproduction experiments, the experi- 
mental charge distribution could be roughly 
fitted with a cutoff in the dispersion integrals of 
the order of 7 pion masses. , 


STRUCTURE THEOREM FOR THE PHOTON 
PROPAGATOR. Herbert M. Fried, Department 
of Physics, University of California, Los Ange- 
les, California (Received February 3, 1959). 


A simple theorem relating the structure of the 
bare and dressed photon commutators and prop- 
agators is derived, and its implication, with re- 


spect to the choice of photon gauge, is discussed. 
} 


Ky, DECAY AND LEPTONIC DECAY OF HYPER- 
ONS. B. Sakita, Department of Physics and As- 
tronomy, University of Rochester, Rochester, 
New York (Received January 19, 1959). 


The partial life time of K — » + vis calculated 
by dispersion techniques assuming that the K 
meson is pseudoscalar and that the relevant Fermi 
interactions are of the V-A type. The results are 
compared with experiments and it is concluded 
that the renormalized axial vector coupling which | 
is responsible for the leptonic decay of hyperons 
is much smaller than the usual universal Fermi 
interaction. 
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